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1. INTRODUCTION 


This document constitutes the Final Report prepared by 
ADCOM, a Teledyne Company, for NASA Electronics Research 
Center, Cambridge, Massachusetts, under contract NAS12-567. 

The report presents work performed by ADCOM during the period 
15 May 1968 to 15 February 1969. 

1. 1 Ob.iect of the Program 

The object of the program conducted during the period 
15 May 1968 to 15 February 1969 was to investigate the effects of 
multipath interference on a NAV/SAT ranging receiver for the 
case where the range modulation signal consists of a BINOR ranging 
code directly phase modulating a 1. 6 GHz carrier, and to design 
an alternate ranging modulation and receiver which is capable of 
improved (and optimized) performance. Finally, the suitability of 
the steerable null technique developed during Part I of this contract 
(and reported in the Interim-Scientific Report dated April 15, 196 8) 
to multipath rejection in the 1. 6 GHz multipath' channel was to be 
investigated. In particular, the following tasks were specified; 

Task I - Modify the multipath channel characterization for 
application to satellite-to-aircraft links operating at 1540 - 1660 MHz. 
Diffuse scattering of the multipath return and the effect of L-band 
antenna characteristics shall be included. 

Task II - Refine the performance analysis of BINOR code 
acquisition and ranging to include the effects of carrier and clock 
jitter and the presence of multipath on acquisition speed and ranging 
accuracy. 

Task III - a. Design a signal for NAV/SAT application which is 
optimized with respect to acquisition speed, ranging accuracy, 
performance in the presence of multipath, hardware complexity 
and reliability. 

b. The basic format of the signal is a carrier 
(1. 6 GHz) which is phase modulated by a subcarrier (320 kHz) 
which, in turn, is modulated by the sum or hard limited sum 
(BINOR) of range ambiguity resolving square waves. 
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c. The parameters to be optimized are: deviation 
of the subcarrier on the carrier; method of modulation and modula- 
tion index for the combined range ambiguity resolving square waves; 
ratios among square wave frequencies; and the method of combining 
the square waves. 

Task IV - Design a receiver for reception and demodulation of 
the signal derived in Task III. Perform a digital computer simula- 
tion of a phase model of the receiver reception process to verify 
the optimization procedure of Task III, to demonstrate feasibility 
of the receiver, and to establish a trade-off between hardware 
complexity and system performance. 

Task V Modify the multipath rejection technique and simulator 

constructed under Item 2 so that it represents the multipath 
channel in the 1540-1660 MHz region and employs the signal process 
derived under Task III. 

1. 2 Organization of the Report and Summary of Results 

The organization of the report follows the order of the 
task statements presented above. 

Section 2 treats the problem of characterizing the multi- 
path channel in the 1. 6 GHz NAV/SAT link. This problem has been 
conveniently divided into three parts, which are: 

1) The physics and geometry of multipath scattering at 
1. 6 GHz off the ocean in the NAV/SAT link. 

2) Characterization of the received multipath signal 
including definitions of, and bounds on, the multipath 
signal parameters. 

3) L-band Antenna characterization in order to estimate 
the minimum antenna multipath-suppression factor 
that can be obtained for grazing angles near 10°. 
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The physics and geometry of the multipath scattering in 
the NAV/SAT link are discussed and the characteristics of the 
multipath interference for both specular and diffuse scattering are 
presented as well as the conditions under which each prevails. 

These characteristics include the relative received multipath 
power, the relative specular multipath time delay and the time 
delay spread for diffuse multipath, the relative specular multipath 
carrier doppler frequency, the fading bandwidth of the received 
diffuse multipath signal, the angles of arrival of the direct-path and 
multipath signals, and the size of the scattering region for pure 
specular and pure diffuse scattering. In addition, the direct-path 
carrier doppler frequency is computed. An interesting phenomenon 
concerning the dependence of the relative specular multipath carrier 
doppler frequency on aircraft altitude fluctuations is also discussed. 
The effects of these small unavoidable altitude fluctuations is shown 
to predominate over the effects of the horizontal aircraft motion 
under certain frequently encountered conditions. This fact has 
apparently been overlooked by some investigators. A discussion 
of the validity of certain results for diffuse scattering which have 
appeared in the literature is also presented, 

The received multipath signal is expressed in a form which 
is both suitable for computing its effects on the performance of the 
NAV-SAT receiver and which contains all the essential characteris- 
tics derived from the physics and geometry considerations. 

Finally, the characteristics of aircraft -mounted L-band 
antennas is presented, Envisioned for the NAV/SAT application 
are three antennas. One on the top of the aircraft for communi- 
cating with satellites at high elevation angles, and two side-mounted 
antennas (one on each side of the aircraft) for communicating with 
satellites at low elevation angles. What is desired for the side- 
mounted antennas is a relative multipath suppression of at least 
6 dB for even the lowest grazing angles which are being considered 
(i. e. , 10°) and including a possible aircraft roll or banking angle 
of up to 1 0°. 

The analysis of the BINOR acquisition and ranging is 
refined in Section 3, to include the effects of carrier and clock 
jitter and the presence of multipath on acquisition speed and 
ranging accuracy in the proposed NAV/SAT ranging system. 
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To this end, calculations for the probability of correct 
acquisition of the BINOR code are made for three cases. In each 
subsequent case additional factors are included in the calculation. 
The three cases are: 

a) Case where effects of multipath, carrier and 
clock loop phase jitter, and video filtering 

of the BINOR code in the receiver are neglected. 

b) Case where effects of multipath are neglected but 
effects of carrier and clock loop phase jitter, and 
video filtering of the BINOR code in the receiver 
are included. 

c) Case where multipath effects as well as the other 
effects are included. 

The third case is, of course, the most complex, and is 
best understood in terms of the modifications required in the 
results of the first two cases. 

From the analyses conducted in Section 3, the following 
conclusions are drawn: 

1) The presence of a small multipath interference will 
not cause a drastic reduction in the total probability of correct 
BINOR acquisition if 

a) at least a 6 dB aircraft antenna discrimination 
against the multipath is maintained 

b) aircraft speeds are on the order of 600 mph or 
more so that the frequency spreading of the diffuse 
multipath component is on the order of 50 Hz or 
more 

c) the grazing angle exceeds 10° and the rms displace- 
ment of the sea surface exceeds 10 cm, so that the 
average relative specular multipath power is no 
larger than approximately 0. 06 5. 
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For the vast majority of the time the grazing angle 
exceeds 2 0° and the aircraft antenna multipath discrimination - 
exceeds 6 dB. The specular multipath component is then negligible, 
and the frequency spreading of the diffuse multipath component is 
more on the order of 100 to 2 00 Hz. The effects of the multipath 
in this case becomes almost entirely negligible during BINOR acqui- 
sition. 

For grazing angles close to 10°, specular multipath may 
become a problem for calm seas at 1. 6 GHz. However, the prob- 
ability of correct BINOR acquisition is still on the order of 9 8% 
on the average for 4 = 5. Analyses of additional worst-case possi- 
bilities that may occur in the presence of specular multipath indicate 
that it is very unlikely that the probability of correct BINOR acqui- 
sition (with 4 = 4) will fall below 95% at the worst-case SNR used in 
all of the calculations. v 

Ranging accuracy is shown to be a function of the rms phase 
jitter in the clock loop. This phase jitter is shown to result from 
the presence of Gaussian noise and multipath components in the 
loop noise bandwidth. During the ranging measurement, in the 
present system, the multipath becomes an important consideration. 
This is so because the BINOR is removed and only square-wave 
clock is transmitted. This increases the signal and multipath 
powers in the clock loop to the point where the multipath, rather 
than the noise, becomes the limiting factor. This appears to be 
especially true when specular multipath is present. However, the 
differential doppler offset between the direct-path and multipath 
carrier frequencies may be sufficiently high compared to the clock 
loop noise bandwidth to cause considerable suppression of the 
specular multipath clock component. This would substantially 
reduce the limiting effects of the specular multipath component on 
ranging accuracy. The diffuse multipath component causes an 
increase in rms range error of approximately 50% over that caused 
by front-end noise alone. The effect of the diffuse component is 
reduced only when the fading bandwidth is much smaller than the 
differential doppler which is unlikely in the present application. 
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In Section 4, the optimum modulation and receiver design 
problem is presented. The discussion begins with an investigation 
of the possible modulation methods by which the ambiguity resolving 
square waves may be applied to the 320 kHz subcarrier. Low devia- 
tion PM is selected. Then the optimum method for encoding these 
square waves is considered. The BINOR technique (hard limited 
sum of square waves) is rejected in favor of a simple sum of square 
waves. Optimum deviations for the square waves on the subcarrier 
and the subcarrier on to the carrier are selected. Then, a des- 
cription of the ambiguity resolving procedure is presented. The 
discussion concludes with a hardware design exercise. This effort 
results in a detailed block diagram of a range data demodulation 
and measurement system appropriate to the optimized ranging 
modulation. 

The ranging modulation which has been designed for 
NAV/SAT applications provides the following characteristics: 

1) No sequencing is required at the satellite transmitter. 
All modulation functions are carried simultaneously. 

2) Signal format provides for convenient generation at 
the satellite. 

3) Ranging accuracy on the order of 10 ft. 

4) 99. 9% prob. of correct amb. resolution in 1 sec. 

5) Carrier suppression of only 2 dB by modulation. 

6) Carrier and clock jitter have negligible effect on 
ambiguity resolution process. 

7) Improved performance in the presence of diffuse 
multipath because all data is averaged over the 
entire acquisition interval. 

8) Receiver demodulator complexity comparable to 
that required for the BINOR (when range measure- 
ment phase meter and system reference generator 
are included in the comparison). 
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9) Receiver hardware lends itself readily to medium 
or large-scale integration. 

and finally, 

10) The signal format can provide simultaneous data 

transmission capability employing straightforward 
communication transmission and reception concepts 
and, implementations. 

These characteristics are desirable for the NAV/SA-T 
application. It is therefore recommended that the modulation 
technique be given further consideration. In particular, modula- 
tion deviations may be readjusted to reduce ranging accuracy 
(which is 3 times better than required) in trade for shorter acquisi- 
tion time and/or higher probability of acquisition and/or higher 
carrier component power. It is noted, however, that these altera- 
tions are not necessary for successful implementation of the mod- 
ulation technique. Further analysis of the demodulation and 
measurement technique will yield improvements in implementation 
beyond those presented in the section. Again, it is noted that the 
complexity of the implementation presented is on a par with other 
proposed systems. It is further recommended that the communica- 
tion capabilities of the signal format be investigated in more detail 
to determine the optimum modulation method and maximum data 
rate capability. 

A digital computer simulation of the receiver digital- 
divide -down carrier acqiiisition aiding PLL is presented in 
Section 5. Analysis of the signal format derived in Task III and 
the corresponding demodulator permitted optimization of the sig- 
nal modulation parameters to be achieved without the aid of 
computer simulation. However, it became apparent that the time 
required for carrier acquisition in the present NAV/SAT receiver 
configuration represented a serious drawback on overall acquisition 
speed. Consequently, a digital-divide-down carrier acquisition 
aiding PLL was proposed as a technique for reducing this contribu- 
tion to the overall acquisition time. The difficulties involved in 
analyzing the operation of this acquisition aiding PLL at low input 
SNR led naturally to the use of digital computer simulation as an 
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aid in the design of this loop. The results of this analysis indicated 
that a first order loop performs better in the acquisition mode than 
a second order Loop, Criteria for the selection of the digital-divide - 
down ratio and the loop gain were developed in terms of the expected 
maximum carrier frequency offset and the expected minimum SNR. 
Acquisition times on the order of 10 milliseconds or less were 
consistently obtained in the computer simulation for input SNR down 
to -10 dB in a 20 kHz IF bandpass filter. It was also found that 
the loop would acquire for input SNR's down to -3 0 dB. 

Finally, in Section 6, the feasibility of steerable null con- 
trol for the rejection of diffuse multipath is investigated. A simula- 
tion experiment was also performed to investigate the ability of the 
steerable null control system to reduce the effects of coherently 
fading diffuse multipath interference. ADCOM's conclusions con- 
cerning the effectiveness of the steerable null control system for 
combatting multipath interference at 1. 6 GHz in a NAV/SAT applica- 
tions are summarized as follows: 

1) Since diffuse multipath interference arises from scatter- 
ing off a much larger area of the sea surface than for the case of 
specular scattering, the multipath interference in the diffuse case 
arrives at the aircraft from a range of different directions centered 
on the specular direction. Consequently, placing a sharp null in 

the direction of the specular point will not eliminate all of the multi - 
path interference. In addition, slight variations of antenna gains 
with directions of arrival of the diffuse multipath components will 
preclude the possibility of achieving the required amplitude gain 
balancing necessary for achieving a true null. At best, only the 
coherently fading diffuse multipath component can be nulled, with 
the incoherently fading diffuse multipath interference being simply 
reduced somewhat. 

2) The presence of incoherently fading diffuse multipath 
interference interferes with proper operation of the steerable null 
control system and may prevent the null from being placed on the 
direction of arrival of the coherently fading diffuse multipath 
component. 
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3) Steerable null control at 1. 6 GHz does not work as well 
at low grazing angles as at intermediate grazing angles (i. e. , around 
30°), even for the case of specular scattering. For the 1. 6 GHz 
NAV/SAT channel grazing angles as low as 10° are being considered. 

4) Steerable null control may not be needed at 1. 6 GHz, since 
antennas can be built at these frequencies which may adequately dis- 
criminate against the multipath if the angular separation between the 
directions of arrival of the direct path and multipath signals at the 
aircraft is greater than about 2 0°. 

A modification of the steerable null simulator (constructed 
previously under this contract) was made to simulate the randomly 
fluctuating carrier phase of a coherently fading diffuse multipath com- 
ponent. The results indicated that the steerable null system could work 
well under such circumstances in that the rapidly fluctuating multipath 
carrier phase prevented the phase of the multipath carrier from dwelling 
too long in a region where operation of the steerable null control system 
became ineffective. 

Since the steerable null system would not be used for the 
situation where incoherently fading diffuse multipath interference was 
dominant, the simulator was not modified to simulate this kind of mul- 
tipath interference. Such an effort could not be justified in light of the 
theoretical considerations itemized above. These considerations are 
based on the results derived and presented in Section 2 of this report 
concerning the characteristics of multipath at 1. 6 GHz. 
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2. MULTIPATH CHANNEL CHARACTERIZATION 
2. 1 Summary 

The problem of characterizing the multipath channel in the 1. 6 GHz 
NAV/SAT link is conveniently divided into three parts. These are: 

1) The physics and geometry of multipath scattering at 1. 6 GHz 
off the ocean in the NAV/SAT link. 

2) Characterization of the received multipath signal including 
definitions of, and bounds on, the multipath signal parameters. 

3) L-band Antenna characterization in order to estimate the 
minimum antenna multipath-suppression factor that can be 
obtained for grazing angles near 10°. 

The physics and geometry of the multipath scattering in the NAV/ 
SAT link are discussed in Section 2. 2. The characteristics of the multipath 
interference for both specular and diffuse scattering are presented as well 
as the conditions under which each prevails. These characteristics include 
the relative received multipath power, the relative specular multipath time 
delay and the time delay spread for diffuse multipath, the relative specular 
multipath carrier doppler frequency, the fading bandwidth of the received 
diffuse multipath signal, the angles of arrival of the direct-path and multi- 
path signals, and the size of the scattering region for pure specular and pure 
diffuse scattering. In addition, the direct-path carrier doppler frequency is 
computed. An interesting phenomenon concerning the dependence of the 
relative specular multipath carrier doppler frequency on aircraft altitude 
fluctuations is also discussed. The effects of these small unavoidable altitude 
fluctuations is shown to predominate over the effects of the horizontal aircraft 
motion under certain frequently encountered conditions. This fact has ap- 
parently been overlooked by some investigators. A discussion of the validity 
of certain results for diffuse scattering which have appeared in the literature 
is also presented. 

The characterization of the received multipath signal is discussed 
in Section 2. 3. The received multipath signal is expressed in a form, which is 
both suitable for computing its effects on the performance of the NAV/SAT 
receiver and which contains all the essential characteristics derived from 
the physics and geometry considerations. 

Finally, the characteristics of aircraft-mounted L-band antennas 
is presented. Envisioned for the NAV/SAT application are three antennas. 
One on the top of the aircraft for communicating with satellites at high 
elevation angles, and two side-mounted antennas (one on each side of the 
aircraft) for communicating with satellites at low elevation angles. What 
is desired for the side -mounted antennas is a relative multipath suppression 
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of at least 6 dB for even the lowest grazing angles which are being con- 
sidered, (i. e. , 10°) and including a possible aircraft roll or banking angle 
of up to 10°. 

2. 2 Physics and Geometry of Multipath 
Scattering in the NAV/SAT Link 

2. 2. 1 Scattering Geometry of the NAV/SAT Link 

The geometry of the NAV/SAT configuration is shown in Fig. 1. 

For the case of specular scattering there is a single multipath ray as 
indicated. The ray which reaches the airplane P from the satellite S is 
the one which scatters off the earth at point M. At point M, the angle of 
incidence equals the angle of reflection. This point is referred to as the 
"specular point". The size of the "specular point" is equal essentially to the 
area of the first Fresnel zone. The concept of the first Fresnel zone is 
useful for situations in which it is small in size compared to the distances from 
it to the observer and from it to the source, and for which the source illumina- 
tes the first Fresnel zone uniformly. This is the case for the NAV/SAT 
application. 

When the scattering off the ocean is diffuse, the airplane P receives 
multipath radiation from the satellite S which has been scattered off the ocean 
from an area centered approximately at the specular point M and which is con- 
siderably larger than the first Fresnel zone. The sizes of the first Fresnel 
zone and the diffuse scattering region will be estimated, for grazing angles 
between 10° and 20°, in Section 2 , 2 , 2 . 5 . 

From the geometry of Fig. 1 several important parameters of the 
multipath channel may be defined and calculated approximately as functions 
of the aircraft altitude, h, and the angle, a, subtended between the aircraft 
and the satellite as measured from the earth center. This angle a may in turn 
be determined from the latitude and longitude separations between the aircraft 
and satellite by the spherical trigonometric relationship 

cos a = cos A a, cos Aa, , (1) 

long. lat. 


This is further illustrated by Fig. 2. 

The multipath channel parameters of interest here are: 

1) The relative specular multipath time delay, t , given by 

ms 


T 

ms 


1 _ 

c 


(R 

m 



( 2 ) 


where c is the velocity of light, and R and are the path 
lengths of the multipath and the direct path respectively. 
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Fig. 1 AC-Satellite 



Note. Not to scale Correct scale is 
R e / R s = 0.15 and h/R e ~ 0.005 


Multipath Geometry 





2) The grazing angle "Y. 

3) The angles of arrival of the direct-path and multipath signals at 
the aircraft <£3 and c£ m as measured with respect to the horizontal 
at the aircraft. The angle <£<3 is also referred to as the satellite 
elevation angle. The complement of this angle, tt / 2 - , is re- 
ferred to as the aspect angle. The directivity angle, 6 , is 
defined as the sum (j^+ 4 ^, and is the angular separation between 
the directions of arrival of the desired direct-path signal and 

the interfering multipath signal. 

4) The relative multipath and direct-path carrier doppler frequencies. 

Of particular interest is the range of grazing angles between 10° and 
20°. The minimum expected grazing angle is 10°, and for grazing angles above 
20 ° the multipath problem is not as severe as for angles less than about 20 °. 

This is so for several reasons based on characteristics of the aircraft-to- 
satellite link and the navigation system implementation. These are, in order 
of importance : 

1) The L-band antenna directivity can significantly suppress the 
multipath signal as a result of the larger directivity angles. 

2) Specular multipath interference, which may present a problem at 
lower grazing angles, is entirely negligible for grazing angles 
above about 20° at 1 . 6 GHz and for normal ocean conditions. 

3) The fading bandwidth, which is a measure of the spectral spread- 
ing of the diffuse multipath signal, is larger at higher grazing 
angles. A large fading bandwidth may reduce, the effects of the 
multipath in the narrow bandwidths of the carrier and clock PLL's 
and in the integrate and dump BINOR code acquisition filters. 

This is so when it corresponds to a frequency spreading of the 
multipath power at these particular frequencies into a bandwidth 
which is larger than the corresponding filter noise bandwidth. 

For the grazing angles in this range, and for aircraft altitudes less 
than about 20 miles, it is a good approximation to assume that the direct -path 
ray S-P and the multipath ray segment S-M are parallel as indicated in Fig. 3, 
and that effects of the earth’s curvature between the scattering point M and the 
point on the earth's surface directly under the aircraft are negligible. In this 
case, the angles V, <£^, and 4 > m are all equal, and given by 
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cot <t> 


cot "V 


(3) 


cot 4>, 
a 


R sin a 
s 


R cos a - (R + h) 
s e 


where R g and R e are the synchronous satellite orbit radius and the earth radius 
respectively, and h is the aircraft altitude. The angle a was defined above and 
in Figs. 1 and 2. Dividing numerator and denominator of Eq. (3) by R g per- 
mits cot y to be written in a normalized form as 


whe re 


cot V 


s in a 


(cos a - cos a ) 
o 



(4) 


(5) 


if h/R e « 1. with R e = 3, 960 miles, R s = 26, 260 miles (for synchronous 
equatorial orbit), and with = 20 miles, 



T. 

0. 1508 and — < 0. 00505 

e 


( 6 ) 


In terms of the grazing angle, y, the relative specular multipath time 

delay, t , is 
ms 


2h . v 2h 
t = — sm y = — 
ms c c 


V 


1 + cot "Y 


( 7 ) 


A curve of sin Y as a function of a is plotted* in Fig. 4 for V between 10° and 
30°. 


The differential carrier doppler, f , is another important para- 

III s 

meter of the multipath signal, and is given by 


* In addition to the approximations of parallel rays and locally flat earth used 
to derive Eqs . (3) and (7), the approximation indicated in Eq. (4) has also 
been made . 


17 




R-5419 


dh 

dt 
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ms 


da 

dt 
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ms 


dT 
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ms 

dt 
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ms 

9h 
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where f c is the carrier frequency (1.6 GHz). From Eqs. (4) and (7), this is 
approximately given by* t 


f 

ms 




cos Y (cot a + cot Y) v 

a 


(9) 


where and v a are the vertical and horizontal (in plane) aircraft velocity 
components given by 


v 


h 


dh 

dt 


and v 

a 


1 da 

(R + h) dt 
e 


1 da 

R dT 
e 


( 10 ) 


The horizontal component of velocity which is perpendicular to the plane con- 
taining the aircraft P, satellite S, and earth center E does not contribute to the 
differential doppler frequency. 

An interesting phenomenon concering the differential doppler, which 
has apparently been overlooked by some authors, is the relative effect that a 
small unavoidable vertical component of aircraft velocity can have on the value 
of f^ even for nominally "constant altitude" flight. This is particularly 
true at lower aircraft altitudes and for intermediate grazing angles where the 
ratio (h/R e ) is small and cot Y is not too large, so that the factor multiplying 
v a in Eq. (9) is much less than 1. To illustrate this point further, the ratio of 
v a to vfr required for the effects of vertical and horizontal (v ) aircraft veloci- 
ties on f ms to be equal, (v a /v] 1 ) 0 , is plotted as a function of a for aircraft 
altitudes of 5, 10, 15, and 20 miles in Fig. 5. From Eq. (9), the ratio (Vq/vj^) 
is, of course, 




1 


cos Y (cot a 


+ cot Y) 


(ID 


* In addition to the approximations of parallel rays and locally flat earth used 
to derive Eqs. (3) and (7), the approximation indicated in Eq. (4) has also 
been made. 

t Added in proof: See Addendum for a revised derivation of f ms - 
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From the curves in Fig. 5, it is seen that this ratio will always exceed ap- 
proximately 32 if the grazing is always greater than 10° (a 572°) and the 
aircraft altitude is always less than 20 miles. For aircraft altitudes less 
than 15 miles and for grazing angles greater than 20° (a ^64°), a much more 
frequently encountered condition, this ratio will be greater than 80. What this 
means is that for aircraft flying at 600 mph and having a v Q component of 
200 mph or less, the vertical aircraft velocity would have to be much less 
than about 6 mph for grazing angles of 10° and aircraft altitudes of 20 miles, 
and much less than about 2. 5 mph for grazing angles in excess of 20° and 
aircraft altitudes less than 15 miles, in order for the relative multipath car- 
rier doppler frequency, ^rns » to be determined primarily by the v a contri- 
bution. Consequently, since aircraft vertical velocities cannot be held to 
less than plus or minus several mph, the random or uncontrollable vertical 
aircraft motion will have a pronounced effect on f ms . Therefore a minimum 
differential or relative multipath doppler offset cannot be counted on, but only 
a range of values for f ms may be estimated, with f ms likely to be anywhere 
in that range with approximately equal probability. The range, of course, 
corresponds to the uncertainty in aircraft vertical velocity, and is centered 
on the more deterministic contribution to f from v Q . 

Finally, the direct-path carrier doppler may be determined iirom 
the geometry of Fig. 1. This is given by 


dr 


c dt 


c C dt 


( 12 ) 


where is the direct-path distance and is the direct-path transmission 

time delay. R is given by (see Fig. 1) 
d 


R , = 


R + (R + h) - 2 R (R + h) cos a 
s e s e 


1/2 


(13) 


so that 
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(15) 
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or f is 
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- f 
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R 

s 

~ (cos a - cos a ) -v + cot Y v 
R , oh a 
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(16) 


It is easily seen that for low grazing angles, Y, and for nominally level flight, 
the contribution of v q clearly dominates that of v^. This is so because, for 
aircraft flying at 600 mph or more, the aircraft will more than likely have a 
component v Q on the order of 50 mph or more (compared to several mph for 
vh)» and cot ^ will be greater than 1. Consequently, except for cases where 
v Q is very small, or where the synchronous satellite is almost directly over- 
head (a = 0), fj may be written as 


f 


d 



i R 
1 s 


sin a v 

a 


(17) 


At the low grazing angles (10° -20°) where the multipath interference will be 
most severe, the direct-path doppler frequency, f^, will be much larger than 
the relative multipath doppler frequency f 

ms 

2. 2. 2 Scattering Physics of the NAV/SAT Link 

The physics of the scattering, in conjunction with the geometry of 
the NAV/SAT link, yields information on 

1) Whether, and under what conditions, the scattering is essentially 
specular or diffuse, and the characteristics of the scattered 
signal in each case. 

2) The received average multipath power. 

3) The fading bandwidth or spectral spreading of the scattered 
multipath signal as received by the moving aircraft. 

4) The spread of time delays of the received multipath signal for 
the case of diffuse scattering. 

5) The size of the first Fresnel zone for the case of specular 
scattering. 

6) The size of the scattering region for the case of diffuse scattering. 
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2.2.2. 1 Character of the Scattered Signal 


The character of the received multipath interference depends on 
whether the surface of the earth at the scattering point M may be considered 
"rough" or "smooth" for the particular wavelength, X, of the incident radia- 
tion, the grazing angle, Y, and the rms surface height variation, <r w . The 
Rayleigh criterion is often used to differentiate between the two extremes of 
"rough" and "smooth". According to the Rayleigh criterion the surface is 
considered "smooth" if 


cr sin Y 
w 


X 



(18) 


and "rough" otherwise. Since most of the earth's surface is sea water, 
and since reflections off the sea surface represent a more serious problem 
than reflections off land, only reflections off the sea will be considered, and 
o- will be taken to represent the rms sea surface height fluctuation. The sea 
surface heights may be considered to be randomly distributed with an appro- 
ximately Gaussian distribution. 


Applying the Rayleigh criterion to the 1.6 GHz frequency, it may be 
shown that the scattering is primarily diffuse for all sea states rougher than 
"slight" at grazing angles in excess of 10°, but that specular scattering 
becomes important for sea states calmer than "smooth" at grazing angles 
below 20°. In particular, the range of grazing angles, Y, at the transition 
point between "rough" and "smooth" scattering surface conditions, as deter- 
mined by Eq. (12), for various sea states are listed in Table 1 below. 

Table 1 


Grazing Angle, Y, for which sin V = — X /cr (X = 18. 75 cm) 

o w 



Approximate 




Sea State 

range of o’ 

00 

range of Y 


Comment 

smooth 

5-10 cm 

28° - 1 3. 5° 


specular scattering becomes 
important below about 20° 

slight 

10-20 cm 

13. 5 ° -6. 7° 


specular scattering is only 
important near 10° 

moderate 

20-30 cm 

Below 10° 


specular scattering is 

rough 

30-50 cm 

Below 10° 


entirely negligible 
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For very "smooth" scattering surfaces only a specular multipath 
return is present. The amplitude and phase of the received multipath signal 
are determined essentially by the electromagnetic properties of the scattering 
surface (i.e. by the Fresnel reflection coefficients for the particular polari- 
zation of the incident radiation and angle of incidence) and the geometry of 
the link, and vary in a rather deterministic manner with the motion of the 
aircraft. As the surface becomes "rougher", either as a result of increas- 
ing grazing angle and/or rougher sea conditions, several things begin to 
happen. First, amplitude and phase fluctuations on the received multipath 
signal appear, as if cophasal and quadrature components having zero -mean 
Gaussian-distributed random amplitudes with unequal variances ’ are 
simply added to the steady deterministic specular multipath component. 
Second, the power in the steady deterministic specular multipath component 
begins to diminish with increasing surface "roughness" as a result of a 
decrease in the mean-square scattering coefficient { | P s | ^ > which is 

given by^ 


< 



e 




sin V 


2 


( 19 ) 


Further increases in surface "roughness" cause further diminish- 
ing of the power in the steady deterministic component of the received 
specular multipath component and cause the variances of the fluctuating 
cophasal and quadrature components, mentioned above, to increase and 
approach each other in value. 

Finally, as the transition from "smooth" to "rough" is completed, 

the power in the steady specular component becomes negligibly small, and 

the scattering is essentially pure diffuse. For this case radiation is received 

from a much larger area of the scattering surface centered around the 

specular point, and multipath components having relative multipath delays 

in excess of t begin to appear, 
ms 

2. 2. 2. 2 Computation of the Relative Received Multipath Power 

The relative received multipath power for the case of pure specular 

scattering is denoted by a 2 , and is given by 

ms 


2 

a 

ms 




r r. 

s A 


( 20 ) 
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where R | is the magnitude of the reflection coefficient of a very smooth 
plane sea; D is the divergence factor to account for the sphericity of the 
earth's surface; { ! ^ } is the mean square scattering coefficient which 

takes account of the slight random roughness of the scattering surface; 


2 2 2 
l R d’ 1’ 2 


is a factor which takes the relative 1/r intensity losses 


into account if necessary; and r s and Fa are the relative antenna gain factors 
of the synchronous satellite and the aircraft respectively, and includes the 
effects of polarization discrimination if any. 


The magnitude of the reflection coefficient, 1 R 0 | , for sea water 
may be obtained for horizontal and vertical polarization from Fig. 11. 2 of 
Beckman and Spizzi chino? which plots | R Q | as a function of the grazing 
angle, \, for various wavelengths, X. Special care must be taken in the 
case where circular or elliptical polarization is used in order to take cross - 
polarization and other effects into account. ® 

9 

The divergence factor D is given by 


D = 


3 R 2 <R, + R 2 ) 

-1/2 

' 2R 2 (R 1+ r/ 

+ R (R. + 2 R ) sin V 
el 2 J 


R e (R 1 + 2 V J 


- 1/2 


(2i ; 


where R g is the earth radius, R 2 is the distance from the scattering point M 
to the airplane P, and (R-^ + R^) is the distance from the satellite S to the 
scattering point M. (See Fig. 1. ) For the NAV/SAT link D = 1. 


The mean square scattering coefficient / I P, 


> 


was given as a 

function of wavelength X, rms sea surface displacement <7 , and the sine of 
the grazing angle V in Eq. (13). For the ratio (tr w sin ^)/\ equal to 1/8, 
which is the Rayleigh criterion dividing line between "smooth" and "rough" 
scattering surfaces { | I ^ ) is e q ua -l to 


<KI 2 > - 


- ( Tr / 2 ) * 


s 0. 085 


( 22 ) 


which is fairly small . 
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is simply given by 


The function F 


'222 

R d ,Rl,R 


F 





< r i 


( 23 ) 


and reflects the fact that the 1/r intensity losses along the two paths are 
slightly different. Since R^ = Rj + 2R£ (i. e. , R^ + 2 R 2 - R^ « R^), this 
factor is essentially equal to 1 for the NAV/SAT link. 


The relative satellite antenna gain factor, r s , is essentially equal 
to 1 since the satellite must provide uniform coverage of the earth's surface 
and the angular separation between the direct-path and the multipath is very 
small at the satellite. The relative aircraft antenna gain factor, r^, can be 
less than 1 as a result of the large angle between the directions of arrival 
of the direct-path and multipath which is in excess of 20°, and as a possible 
consequence of polarization discrimination. 

2 

In estimating a reasonable worst-case value of a the following 
,. , ms 

assumptions are made: 

= 1 

worst-case 

■2) (D 2 ) = 1 

worst-case 



3) 



worst-cas e 


e 


-(0. 37-n) 2 


0. 2583 


Corresponding to "Y =10° and a =10 cm. 

w 


4) 





worst-cas e 


1 


5) 


r = 1 

s 


6) (r.) , = 0. 25 (-6 db) 

A worst-case 


For these worst-case values (a ) 

ms worst-case 


0. 06457. 
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For the case of pure diffuse scattering, the average relative re- 
ceived multipath power is denoted by 20 -.^, and may be computed approxi- 
mately ^ 0 under certain conditions. These conditions are: 

1) The undulations of the scattering surface may be described by 

a two-dimensional Gaussian distribution for which the rms 

surface displacement is cr and the correlation length is L. 

w 

2 ) The ratio 2cr / L is small (less than 1 / 1 0 ) 

w 

3) The ratio (cr siny)/\ is greater than 1/8. 

w 

4) The autocorrelation function of the surface fluctuations is an 
analytic function. 

2 

The first condition is a fairly reasonable one in that measurements 
indicate that it is approximately satisfied in practice, and it also simplifies 
the mathematics. The second condition is equivalent to requiring that the 
slopes of the waves be small. This requirement is essential for the Kirchhoff- 
Huygen approximation to be valid, ^ and for shadowing and multiple scattering 
effects to be unimportant. In addition, the rms slope must be sufficiently 
less than the grazing angle, y, which would pertain for a perfectly smooth 
scattering surface in order that the reflection coefficient, R Q , used in the 
determination of the average received relative diffuse multipath power may 
be evaluated at the grazing angle y. In this regard it is also important that 
I R 0 I not vary drastically with V for those values of y for which the approxi- 
mate computation is made. In particular, the approximate computation will 
predict incorrect results for vartical polarization near the Brewster angle. ^ 
The third condition is the Rayleigh criterion condition for pure diffuse 
scattering. 

If these conditions are well fulfilled, most of the received diffusely 
scattered multipath power will come from the region surrounding the specular 
point, and the method of steepest descent may be used in performing the 
integration required in summing the powers contributed from each elementary 
area of the scattering surface. The result for 2c ^ obtained by Durrani and 
Staras, iU with appropriate modifications to take the characteristics of the 
NAV/SAT link into account*, is 


* The computation by Durrani and Staras was for an orbiting satellite having 
an altitude in excess of 100 miles and for grazing angles which could extend 
down below 10°. We also include the aircraft relative antenna gain factor 
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( 24 ) 


2cr 2 = |R | 2 D 2 r Q (h, V ) 
m 'o' A 

2 2 

where j R 0 | . D , and are defined as before, and Q (h,Y) is a function 
essentially of the aircraft altitude and the grazing angle {although Durrani 
and Staras define it in terms of the orbit height, the aspect angle, and the 
angle between the aircraft and specular point as measured from the earth's 
center). The function Q (h, Y) is approximately equal to 1 except for large 
altitudes and small grazing angles. For the NAV/SAT application it may be 
approximated by 1. Again is approximately equal to 1. The worst-case 
value of | R q | ^ for both horizontal and vertical polarization is approximately 
-2 db. 1-3 The aircraft antenna relative gain factor, r^> is again taken as 
0. 25 (-6 db), so that 2cr ^ is, in the worst-case. 


2cr - 0. 16 

m / worst-case " 


2. 2. 2. 3 Computation of the Fading Bandwidth 


(25) 


A non zero fading bandwidth arises because of relative motions of 
the transmitter and receiver, and/or the scattering surface. If the trans- 
mitter and receiver were stationary with respect to each other and the 
scattering surface, and the only motion was that of the ocean waves, the 
fading bandwidth would be relatively small. However, the high aircraft 
velocities result in very much larger fading bandwidths. 


The fading bandwidth may be determined by computing the auto- 
correlation function of the received diffusely scattered multipath signal for 
the case where carrier only is transmitted. Since the decorrelation time, 
t q , is much smaller than the periods of ocean waves, the sea may be con- 
sidered as static in the computation. This autocorrelation function has been 
computed by Durrani and Staras 10 for aircraft motion which is coplanar 
and perpendicular to the plane of the aircraft, satellite, and earth center. ^ 


t Small vertical aircraft velocities do not contribute significantly to the 
fading bandwidth (i. e. to further spectral spreading of the received 
diffuse multipath interference) if the aircraft velocity is on the order 
of 600 mph or more and is in level flight, and for grazing angles in 
excess of 10°. 
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In both of these cases the autocorrelation function is found to have a 
Gaussian shape exp (- ( T / T 0 )^)' For small aircraft altitudes (i. e. less 
than 20 miles) and for grazing angles in excess of 10°, the correlation 
times for these two cases are essentially equal, and given by 

-1 

(26) 

where v/c is the ratio of the aircraft velocity to the velocity of light, and the 
other quantities have been defined previously. 

The relative power density spectrum of the received multipath for 
the case where only a carrier is transmitted is 
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(28) 


and w is measured from the received multipath carrier frequency (including 

specular doppler offset). For v = 600 mph and (2cr /L) = 1/10, B is 

w m 

B = 286 sin "Y Hertz (29) 

m 


which for \ = 10°, yields B =50 Hz. 

m 

2. 2. 2. 4 Computations of the Sizes of the Scattering Region for Specular 
and Diffuse Scattering and the Multipath Time Delay Spread 

For pure specular scattering, the primary contribution to the 
received multipath signal comes from the first Fresnel zone. This is an 
area of the scattering surface centered at the specular point such that the 
difference between the path distances of aircraft received radiation scattered 
off any point in this area and of aircraft received radiation scattered off the 
specular point is less than a half-wavelength of the radiation. A perfectly 
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smooth and flat area which is illuminated uniformly with radiation is assumed. 
For the NAV/SAT application, this area is approximately elliptical with 
semimajor and semiminor axes which are denoted by AX g and Ay g respectively 
as shown in Fig. 6, and which may be determined from the equations 
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(30) 


(31) 


(32) 


(33) 


Expanding the square roots in Taylor series and retaining terms up to 

Ax ^ anc i Ay^, Ax and Ay are approximately 
s ss s 


Ax 

s 


1 

sin Y 



(34) 



(35) 


The maximum value of Ax^ and A y s occur for the maximum aircraft altitude 
(20 miles) and the minimum grazing angle (Y =10°) and are approximately 
Ay s =; 186 meters and Ax^ 1070 meters. 


The size of the scattering region for pure diffuse scattering may be 
obtained as follows: The distances Ax^ and Ay^ (see Fig. 7) for which the 
diffuse multipath power reflected per unit area drops to 1/e of the diffuse 
multipath power reflected per unit area at the specular point M are computed 
and taken as representing size of the scattering region. The variation in 
diffuse multipath power reflected per unit area over the scattering region is 
a result primarily of the variation of the term^ 


tan A 



( 36 ) 
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Point 



To Satellite 


or 


Fig. 6 Configuration for Computing the Semimajor and Semimin 
Axes of the First Fresnel Zone Ellipse 
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Point 



To Satellite 


Fig. 7 Configuration for Computing the Time Delay Spread of 
the Diffuse Multipath Scattering 
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where here 3 is the angle made by the bisector of the incident and scattered 
rays with the z-axis. The quantity tan^p is a function of the location of the 
unit scattering area with respect to the specular point; the grazing angle, y , 
and the aircraft altitude, h. If Ax and Ay are the position coordinates with 
respect to the specular point of the unit scattering area under consideration, 
tan^p is approximately given by 

2 1 
tan 3 * 

4h 

to second order in Ax and Ay. Then Ax and Ay are 

d d 


. 4 
sm y 


Ax + Ay 


(37) 
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This may be rewritten as 
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(41) 


where Ax s and Ay s are the semimajor and semiminor axes of the first 
Fresnel zone ellipse. Note that the factor (sin y) in Eq. (41) means that the 
diffuse scattering region will be much more elongated in the x-direction than 
is the first Fresnel zone. Considering the locus of points Ax, Ay for which 
the diffuse multipath power reflected per unit area drops to 1/e of what it is 
at the specular point, it is noted that these points correspond to different 
relative multipath delays. Since the maximum spreading in relative multi - 
path delay for this locus occurs for Ay = 0, Ax = Ax^, the relative multipath 
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delay difference (T md - T ms ) shall be taken to be a measure of the multi - 

path delay spread. Here T m( j is, of course, the relative multipath delay for 

the scattering point at Ay = 0, Ax = Ax^. The relative multipath time delay 

spread, (t - t ), is then 
md ms 
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md 
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ms 
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where 
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If Ax^ is sufficiently smaller than (h/sin V), d may be expanded to second 
order in Ax^, and ( T m( £ - T ms ( i s fcmnd to be approximately given by 
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To check the validity of this approximation consider the require- 
ment that Ax^ be sufficiently less than (h/sin Y). 


Ax = 
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sin Y 
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( 46 ) " 
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For Y i 20° and | ~)~ Ys - 


2 

sin V 


2 g 

w 

L 


s 0. 390 


(47) 


which is sufficiently less than I to conclude that the approximations will 
probably be valid to within about 10% error. However, for V = 10°, and 
(2<rw/D * 1/10, 
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sin Y 


2a 
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L 


-1.15 


(48) 


which is NOT less than 1. Consequently the approximations lose their 
validity as the grazing angle decreases and the seas become choppier. Not 
only does the expansion of Eq. (43) to second order in Ax^ become question- 
able, but so also does the expansion of tan^fl (Eq.(37) ) to only second order 
in Ax and Ay, and the neglection of the dependence on Ax and Ay of the 
other factors which enter into the expression for the diffuse multipath 
power reflected per unit area. In addition, the validity of the conditions 
which must hold in order for the development of the diffuse multipath com- 
putations performed by Durrani and Staras^® must be called into question 
at the low grazing angles. Unfortunately these authors apparently failed to 
realize this and extended their approximate results down to very low grazing 
angles for which all the assumptions and conditions used to derive their 
results break down. This leads them to the totally absurd conclusion that 
the scattered multipath power received by the aircraft (or orbiting satellite 
in their example) becomes negligible at very low grazing angles. What 
happens, of course, is that the diffuse multipath power drops and the 
specular multipath begins to dominate at the low grazing angles. In addi- 
tion, shadowing and multiple scattering effects begin to play an important 
role at low grazing angles. These factors were apparently all ignored by 
Durrani and Staras. 


In any event, it is reasonable to conclude that, for grazing angles 
between 10° and 20° for which the Rayleigh criterion for diffuse scattering 
is satisfied, the relative multipath time delay spread, ( T m d ” T ms^> * s on 
the order of 10% to 50% of the relative specular multipath time delay, T ms 
This delay spread must be compared to the period of the highest modulat- 
ing frequency of the transmitted signal to determine whether a single mod- 
ulation time delay, T ms , may be used in describing the received multipath 
signal, or if the received multipath signal must be described as the sum 
of various multipath signals each having a different relative time delay. 
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For the NAV/SAT application the period of the highest modulating 
frequency (320 kHz) is 3. 125 Msec. The maximum relative specular 
multipath time delays, T ms , occur for the maximum aircraft altitude. For 
an aircraft altitude of 20 miles and a grazing angle of 20°, T ms is 73. 5 
M sec, whereas at an aircraft altitude of 5 miles and a grazing angle of 10°, 
T ms is approximately 9. 4 ^ sec Assuming minimum and maximum time 
delay spreads of 10% and 50% of T ms , the range of time delay spreads 
will be from approximately 1 to 35 Msec. Since this is not a negligibly 
small fraction of the period of the highest modulating frequency a sum of 
multipath components having different modulation delays must be used in 
describing the received diffuse multipath signal- This is described fur- 
ther in the next section 


2. 3 Characterization of the Received Multipath Signal 

The received multipath signal is described in terms of the received 
direct-path as follows: The received direct-path signal may be written as 

e (t) = \Iz C sin {to t + d ( t ) + 6 ( t ) } (49) 

d ' c 


where C is the received direct-path signal power, w c is the transmitted 
carrier frequency in radians per second, d(t) is the direct-path carrier 
doppler phase 


d(t) 



2-nf (t)dt + <j> 
d d 


(50) 


(where f^(t) is the instantaneous direct-path carrier doppler frequency and 
4^ is the initial carrier phase), and 9(t) represents the received phase 
modulation, the doppler on which may be totally neglected. For the case 
of BINOR ranging, 9(t) is 


0(t) - p b 13 (t) (51) 

where 0 = 1.15 radians and b^(t) is a BINOR code waveform generated 

by hard limiting the sum of 1 3 coherently-generated square waves, the 
highest frequency square wave being at 320 kHz and the other square waves 
being at a frequency 2*^ times 320 kHz k = 1, 2, 3, . . . . , 12). 
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The received multipath signal is then written as the sum of a pure 

"steady’’ specular component, e (t ), and a diffuse component, e (t). 

ms md 

e (t) = yJ2C a sin {u t + d(t) + d (t) + © (t - t )} 

m v ms c ms ' v T m« J 



E 


m 


x (t, t ) sin {u> t + d(t) + 
cm c 


d (t) + e ( t - T )} 

ms m 


+ 



m 


x (t, t ) cos {u< t + d(t) + d (t) + 6( 
q m c ms 


t-T )} 
m 


(52) 


where a is the mean relative specular multipath power; d (t) is the 
ms , , , , , . ms 

relative specular multipath carrier doppler phase 


d (t) 
ms 



2tt f (t) dt + 
ms 


ms 


(53) 


(where f (t) is the instantaneous relative specular multipath carrier dop- 
pler frequency, and 't’ms i s the initial phase); T ms is the relative specular 
multipath time delay; x c (t, r m ) andx^(t, x m ) are independent zero mean 
Gaussian random variables having identical variances, anc ^ T m is 

the relative time delay associated with the diffuse multipath components. 
The sum over the various relative multipath delays means that delay inter- 
vals &T m are being considered such that &T m is at least as large as the 
period of the carrier, but is much smaller than the period of the highest 
modulating frequency component in 0(t). The relative received diffuse 
multipath power is 


1 ) 

2 ) 

3) 

for all 


2cr =2 y a ( T ) 

m t—J m m 

T 


(54) 


m 


Additional properties of x (t, t ) and x (t, t ) are: 

cm q m 


( x (t, t ) x (t + t . t + At )) = 0 for all t and A r 
c m q mm' m 


(x (t, T ) X (t + T. T )) = (x (t, T )X (t+ T, T )) = 4 >(t, T ) 

' c m c m'^qmq m ' m 


(x (t, T ) x (t + TjT + At )) = (x (t, T )x (t+T, t +At )) = 0 

N c me m m /N qmq m m ' 


t and At > 6t . 

m m 
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The first property expresses the statistical independence of 
T m ) and x q(t + T « T m + A-r m ) for all t and A r m . The second property 
indicates that the autocorrelation functions of x c (t, T m ) and Xq{t, T m ) are 
equal, and finally, the third property states thatx c (t, T m ) and x c (t + t, 

T m + Ar m )> and Xq(t, T m ) and xq(t + t, T m + A Tm ) are also statistically 
independent for all t and for all A r m greater than & Tm . 

For purposes of analysis of the NAV/SAT receiver performance 
in the presence of multipath, the diffuse multipath signal may be written 
as 

e (t) ~ V2C x (t) sin { u t + d (t) + d (t) + 6 (f-r )t 

md c 1 c ms ms ’ 

+ C x (t) cos | oo t + d (t) + d (t) + 6 (t-T )l , rr , 
q 1 c ms ms > (55) 


where x (t) and xq(t) are statistically independent zero-mean Gaussian 
random variables having identical variances given by c r m ^. The power 

density spectra of x (t) and x (t)are identical and equal to 

c q 


$>. 


m 


(oj) 


2<j 


m ) 


(f/B )‘ 
m 


-s/JT b 


(56) 


L <f> 

where 2cr m and B are given by Eqs. (24) and (28). Note that m (co) is 
l/2 that of ^md(“) given in Eq. (28). This is because the total power in x (t) 
or x (t) iscr m 2 ; while the total relative diffuse multipath power is 2c . 


This representation of the diffuse multipath signal, e mc ^ (t), over 
estimates the undesired effects of the diffuse multipath component on the 
performance of the NAV/SAT receiver if worst-case values of T mg are used. 
Consequently, it will be used in the Task II analysis. This is so because 
the representation using a single relative multipath time delay is a more 
coherent interference than the representation which considers the entire 
spread of delay. 
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2. 4 Aircraft Antennas for Communication Links 

Employing Navigation Satellites 

2.4.1 Introduction 

The selection of suitable aircraft antenna types for use with naviga- 
tion satellites is discussed in this section. Both subsonic and supersonic 
aircraft applications are considered in choosing suitable types. The 
difference in speed between these two aircraft models, however, is not such 
as to allow the use of highly protruding antenna types at the lower speed. 
These types of antennas offer greater flexibility in pattern synthesis than 
do the flush mounted types most suitable for supersonic use. Thus the choice 
must ultimately be made from a group of antennas which are either flush 
mounted or else protrude from the aircraft surface without disturbing the 
aerodynamic integrity of the aircraft. This latter requirement implies 
that only low profile types can be used. The conical spiral antenna, for 
example, is ruled out by the low profile requirement even though it offers 
many advantages in flexibility of pattern synthesis. 


^ There has of course been substantial previous work in this area, 
and commercial types supposedly optimized for satellite use have 
already appeared.^ This previous work has attempted to optimize a 
design tradeoff between multipath rejection, achieved by suppressing radia- 
tion on or below the horizon, and substantial hemispherical coverage 
directed at good reception from several satellites which, it is felt, improves 
system accuracy in some instances. 


As far as can be determined at present, simultaneous reception of 
signals from multiple satellites is not contemplated as a standard opera- 
tional regime. According to the best information available concerning 
planned systems, satellite transmissions will be processed in sequence 
even when more than one satellite is visible at any given time. Thus the 
former rationale for the choice of a nearly hemispherical radiation pat- 
tern no longer applies, and other considerations are used to determine a 
desirable radiation pattern and an antenna design capable of producing the 
desired pattern. 


In view of the analytical work carried out under the present phase 
of the contract, and in view of supporting data which are now available, it 
appears that the most important determinant of the radiation pattern (and 
the antenna design suitable for the synthesis of this pattern) is the require- 
ment for multipath suppression from angles near to and below the horizon. 
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Table 2 


Aircraft Antenna Characteristics 


F requency 
Response 


1540-1660 MHz; Bandwidth exceeds BINOR require- 
ments but is chosen to be compatible with projected 
ATS communications experiments 


Polarization 


Circular 


VSWR 


1.50 maximum 


Axial 

Ratio 

Radiation 

Pattern 


Gain 

Beam 

Squint 

Sidelobe s 


Powe r 

Handling 

Capacity 


< 1 dB on axis 
5 2 dB at ±45° from axis 

Axial, circularly symmetric pattern. 
Beamwidth 

3 dB 80 - 90° 

10 dB 150 - 160° 

~ 5 dB on axis referred to a C. P. isotrope 


<5° 


<-15 dB over ground plane. Not seriously 
modified in operational environment 

Receive: no requirement 

Transmit: 200-500 watts average power 


2.4.2 Aircraft Antenna Requirements 

The requirements for the antenna element are listed in Table 2. 

The assumption is made here that the aircraft is in nominally level flight. 
The effects of aircraft motion will be considered later. The influence of 
the aircraft environment on the desired free space pattern will also be 
considered. Some of the tabulated antenna design parameters are 
discussed in more detail in this section. The radiation pattern is discussed 
in the following section. 
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The operating frequency has been chosen in accordance with 
present plans for L-band NAV/SAT operation. The bandwidth' exceeds that 
required for the BINOR waveform which occupies less than ±5 MHz band- 
width but has instead been chosen equal to the bandwidth of the projected 
L-band ATS-E and later ATS experiments. The required 100 MHz band- 
width is not a difficult design objective and is easily met by several types. 

Circular polarization has been chosen as usual for this application. 
Partial depolarization, or reversal of handedness of polarization, occurs 
after reflection from the ocean surface and the reverse polarized signal 
is not accepted by the circularly polarized aircraft antenna. This may 
provide additional multipath rejection capability in addition to that given 
by the antenna radiation pattern itself. 

Another element of the polarization requirement is the axial ratio, 
i. e. the ratio of major to minor axis of the (in general) elliptically polar- 
ized antenna. This is a measure of the deviation from exact circular 
polarization. For almost every antenna type which produces an axial 
pattern, circularity is best on the antenna axis and falls off monotonically 
with angle. The axial ratio is a measure of polarization loss and this loss 
will depend on the relative orientation of satellite and aircraft antenna axes. 
A careful power budget for such a link should take into account the polari- 
zation loss. A maximum axial ratio of 2 dB within a 90° cone has been 
specified. 

The voltage standing wave ratio (VSWR) is specified as l. 50 maxi- 
mum and this specification is easily realizable with several antenna types 
over the specified frequency band. 

Beam squint has been specified to be less than 5°. This specifi- 
cation is included here because of the importance of maintaining pattern 
symmetry with respect to the antenna axis so that deviations from the 
planned antenna angular coverage near the horizon will not occur. Beam 
squint arises when balanced antennas (e. g. spirals, dipoles) are fed from 
unbalanced lines. In making the transition from an unbalanced line (e. g. 
coax) to a balanced antenna structure, there exists on the balanced 
structure an in-phase component of current as well as the normal anti- 
phase or transmission line mode currents, unless^ an extremely careful 
transformation from feed to radiator is made. This in-phase drive 
excites higher order modes of radiation from the antenna which add 
vectorially with the normal radiation to produce a squinted main beam. 

The allowable beam tilt is specified to be 5° maximum. 


41 



2. 4. 3 Radiation Pattern Characteristics 


As has been mentioned the main emphasis in this work has been 
the control and suppression of multipath effects at the aircraft terminal 
due to reflections from ocean surfaces. The aircraft antenna radiation 
pattern is the main instrument used in effecting multipath suppression and 
will be specified here as carefully as possible. Again, the aircraft is 
considered to be in nominally level flight and the influence of the aircraft 
structure will not be considered in detail. It will be assumed that the 
radii of curvature of portions of the aircraft near the antenna are large 
compared with the carrier wavelength. This specification is not sufficient 
to eliminate any influence by the aircraft on the radiation patterns and, 
infact, in a detailed antenna design these effects cannot be disregarded. 

In carrying out a detailed design it is necessary to simulate the presence 
of the aircraft using either a mock-up of neighbouring metallic surfaces, 
or aircraft and antenna scale models in conjunction with frequency 
scaling . 


Table 2 gives the specification of the radiation pattern. The 
desired pattern is a beam pattern with circular symmetry about the beam 
axis. The 3 dB beamwidth is set at 80-90° and the 10 dB beamwidth at 
150-160°. The sidelobe level is specified as no greater than -15 dB when 
patterns are taken over a ground plane and this specification should not be 
seriously modified by the aircraft environment. The pattern is illustrated 
in Fig. 8 for an antenna mounted on the dorsal surface of a modern jet 
aircraft. 


The radiation pattern shown is not an ideal pattern for the satellite 
aircraft link since signals from 10-45° above the horizon are suppressed 
from 3 to 10 dB. It is however well-suited for multipath suppression since 
signals from +10° are suppressed by at least 10 dB, and for signals from 
satellites at higher elevation angles multipath rejection will be greater 
than 15 dB due to antenna directivity alone. Suppression due to polariza- 
tion reversal of the multipath signal will add to this 15 dB figure. In 
addition, this radiation pattern can be implemented in a flush mounted 
configuration. 

The gain has been specified as a nominal 5 dB referred to a 
circularly polarized isotropic antenna. This specification is not critical 
and must be consistent with the specified pattern shape. The precise 
antenna gain can be determined from the full radiation pattern taken either 
on a scale model or with a mock-up of nearby aircraft surfaces as 
described earlier. 
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2. 4. 4 Selection of Antenna Configuration 


The performance characteristics of Table 2 constitute a 
difficult specification for the user aircraft antenna. Several antenna types 
do have the potential for meeting this specification if no constraint is 
placed on protrusion of the structure above the aircraft surface. However 
supersonic aircraft users cannot tolerate appreciable protrusions above the 
surface because of frictional heating effects and aerodynamic disturbances. 
This factor eliminates most candidate antenna designs which are not flush 
mounted. In particular, spiral antennas mapped onto conical structures 
are eliminated unless they can form part of a vertical tail member and 
maintain an unobstructed field of view. The advantages of the conical 
spiral are the number of design parameters available, including cone angle 
and spiral pitch, and the freedom from the requirement for a ground plane 
which lessens the dependence of the pattern on the aircraft structure. 
However, if this and other medium and high profile types are eliminated, 
the choice of. a suitable element featuring good multipath rejection in 
accordance with the specification of Table 2 is made from a list of flush 
mounted or low profile types, similar to that considered in the Interim 
Report of the present contract. These include cavity backed planar spirals 
(Archimedean or logarithmic i. e. ), slotted flush mounted dipoles, curved 
dipole turnstiles and curved dipole crossed-slots. None of these types can 
be ruled out categorically when only L-band operation is considered since 
maximum protrusion above the aircraft surface is 2-2 1/2 inches for 
operation at 1600 MHz. The choice is dictated instead by considerations 
of radiation patterns. With the exception of the planar spirals, all of the 
types mentioned above produce hemispherical or near-hemispherical 
patterns (at least 0 dB over a 160° cone, referred to a C. P. isotrope). 

Such a pattern puts too much power on or near the horizon and all types 
producing such near -hemispherical patterns must be discarded. 

Thus the cavity-backed planar spiral antenna emerges as the 
recommended choice for the SST - satellite link with good inherent multi- 
path rejection capability. The choice between logarithmic and Archimedean 
is made on the basis of design and production costs. The Archimedean 
spiral is marginally simpler to design and fabricate and is recommended. 

A cost advantage might be obtained, however, from manufacturers with 
more experience in the design of logarithmic spirals and this latter ap- 
proach is clearly acceptable. 

Performance capabilities of the Archimedean two wire spiral 
antenna are compatible with the specification of Table 2. Most applica- 
tions of the Archimedean spiral, in particular ECM use, have capitalized 


44 



on the wide bandwidths inherently available with this type. Manufacturers' 
catalogues show spirals with specifications maintained from 2-11 GHz. No 
commercial type appears to meet the present requirements completely but 
manufacturers see no deterrent to meeting the present specification com- 
pletely. The theory of operation of this antenna is given in Section 3 . 4 of 

the Interim Report. The theory and applications are discussed in detail 
19 

by Kaiser- 7 

2. 4. 5 Use of Switched Antennas 


The formulation of the desired antenna pattern assumed that the 
aircraft was in level flight and that the satellite being utilized was essen- 
tially overhead. The effects on the multipath rejection capability of the 
aircraft antenna system of satellites at elevation angles well below the local 
zenith were not considered. Similarly when aircraft banking takes place 
during turns, the antenna being used in level flight operates in a less than 
optimum way. Thus the appearance of the satellite at low elevation angles 
increases the power arriving from a sector close to the horizon and thus 
potentially increases the multipath power received by the antenna; depending 
on the antenna pattern and aircraft orientation. 

The addition of two antennas identical to the one specified in 
Table 2 can overcome most of the deleterious effects caused by the satellite 
appearing at low elevation angles relative to the aircraft. These additional 
antennas would be arranged circumferentially around the aircraft at the 
same longitudinal station as the primary antenna. This location would be 
chosen on the basis of adequacy of ground plane, absence of shadowing ob- 
structions, etc. The situation is illustrated in Fig. 9. This system would 
be used as follows. When the satellite is at high and intermediate elevation 
angles antenna 2 is used in order to suppress multipath pick-up as much as 
possible. It should be noted that the angles of arrival at intermediate eleva- 
tion angles are such that the direct signal will also be attenuated somewhat 
due to the pattern shape, and this should be taken into account in a detailed 
power budget. At satellite elevation angles below about 43°, antenna 1 or 3 
(depending on aircraft heading) is used in order to maximize gain in the 
satellite direction. The appropriate antenna is selected by a simple switch- 
ing arrangement. 

In view of the low cost of the recommended basic antenna this 
approach appears to be an inexpensive and practical solution to the problems 
presented by the appearance of the synchronous satellites at various eleva- 
tions with respect to the aircraft. And, in fact, even more elaborate ar- 
rangements of antennas and switching equipment can be envisaged in order 
to optimize multipath rejection by a simple, inexpensive antenna element. 
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Another possible configuration would use non -identical antennas. 
Primary coverage of a large overhead sector would be obtained using a 
spiral of the type discussed previously. Coverage of satellites at lower 
elevation angles would be obtained by the use of simple horn radiators 
(incorporating suitable arrangements for circular polarization) a few wave- 
lengths in aperture width, producing a pattern with a 20-30 degree half- 
power beamwidth. This approach is practical at 1. 6 GHz because of the 
aperture width. Switching arrangements would be similar to those men- 
tioned above. 



P-S4J1 


Fig. 9 Switchable Antenna Patterns to Compensate for In-Flight 
Maneuvers 
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3. BINOR ACQUISITION ANALYSIS 


3. 1 Summary 

Calculations for the probability of correct acquisition of 
the BINOR code are made for three cases. In each subsequent case 
additional factors are included in the calculation. The three cases 
are: 


a) Case where effects of multipath, carrier and 
clock loop phase jitter, and video filtering 

of the BINOR code in the receiver are neglected, 

b) Case where effects of multipath are neglected 
but effects of carrier and clock loop phase jitter, 
and video filtering of the BINOR code in the 
receiver are included, 

c) Case where multipath effects as well as the other 
effects are included. 

The third case is, of course, the most complex, and is 
best understood in terms of the modifications required in the 
results of the first two cases. 

Ranging accuracy is shown to be a function of the rms 
phase jitter in the clock loop. This phase jitter is shown to result 
from the presence of Gaussian noise and multipath components in 
the loop noise bandwidth. 

3. 2 System Model 


A model of the NAV/SAT receiver is shown in functional block 
diagram form in Fig. 10. The functions of the various blocks are 
described in terms of their signal processing functions. These are most 
easily discussed in terms of case (a), in which multipath, and effects 
on acquisition of the BINOR code of phase jitter on the extracted carrier 
and clock loop reference signals are ignored. In this case the input RF 
signal may be written as 


'RF 


1 ! 2 

(t) = { 2C ) sin Tco t + d(t) + /3b. 

L c 13 


(t)]+n RF (t) (57) 
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where 


C is the total received signal power; 

9 

is the carrier frequency and is 2n(l. 6X10) rad/ sec. ; 
d(t) is a slowly varying Doppler phase; 

b i ^ ( t ) is the BINOR code. It is generated by hard limiting 
13 coherent equal amplitude square waves, each of 
which has a frequency 1/2 that of the next higher 
frequency square wave. The highest frequency square 
wave 'serves as the /'clock" and has a frequency, o> 
equal to 27 t(320 x 10°) rad/ sec. The frequency of the 
lowest frequency square wave is 2^(78.125) rad/ sec.; 

jS is the phase deviation of the BINOR code on the RF 
carrier and is given as 1.15 radians; and, 

n^ F (t) is "white" Gaussian noise having double- sided noise 
power density denoted by <&/2. 

The worst- case ratio C /$, including a 6 dB margin, is given 
as 35. 7 dB-Hz. 

The function of the carrier tracking loop is to acquire and track 
the carrier component of ej^pft), and use the extracted carrier as a 
reference signal to phase demodulate e^jr{t), thereby producing the 
video signal e v (t) given by 

e v (t) = (C) 1/f “(sin/S) b 13 (t) + n v (t) (58) 

where n v (t) is "white" Gaussian noise of double- sided noise power density 
equal to $/ 2. 
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This signal is bandpass filtered to extract the clock fundamental 
frequency, u This differs from an alternate configuration in which 
the clock loop is driven directly from the output of the lowpass filter. 

A careful analysis has shown that the phase variance for the alternate 
configuration is double that given by Eq.( 61 ) below. Thus, by including 
the bandpass filter, performance is improved by 3 dB and approaches 
within 1 dB of the performance previously, but incorrectly assumed for 
the alternate configuration. 

The output of the bandpass filter may be written as 

e BP (t) = (C)^ 2 sin^p 13 (^) sin w^t + n Bp (t) (59) 

where p^is the correlation coefficient between bj 3 (t) and ci Q (t), the 
highest frequency square wave. (The correlation coefficients between 
b 1 3 (t) andeach of its component square waves cj£ k (t) (k = 0, 1, •••, 12) are 
equal to p 13 . ) The factor 4 / n is the relative amplitude of the first har- 
monic of cJ 2 0 (t):‘ ngp(t) is Gaussian noise which has a flat noise power 
density (double -sided) of $/2 over a bandwidth at least as wide as the clock 
loop noise bandwidth B c< , centered around u c ^ (and -w ^). 

The function of the clock tracking loop is to acquire and track 
the clock frequency u> c j?. The presence of ngp(t) results in a phase jitter 
on the extracted clock reference signal, <|) c j? n j which is approximately 
Gaussian distributed and has a variance given by 



( 60 ) 


where 2 B ^ is the double- sided noise bandwidth of the clock loop and is 
determined by the requirement that the loop SNR is 10 dB at the worst- 
case value of (C/$): 


(snr) l 



2 . 2 
p 13 Sm 3 



10 


(61) 


That p 13 (4/tt) is the correlation coefficient between b 13 (t) and sin w c ^t 
may be seen by computing (b 13 (t)ci (t + r)) which is found to equal 
p 13 (c^ 0 (t)ci 0 (t +r)) for all r. <b 13 (t)sin w c £ k t> = 4 /tt p 13 is also true 
for c-fi j(t), but not for the other ci k (t) (k 2 2). 


50 



With (C/4>) 0 = 35. 7 dB-Hz, and p-^ 3 = 0. 225 586 . . . , the required value of 
(2B c _g) is approximately 12. 8 Hz. With this noise bandwidth the mean 
time to acquisition of the clock loop is approximately 0. 4 sec. 

Combining Eqs.(60V and(61)the phase variance may be written as 


/, 2 \ J_ Wp 

\ c 2 n/ 20 (C / 4) 



(62) 


during acquisition of the BINOR. During the range measurement the 
BINOR is turned off and only the clock square wave is transmitted (with 
phase deviation 8). If an AGC system is used to decrease the signal 
gain in front of the clock loop in order to maintain the same loop noise 
bandwidth and damping coefficient, the loop SNR is modified by replacing 
Pl3by 1 in Eq.(61) above. The ^result is that (SNR) L becomes equal to 
approximately 23 dB, and( 


,± I JC vl 

<‘h> 



(63) 


The variance of the clock loop phase jitter will be related to the rms 
range error in a later section. For purposes of acquisition analysis 
in case (a), the phase jitter on the extracted clock reference signal 
is ignored. 

The phase shifter serves to compensate for phase shifts introduced 
in the video filter in a manner which will be explained in the discussion of 
case (b). Since the video filter is ignored in case (a), the effect of the 
phase shifter is also ignored, and the signal denoted by eppW in Fig.10 
is simply e v (t). 

The BINOR code acquisition system serves to generate a square 
wave whose period is the same as that of the BINOR code and which is 
in phase with the lowest frequency square wave of the BINOR code. This 
is accomplished by successively frequency dividing down the extracted 
clock reference and comparing the phase of the resulting square waves with the 
phase of the corresponding square wave component of the received BINOR. 

The first divided down square wave is either in phase or 180° out of phase 
with the corresponding square wave component in the received BINOR. 
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The phase is determined by multiplying the received BINOR by this 
divided down square wave, integrating the resulting product for a 
time AT (where T is the period of the BINOR code and A is an integer), 
and sampling the output of the integrator. If the sampled output is 
positive, the square wave is assumed to have been in phase with the 
corresponding square wave of the BINOR code, whereas a 180° phase 
relationship is assumed if the sampled output is negative. If the square 
wave was found to be in phase, it is frequency divided by two and the 
above process continued. If the square wave was found to be out of phase, 
it is first complemented and then frequency divided by two and the above 
process continued. If all decisions are made correctly the final frequency 
divided down square wave will be in phase with the BINOR code. This 
square wave is used in the range measurement system. 

3. 3 BINOR Code Acquisition Analysis - Case (a) 

The purpose of an acquisition analysis is to determine the 
probability of correct acquisition in terms of the various parameters 
involved and to find out which value of A must be used in the integrate 
and dump operation to insure at least a 99. 9% probability of correct 
decision for the worst-case C /$. 

On the assumption that the first (k - 1) decisions were made 
correctly, the output of the integrate and dump circuit at the end of the 
kth square wave acquisition is for case (a): 


= ±(s,. + 


, AT 

V = (/ 


(C) 1 / 2 sin/3 


b 1 3 (t )c^^(t)dt 


+ J n y (t) ci k (t) dtj (64) 

o 

The "+" indicates the in phase case, and the'"-" the out of phase case. 
c^ k (t) is the kth frequency divided down square wave and is assumed in phase 
with the corresponding BINOR square wave component . The noise component 
of e^ is denoted by n^. and is Gaussian distributed with variance cr^ . 

The signal component of e^ is denoted by s^. and is 


1 1 2 

AT(C) sin/3 



AT 

I b 13 (t)ci k (t)dt 

o 


(65) 
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The integral in brackets is the correlation coefficient between b^ (t) and 
ci, (t) and is equal to for all k, so that 


s^ = iT(C ) ly/2 sin£ p^ = p 


( 66 ) 


2 

for all k. The variance a, is 

k 

IT IT 


°l~- / dt l / d k C VV C VV ("vV'VV) (67) 

0,0 


Since n^(t) is assumed "white", 


(w W) = I 6 <1 1 ■ V 


( 68 ) 


where 6(t) is the unit impulse function. In this case 0^ easily seen to 
be given by 



(69) 


for all k. 


The probability of correct acquisition of the kth divided down 
square wave, given that all preceeding (k - 1) divided down square waves 
were correctly acquired is 


( p > = /” p<v 

k - s k 


(70) 


where P(^) probability that the noise n^ lie between rfy. and 

77 ^ + d77j_. The above is simply the probability that the sum 
be greater than zero, since only then will the sign of e^ be the same 
as the sign in front of the brackets in Eq. (64). For a Gaussian distributed 
n k» P(i7k> is 


p<y, =, 


1 

V2tt tr. 
k 




(71) 
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so that 


» k ok 

(P rK = / p( V d \ = / p( V d \ = J P( 7 } k )d \ + / P( V d \ 


-s, 


- 00' 


V ~ 'k j ' V 'k 

- CO O 


'\+\iT J [ 1 + erf M <72 > 

0 

where C = \l and = s yj J2 a k = /i/,/2 0. The total probability 

of correct acquisition, (P c )rp, is 


12 

(P ) = TT (P ) 

c c 

T k=l k 


\{ 1 + erf(M/V2ff)| 


(73) 


for case (a). 


The ratio J2 a is 



M/v^ CX = p 13 sin/5 y^T 'I~ 

(74) 

At the worst 

case C/$, this is 



(/j/a/2 a) = p sin p ./IT -/(P") 

O io St' o 

(75) 

■j^ 

with the correlation coefficient given by' 1 ' 


P 13 

( 2 ) ( 4 ) (g) (g) ( 1Q ) ( 12 ) - 0.225586... 


sin/J 

= sin(l. 15) = 0.912764 


Vt 

- 1 / 2 

= (78.125) - 0.113137 


(C/*) o 

- 35. 7 dB-Hz = 3715. 35 



- 60.9537 



(HJJ2 0) = 1. 41996 JT ~ 1.420 VT 

o 

(76) 

Stiffler. 
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0. 001 0. This is 


For (P ) - 99. 9% - 0. 9990, (P) - 1 - (P ) = 

approximately equal to 

(P ) « 12 -kl - erf(u IJ2 a) 1 - 0.0010 

e T 2 L ^ o' 

or 

[ 1 - erf (ji/V2 a) ] » 0.00016667 = eric{y.l J2, a) Q (77) 

By using tables for erfcx (c. f. Abramowitz and Stegun, Handbook of 
Mathematical Functions, Dover Publications, 1965, p . 316), and using 
linear extrapolation, the required value of {\ii <J2 cr) Q is found to be 

^ 2. 663 (78) 

o 

min 



which agrees with results obtained for this case by others. 
For f = 3 and 4=4, (pi J2 cr) Q equals 


(tStf), 


2. 459 4 = 3 


2. 840 


(79) 


Consequently the required value of 4 in this case is 4 = 4. 


3 . 4 BINOR Code Acquisition Analysis - Case (b) : 


In this case the additional effects of carrier and clock loop phase 
jitter and video filtering are included. Under the assumption that both 
carrier and clock loops may still be treated as linear systems in the 
tracking mode (corresponding to the assumption of small phase jitter 
and/ or tracking error), the carrier loop phase error consists of the sum 
of a tracking error, denoted by and a Gaussian distributed phase 
jitter, denoted by <^ , which results from the presence of Gaussian noise 

at the input to the carrier -tracking loop. In a second-order loop, the 
tracking error is 




ct 


_ !<t> . “c rtt) 


2 c 2 

00 00 
n n 



( 80 ) 
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where 


CO - 

n 



( 81 ) 


and 2B C is the double -sided loop noise bandwidth (2B C = 50 Hz) and C c is 
the damping factor (C c ~ l/*/ 2). The factor r(t) is the acceleration along 
the aircraft-satellite line -of-sight vector in ft/sec^, and "g" is the 
acceleration due to gravity (g = 32. 2 ft/ sec^). With f =1.6X10® Hz, 

<t>ct is 


♦=t • °- 14 %) 

The variance of the phase jitter 4> is 


(82) 


<t> ■ (D 


2 a 
cos p 


B 


(C/4) 


cos 


2 0(§) o (c/ *> 


(83) 


(C/4) 

= t°- 04033 1 

The sum of <f> c ^ + <f> cn is denoted by <J) C . The video output, e y (t). 


is now 


e (t) = (C) 1 ^ 2 sin 0 cos (j) (t)b (t) + n (t) 
v T c 13 v 


(84) 


and the output of the bandpass filter at the input to the clock loop is 

e „ D (t) = (C) iy 2 sin£ cos (j> (t) - p sin co t + n (t) (85) 

BP c 7r 13 cJl BP 


If 4> c (t ) is small, its effect on the clock phase jitter through the 
term cos<j> c (t) in Eq. (85) may be neglected, so that \is to a first 

approximation unchanged from the results given in Eqs. (62) and (63). 
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The output of the lowpass filter, whose impulse response is 

h LP U) iS 

p (t) = < C )^ 2 sin £ cos + c (t)|b 13 © h Lp j + (n^ 0 h Lp ^ (861 

where 0 denotes the operation of convolution. 

The output of the k^* 1 integrate and dump circuit (Eq. (8)) is 
modified to 



where the effects of carrier loop phase jitter, video filtering, clock loop 
phase jitter, and the phase shifter delay 


t 

o 


<t> 


ci 


<JL> 


ci 


have been introduced. Since these effects are small, they may be treated 
as perturbing influences in the calculation of (P c )^ and < P c>T- in that only 
their dominant or lowest order effects need be considered. 

2 

The noise is still Gaussian, although its variance cr^ is slightly 
modified as a result of the combined effects of video filtering and multi- 
plication by the kth frequency-divided-down clock. The effect of clock 
jitter on the computation of is negligible and may therefore be neglected. 
This is shown in Appendix A where a 2 is written as 

222 

a k ■ a r k (ss) 

2 

and where 7 k is seen to be slightly less than 1 for the case where h ^(t) 
is a single -pole lowpass filter whose 3 dB bandwidth is made equal 
to the fifth harmonic of the clock. That is 

h Lp (t) = a e u^t) (89) 
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where a - 5w c ^,. A reasonably good approximation to (P c )j t may be 
obtained by using the expected value of s^ in Eq.(72). This over 
estimates (P )^ somewhat, but only slightly. A more exact treatment 
would involve computing the expectation value of (P c )v. One way of 
doing this is to expand (P c )^ in a Taylor series about (P )^ evaluated 
at the expectation value of s^. This would require knowledge of the 
higher order statistics of s^, and begins to get quite involved, even 
for second order statistics. Consequently it was felt that since <f> c (t) 
and <J> Cjg (t) are small, a more exact treatment would not contribute 
anything of significance in relation to the effort which would be required. 


The expectation value of s 


k 


is shown in Appendix B to be given by 



lie 


1 2 
2 a cn 


cos 


*ct E 






(90) 


where 


and 


V T) 


V 13 w 


- iT(C) 1 / 2 p 13 sin)3 

(91) 

-a 2 = E, fcos <b \ 

1 cn cb \ cnl 

cn ' ’ 

(92) 

o 2 = 4 2 \ 

cn y cn/ 

(93) 

/" d ch «)* (T-n 

-oo k 1 3 

(94) 

i7 / fT b 13 ,t)c V t + T)dt 

o 

(95) 


The phase shifter time delay, t Q , is chosen so as to maximize 
Xp. when = 0. This choice of t Q yields the maximum /s^N possible 

and also minimizes the effects of <t 5 c ^ n on acquisition. ' 


Since t Q is small, and since | <f> c ^ | is less than 90° with very 
high probability for values of (C/9) greater than (C/4>) 0 , it is only neces- 
sary to compute Xp.( r ) for | T | < T ^,/4 ” 2 tt/4w c ^. For h-^p(t) as given in 
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Eq. (89), transient behavior may be considered to have almost completely 
died out after a time t = T Cj g/4. Consequently, it is only necessary to 
compute </j for I 7 < T c 8.1 2 - In this ran g e * 11 ma y be shown" 

that 0 cjPkbl 3 (T) is given by 


T 


^ ci k b 13 (T> ' Pl3 \ " (1 " Sk> (T 


cl 


/ 2 ) 


for T < 


cl 

2 


(96) 


where a is the ratio of 
Table 3 nelow. 




/ 2 ) toil) 



( 0 ), and is given in 
13 


Table 3 


VALUES OF a AND (1 - a ) AS A FUNCTION OF k 
k k 


k ' 


(1 - a k > 

1 

0 

1.0 

2 

1008/ 1848 - 0. 545454. . . 

0. 454545. . . 

3 

1512/ 1848 - 0. 818181. . . 

0. 181818. . . 

4 

1736/ 1848 = 0. 93939. . . 

0. 060606. . . 

5 

1820/ 1848 = 0. 984848. . . 

0. 015151. . . 

6 

1844/1848 = 0. 997835. . . 

0. 0021645. . . 

7-12 

1. 0 

0 


For the single -pole lowpass filter whose impulse response is 
given in Eq. (89), Xk( T ) * s approximately 

T 

X k (T>3 1 + a T ~ l 2 (1 ~ a k )(QT ' *> --f i £T<0 (97) 

cl 


This was done by computing the BINOR code b^ bit by bit and evaluating 
i/Jcikb]^ 1 ") a i 7 ~ 0 and 7 = ±T c ^/ 2. This calculation of the BINOR code 
b ^3 and ^cj^bi 3 ^) was facilitated by making extensive use of readily 
derived properties of BINOR codes and the functions ^c.gvbi oCO. A dis- 
cussion of these calculations is somewhat lengthy and is therefore not 
presented. 59 




1 


ci 


1 -CtT 

(1 - a )[(oT - 1) + e ] 


° T J 4 


T 


0 < T < 


ci 


( 98 ) 


where 


and 


a - 5w 


1 Or 57 r 


T ci> V 2 


r = t - <b „ , 

o *ci ' ci 
n 


( + <*o ■ 


(99) 


(100) 


Introducing Eqs. (99) and(100)into the above expressions for X^( T )> 
Eqs. (97) becomes 


( T ci v 

“oT" ’ ‘o) ’ 1 + 5T ( 1 " ' l ) ' * U ’ 


( 101 ) 


J) „ > <t „ 
m JL 

n o 


and 


X k( — ^ 


CO 




0 = 1 + 5^ (1 ' a k ) ( 5 *ci 0 " X ) '^ a ■ 


-^ (1 - a k ) ( 5 ^ - I+e % n ) + !»-^.)t. 


K ' k r cj£ 


*ct < Kl 

n o 


n 

( 102 ) 


With <|> CJ g = 0, <(> c i 0 is chosen so that X k ( -t Q ) = (X^niax- 
Differentiating Eq.(102) with respect to <j> c ^ Q and setting the derivative 
equal to zero, and putting <f> c jg n = 0, <J> C ^ 0 is found to be 

+ci = ? fn2 <103) 

o 
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and (v ) is 
k max 


(X ) = 1 - r- (1 - a, )in2 

K max D77 k 


With £ k = it (1 - V- and0 L = (+CJ )• l \ 

n \ n/ r ci 
n 


(104) 

{Kl ) 


/W 


of Eq.(90) is shown (in Appendix C) to be given by 


n' 


n o 


1 

\ 

) /“> , 

1= 1 - 1 

/ ci 

/ k 


/ 


- (1 - in2)erf/ 


in2 


2 5a 


c i 


n' 


n o 


ci 


in 2 


/oa 


c i 


V 2 0a 


n I 2 


ci 


n' 


+ e 


V 2 . / k/ff " 

r A , 

XT 2 / 2 


/5 a 


ci 


in2 


1 “ Gl'fl ‘ 

\V 2 V2 5a ci 


Alternatively, ^ X ^ 


ci 

n 


{KrKt ) 


n o 


/“ | J may be written as 

>-« Xj 


n> 

( 105 ) 
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<p 


ci 
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( \ 


\ k 

\ 

n o J/ 


/5 a 


ci 


\f 2 5a 


' (X k max ~ f k| 


C^n) 

+ e^| 


n I 2 


ci 


n/ 


a/ 2" / tTa 


5 a \ l 


1 + erf 


in2 


ci 


t/ 2 " 5 a ~ \/2 
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ci 
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- erf i 
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- in2 


1 - erf/ 
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The probability of correct acquisition of the k^ 1 frequency- 
divided-down square wave component, given that all previous (k - 1) 
components were correctly acquired, is now given by 


( p >, 

c k 


1 

2 ' 


1 + erf 



r<x, ) 

k max 

y,. 



(107) 


and, again, the total probability of correct decision, (P ) , is, 
i c 1 

(P c»T = ft< P c ' k < 108 > 

k-1 


It is interesting to compute the ratio of (X^ 5 ) rnax to since the 
deviation of this ratio from 1 is a result of the lowpass filter h Lp (t). 

The effect of the filter was seen to (1) round the peak of the function XjJ 1 ') 
and reduce its maximum from 1 to (X^) max > and (2) reduce the variance 
of the noise component at the output of the integrate and dump circuit by 
a factor The rounding of the peak of Xv( T ) reduces the effects of clock 
jitter on acquisition, since, if the phase shifter phase <f) c ^ o is chosen 
properly, the first order effect of <j> CJ g n on acquisition probability is 
zero (corresponding to r(CT c ^ n ) having zero derivative at = 0). 

Consequently, if the ratio of (Xk) max to y k is not too much less than 1, 
the effect of the lowpass video filter hpp(t) is to improve the proba- 
bility of correct acquisition by reducing its sensitivity to small clock 
jitter. The values of this ratio are given in Table 4, where it is seen that 
this ratio is indeed very close to 1, and even exceeds 1 slightly for k > 4. 

Curves of (P c )p computed from Eq. (73) and Eqs.(108) and(l07) 
with <j> c £ = 0 are presented in Fig. IE as a function of (C/$)/(C/$) 0 . In 
these curves it is seen that a value of jC - 4 is required if a 99. 9% proba- 
bility of correct acquisition is to be achieved with a 6 dB margin allowance 
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(^ n -^ 0 )) 




<i. e., for (C/$) Q = 35. 7 dB -Hz). If the margin allowance is reduced to 
5 dB, the required acquisition probability may be achieved with £ = 3. A 
3 dB margin allowance will permit the required acquisition probability to 
be achieved with £ = 2. If no circuit margin is required, £ = 1 is sufficient 
to achieve 99.9% acquisition probability. 


Table 4 

VALUES OF (x, ) /y, FOR k = 1, 2, ••• , 12 

k max' k 


k 

(X ) 
k max 

y k 

(x kWA 

1 

0. 955873 

0. 967646 

0. 987834 

2 

0. 979942 

0. 983956 

0. 995921 

3 

0. 991977 

0. 992011 

0. 999966 

4 

0. 997326 

0. 996013 

1. 001318 

5 

0. 999331 

0. 998009 

1. 001327 

6 

0. 999904 

0. 999005 

1. 000900 

7 

1. 00 

0. 999502 

1. 000498 

8 

1. 00 

0. 999751 

1. 000249 

9 

o 

o 

0. 999875 

1. 000125 

10 

1. 00 

0. 999938 

1. 000062 

11 

1. 00 

0. 999969 

1. 000031 

12 

h-» 

O 

O 

0. 999984 

1. 000016 
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(Pc ) T (P« ) r 



Fig. 12 Correction to Acquisition Probability 

Resulting from Clock and Carrier Jitter 
and Video Filtering (0^ = 0) 
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3. 5 BINOR Code Acquisition Analysis - Case (c) 

3. 5. 1 Characterization of the RF Multipath Signal 

In this section a specular and a diffusely scattered multipath 
return is added to e^p(t) and their effects traced through the NAV/ SAT 
receiver. 


The specular multipath component of e^p{t) arises from coherent 
specular scattering of the transmitted signal from a very small region of 
the sea surface consisting essentially of the first Fresnel zone. This 
component of ep^p(t) may be written as 


e 



ms 


V~2C 


a sin[oj t + d(t) 
ms L c 



(109) 


where a ms is the relative amplitude of the specular multipath component, 
is the relative delay of the specular multipath component, and ^ mS (t) 
is the relative specular multipath Doppler phase given by 


d (t) = S'* co (t)dt + ij) 
ms J e 


-00 


ms 


ms 


( 110 ) 


Here u> ms is the relative specular multipath Doppler frequency, and ^ 

is the initial phase of d (t). 

ms 


The diffuse multipath component of e^jr(t) arises from incoherent 
scattering of the transmitted signal from a much larger region of the sea 
surface centered on the specular region. If the modulation delays are 
taken to be essentially the same as for the specular scattering, this 
component of epjp(t) may be written as 


e RF (t) md 



+ V 2C x^(t)cos[u> c t + d(t) 



( 111 ) 


where x Q (t) and Xq(t) are independent Gaussian random variables having 
identical variances denoted by cr^, and r ms is the delay associated with 
the specular multipath component. 


Additional properties of the random functions x (t) and x (t) are 

c q 

1) The cross correlation function /x (t)x (t + T)\ is 

zero for all T, and \ c 3 7 

2) The autocorrelation functions ^x c (t)x c (t + T)^ and 
(x (t)x (t + t)) are equal and may be denotea by 

<l> V). q 
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The autocorrelation function cj> m (T) and its Fourier transforms, 
$ m (cL>) depend on such factors as the aircraft velocity, the sea state, the 
grazing angle, the wavelength of the incident radiation, and the aircraft 
altitude. 


The presence of the multipath causes increased carrier and 
clock loop phase jitters and introduces extraneous components at the 
outputs of the integrate and dump filters. These can all affect the proba- 
bility of correct acquisition of the BINOR. The multipath interference 
contains both deterministic and random elements: This is so because 
it consists of the reflection of a deterministic signal off a random rough 
surface. It should be expected that both of these features of the multipath 
interference will play a role in determining its effects on system per- 
formance. This is indeed the case. 

3.5.2 Treatment of the Probability of Correct Acquisition 

In estimating the effects of multipath on system performance it 
is necessary to compromise between the desire for a complete statistical 
description of these effects and the effort required to attain such a 
description. In discussing the probability of correct acquisition of the 
BINOR code, consideration is limited to a computation of the expected or 
average probability of correct acquisition obtained by averaging- over 
appropriate random variables of the multipath interference. While this 
is not a complete statistical description of the probability of correct 
acquisition, it is still the most useful single statistical parameter in the 
case where the effects of multipath are small. It is also the simplest 
parameter to obtain, involving as it does no higher than second order 
statistics of the random multipath signal parameters. Confining the 
analysis to a computation of the average probability of correct acquisition 
is useful if the computation predicts that the effects of the multipath are 
small even when worst-case bounds are placed on the statistical parame- 
ters of the multipath interference, and if the objective is to determine 
whether or not the presence of a small multipath interference presents a 
serious problem. 

In addition to the average probability of correct acquisition, 
additional worst-case possibilities are considered and estimates made 
of the frequency of occurrence of these worst-case conditions. The 
probability of correct acquisition under these additional worst -case con- 
ditions are computed. 
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3. 5. 3 Effects of Multipath on Carrier Loop Phase Jitter 

The contribution of the multipath to the phase jitter on the 
carrier -tracking loop is now computed. The carrier component of 

e RF (t> iS 

e (t) = V 2C cos/3 sin[oj t + d(t)l + n (t) 
c L c J c 

+ V 2C cosj3 a sinlw t + d(t) + d (t)l 
ms L c ms 3 

+ V2C cos£ x (t)sinlu t + d(t) + d (t)l 

r L r mS J 


+ V 2C cosiS x (tjcoslw t + d(t) + d (t)l 
q L c ms 3 


(112) 


and where n (t) is Gaussian noise of power density 2 (double -sided) 


over a bandwidth much larger than the carrier loop noise bandwidth. 

The phase jitter of the carrier reference signal extracted by the 
carrier loop is now 


<t> (t) = <|> + <J> (t) + <j) (t) + <j> (t) 

T c ct cn T cs cd 


(113) 


where it is assumed that the loop is a linear system so that <j> c ^ and <|> cn (t) 
are essentially the same as for case (b). The phase jitter resulting from 
the specular multipath component is 4 ) cs<t)-> while that resulting from the 
diffuse multipath component is <J) C( }(t). Using the phase model of the loop, 
the input phase, from e c (t), neglecting n c (t), is 


4>^(t ) = tan 


a sinfd (t)l + x (t)sin[d (t)l +x (t)cos[d (t)l 
ms 1 ms 3 c L ms 3 q 1 ms 3 

1 + a cosfd (t)l + x (t)cos[d (t)l -x (t)sin[d (t)l 
mS Lrns J r. L m s J n L ms J 


ms 


(: 


14) 


which to first approximation is 

(t.(t) » a sinfd (t)l + x (t) sinfd (t)l + x (t)cosfd (t)l (115.) 
T i ms a ms 3 c a ms 3 q L ms J 

Higher order terms are assumed negligible as are the effects of the 

amplitude fluctuations of e c (t) . Assuming that the relative specular 

multipath Doppler frequency, changes only slowly with time, <{> cS (t) 

may be written as 

i (t) = | H (to (t)) [ a sin 

cs c\ ms / 1 ms 


d (t) - 
ms 


<t> /w (t) 1 ) 
c( ms J 


(116) 
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where | H c (u>) | is the amplitude of the loop transfer function and <j> (w) is 
the phase of the loop transfer function H c (w). That is, H c (u) is given by 


H (u) = |H (w)|e 
c 1 c 1 




(117) 


The phase jitter <f> cc [(t) is a random variable which is approximately 
Gaussian distributed (if the loop is treated as a linear system) and whose 
variance is denoted by cr^. If the spectral broadening of x c (t) and x q (t) 
is much greater than the carrier loop noise bandwidth B c , and if the 
power spectral density of x c (t)sirTjd ms (t)] + Xq(t)cos[d rns (t)] is approximately 
constant over B c , the variance cr^ may be written as 


cj 2 = $ <u )2B (118) 

cd m ms c 

where ^ > rn ( w 1 - ns ) is the power spectral density of x c (t)sin[d mS (t)] + 
x q (t)cos[d ms (t)] at ui = 0. 

Although only limited experimental data is available, there is 
some theoretical support to indicate that for aircraft velocities on the 
order of 600 mph the spectral broadening of the diffuse multipath com- 
ponent can be significant at L band, and that the power density spectrum 
of x c (t) and Xq(t) is Gaussian in shape. That is 

d 2 -<f/B ) 2 

V“ )= ?rf- e m <U9) 

m 


where the bandwidth B m depends on aircraft velocity; the sea state; the 
grazing angle; the wavelength of the incident radiation; and to a lesser 
extent, the aircraft altitude. The total power in x (t) and x (t) is a? 


Introducing Eq.<119) into Eq.(118), the variance is written as 


'nr 


CT 


2 

cd 



-<f / B )' 

ms' m 


( 120 ) 


To determine the relative importance of diffuse multipath in the carrier 
loop, is compared to cr^. 


69 



2 

a 

Ctl 


( 121 ) 



2 2 

The ratio of a , to a 

cd cn 


is less than or equal to 



2 


a 

cn 


2 

9 CT 

< ^ m 

= 4 * B 

m 


2 

cos 



( 122 ) 


at the minimum(C /$). Note that 2a m is the ratio of the received diffuse 
multipath power to the direct path power. An approximate calculation by 
Star as 1 ^indicates that for angles not too close to grazing, the scattering 
process itself discriminates against the multipath by approximately 2 dB. 
Allowing for an additional suppression of the multipath by the aircraft 
antenna of 6 dB, 2o 2 is -8 dB or 0. 16. In order that the variance of the 
diffuse multipath phase jitter component be less than that of the Gaussian 
noise, B must be greater than 


B > 
m 


(2a 2 ) 

m 


cos 2 /S = 56 Hz 


(123) 


B m has been estimated to be at least as large as this for grazing angles 
in excess of 10° and for aircraft speeds in excess of 600 mph. Con- 
sequently, the diffuse multipath contribution to the carrier loop phase 
jitter is on the same order as that of the noise at minimum(C /$>). There- 
fore, it doubles the effective noise power. At(C / $) ratios much higher 
than minimum, it may form the system noise floor. 


If the relative multipath Doppler frequency w ms lies within the 
carrier loop noise bandwidth, an upper bound on I ^cs (t > I is a ms* and an 
upper bound on the rms carrier phase jitter from specular multipath is 

Vv 2 a ms" 

6 

There is theoretical and experimental evidence to indicate 
that for a Gaussian distributed sea surface with an rms deviation of 10 cm, 
and for A - 18.75 cm, and y ^ 10°, the average relative specular multi- 
path power, is approximately 0. 065 if a 6 dB multipath suppression 

factor resulting from aircraft antenna directivity is included. 
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3.5.4 Effects of Multipath on Clock Phase Jitter 

The effects of multipath on clock phase jitter is now investigated. 
The video output, e^(t), 


e (t) sin/3 cos <f> (t)b (t) + n (t) 

v c 13 v 


+ (C) 1/2 sinfi a ms cosfd^t) - * c <t)]b 13 <t - r^) 

+ (C) 1 / 2 sinjS cosfd (t) - <f> (t)lx <t.)b (t -r ) 

L ms T c J c 13 ms 

-(C) 1 / 2 sin0 sinfd (t) - <j> (t)l x (t)b (t - T ) 

L ms T c 1 q 13 


(124) 


m s 


and the output of the clock loop bandpass filter becomes 
e Bp (t)=(C )^ 2 sin£ P 13 (^)cos ^Msin^t] + n Bp (t) 

+ (C ) 1/2 sin^ P 13 Aa ms c°s[d ms (t) - + c (t>] sinp^t - * c / ms ] 

+ <C > 1/2 sta ^13 ( f >COS [ d m S ,t) - ' “ c / ms l 


-(C ) ly/2 sin^ P 13 (-—) sin 


d (t) - 56 (t) x (t)sin 
ms r c q 


<*> .t - O) T 


q i d cf ms 

The phase jitter of the clock reference signal extracted 
by the clock-tracking loop is now 


(125) 


<l>ci (t) = k-o + Ko (t) + 


ci 


n 


cl 


ci 


(126) 


where <{> ^ (t) is assumed to be essentially unchanged from the situation 
of case (b f! This means that the clock loop behaves in a linear fashion 
and that effects of carrier lock error, <j> c (t), on the noise phase jitter of 
the detected clock signal <j> CJ g (t), are still assumed negligible. 


m 


The phase jitter component if> c ^ (t) results from the presence of 
the specular multipath clock component, and is approximately equal to 

<f> „ (t)»-|H (u ) | a cosfd (t) - <j) (t) - if) (u )] sinfu T ] 
c £ c£ ms ms L ms T r> T r-t 1 L r-c mc-i 


ci ms 


c £ ms J 


(127) 
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where |H c ^(u rnS )| and -|i c |(u ms ) are the amplitude and phase of the closed 
loop response of clock loop at the relative Doppler frequency of the specular 
multipath carrier component, ^ ms . Again, only the lowest order contri- 
bution to <b ^ (t) is considered. 

T ci 

S 

The diffuse multipath phase jitter component, ^^(t)* is deter- 
mined by noting that the diffuse multipath component of the "phase input" 
of the loop is approximately 



-cos[d (t) - <f> (t)lx (t)sinlw t 1 

L ms r c J c L ci ras J 

+ sin[d (t) - (t)lx (t)sin[w r 1 

L ms c J q L ci ms J 


(128) 


where the same type of approximation has been made here as was made 
for the carrier loop. The autocorrelation function of <Jr^(t), neglecting 
carrier lock error, <j> (t), is approximately 


(W* + T) > 


• 2 r 

sin [ 


W .T ] COSf 

ci ms J L 


Tl < j> (T) 

ms J m 


(129) 


and its power spectrum is 


(u) = sin 2 [w T ]^[4> (w 
id L ci ms J 2 L m 


-ID ) + 4* (<l) + to 

ms m ms 


)] 


(130) 


Again it is assumed that the power spectrum $ m (“) is spread out in fre- 
quency by an extent exceeding the clock loop noise bandwidth, B c ^, and 
4>.^(w) is assumed to be essentially constant over the clock loop noise 
bandwidth. The spectrum ^’ rn (w) is considered to be approximately Gaussian 

with bandwidth B such that 
m 


4? <w) 
m 


m 


-<f/B )‘ 

' m 


pr B 


m 


(131) 


The phase jitter 
variance, aP , 

ci d 


<J) i (t) is approximately Gaussian distributed with a 

approximated by (0) 2B So, from Eq. (130). 

id ci 
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CT 

m 



-(f / B Y 

ms m 
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ci 
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An upper bound on o ^ may be obtained by setting f^ ^ = 0, and 
replacing sm*[u r ] by l. c d 


ms 


2 , ^ m _ / 2 \ 

a c£ ' 7 7 r B c£ ( a ci 

d . m \ d 


(133) 


'max 


given by 


This is to be compared with a at minimum(C /$) which is 

n 


B 


ci 


C ^n p 2 u sinV— W 


a 


(134) 


With B c ^ selected, as before in Eq. (6) so as to make a 


inimum(C /$), B ci is 6.4 Hz, and the ratio of ^cl^^ax^ 0 CT cj£ iS 


0. 05 at the 


rm 


O’ 


c £ 


/ max 


2a B o 

m c£ 


2 Jn B 2 

r c^ m a ci 

n n 


(135) 


For the upper bound on the mean square clock phase error introduced by 
diffuse multipath, to be less than that introduced by the noise, 

at minimum(C / $), B m must be greater than approximately 12 Hz. 


c£ 


n 


The effects of the specular multipath on the clock loop are not 
likely to be as serious as with the carrier loop for two reasons. First, 
the factor sin[w CJ gT ms ] ^ ^c£ reduces its effect somewhat, particularly 
for w c ^ T mS near multiples of tt, and second, as a result of the smaller 
noise bandwidth of the clock loop, the relative multipath Doppler frequency 
go , is more likely to lie outside the clock loop noise bandwidth. Again, 
specular multipath should be almost entirely negligible except possibly 
for calm seas at low grazing angles. 


During the ranging measurement the B1NOR is replaced by just 
the square wave clock. This has the effect of setting p ^3 equal to 1 in 
Eq. (134). Equivalently, the result in Eq.(135) is modified by dividing the 
ratio of rn3: J’ 7 ‘^ii n ^13 to C'ktain ratio which is valid during 

ranging. It is seen therefore, that the diffuse multipath could very well 
be the dominant source of noise during ranging. This will be true unless 
B m i s greater than about 450 Hz, which appears to be unlikely. The rms 
ranging error in the presence of multipath could be significantly larger 
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than the error computed on the basis of Gaussian noise alone. More will 
be said about this in a separate section. 

3. 5. 5 Effects of Multipath on B1NQR Acquisition 


Proceeding to the question of BINOR acquisition probability, the 
output of the kth integrate and dump filter, assuming correct acquisition 
of all previous k - 1 square waves, is 


e k " 1 (s k + V 


(136) 


where s, and n, are now 
k k 


s k =f* T dt(C) ljf2 sin 0 cos < l ) c ^ t )( b 1 3 ® h LP ) c£ k f 1 “ 1 
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c SL 


( b 13 ® h Lp) t 
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, J ^ dt/n (§) h \ 

k j LP J t k y o » cl 


+ f tT dt(C)^ ^sin£ cos[d {t) - <b (t)]x (t) 
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Examination of n k reveals that n^ consists of a term arising 
from the presence of Gaussian noise and two terms arising from the 
presence of diffuse multipath. All of these terms are Gaussian distri- 
buted, and consequently may be treated as a single Gaussian random 
noise component whose variance is simply the sum of the variances of 
the separate components. 

Examination of s^ reveals that it consists of two components. 

The first is the desired signal component modified by the presence of 
carrier and clock phase jitter. If it were not for these phase jitters, 
this component of s k would be a deterministic constant. Instead, it is 
now a random variable. (Although not a Gaussian random variable. ) The 
second component of s^ arises from the presence of the specular multi- 
path component, and is also a random variable. It is random not only as 
a result of <f> c (t) and <() c ^(t), but also because d mS (t) is random (<|> mS is 
uniformly distributed from 0 to 2tt), as is a mg . 

For a given s , the probability of correct acquisition is simply 
K. 




(139) 


If the probability of obtaining each given were known, it would 
theoretically be possible to obtain the probability distribution of (V {s 
This would be a useful thing to know, since it contains complete 
statistical information. The statistics of s^ are extremely difficult to 
obtain, however, because of the way the random variables <j) c (t), <t> CJ £(t) 
and d mg (t) enter into s^. That is, statistics of their behavior over a 
period of time equal to -2T are needed. In addition, these variables enter 
into s k through the carrier jitter term cos <|> c (t), tiie clock jitter term 
c S. (t - t Q + (J> c (t) /co CJ g ) , and the relative carrier Doppler term coS [d ms (t) - 
(j> c (t)] , which further complicates matters. These functions multiply 
the filtered BINOR as well as each other (as indicated in Eq.(138)) and 
are then integrated over a time iT in order to form s^. Consequently, 
the effort involved in obtaining a generalized distribution of (P c { s g}jg 
would be considerable at the very least. Fortunately, however, much 
useful information can be obtained with much less effort by recognizing 
that on the average, and for the majority of the time, the effects of the 
perturbations (i. e. , the phase jitters, <f> c (t) and <j> c j?(t), and of the 
specular multipath) are small.' In fact they must be small if the system 
is to work properly. What is of interest primarily is the possibility of 
these perturbations causing the total probability of correct acquisition to 
drop from the required 9 9. 9% to, say, 98%, rather than from 99. 9% to 
say 80%. In other words, the following approach is justified; 1) assume 
that the effects of these perturbations are small, 2) develop an analysis 
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and determine answers on this basis, and 3) see if the results are 
consistent with this assumption when reasonable values for phase jitter 
variances, and multipath parameters are introduced. If the results are 
consistent with the assumptions then the results are valid. If not, then 
the system must be modified to reduce the perturbations and their effects 
so that they will indeed be small. With this in mind, (sj^is written as the 
sum of the average value of s^ plus the deviation from the average. 

< s k> + ( s k'< s k>) (140> 

The deviation from the average is assumed to be much less than the 
average, and is the random component of s^. Introducing Eq.(140) into 
Eq.(139) and expanding (P c {sjj)k i n a Taylor series about (s^) , the 
following approximation is obtained 



±i(_ 2_ -< S k > 2 / 2 <\>Y 2 3 k \( S k ~< S k>) 2 

' 2!2 r Km 

to second order. 

pThe average of s is 

K 



(141) 


(142) 


Here orders of statistical averaging and time integration have been 
interchanged, and (J) c and <f> c ^ have been assumed to be statistically 
independent. The specular multipath component of s^ averages to zero 
as a result of the uniform distribution of d ms . The lowest-order con- 
tribution of the specular multipath component of is discussed further 
below. 
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To a good approximation. 


E , (cos h ) is 
c 


E I (cos <J) ) » cos 
© c 

' n 



The expectation value 


( 143 ) 


E 


Ki\ 


(>>(■• - b) 


is computed by realizing that the results of Appendix C need only be 
modified by replacing ^ q by <f> ci - <|> ci , g and a ^ by a cje = < CT ci n +a c^ > 
and then averaging over <p c ^ . That is. 
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(144) 


For the noise and multipath levels which pertain to the NAV/SAT 
application, <j> c ^ /\[2 o is less than 1, so that Eq.(144) maybe expanded 
in a Taylor series about <|> c j) = 0 prior to averaging. This is derived in 
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Appendix D to second order in 



and is 
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where r<a j is defined by Eq.{106), and 
cjk 


i , v ,2 1 2 . 2 . , ^ /, 2 N 

H „(oj ) —a sm [w T 1 =E, | . 

cl! ms 2 ms *- cl ms-* <p ^ l r ci 

c S/ 


(146) 


Of special interest is the leading term in (P c { s k})k‘ This provides 
an estimate of the effects of the perturbations on the probability of correct 
acquisition. It is also the leading term in a computation of the expectation 
or average value of (P c { s k})k» and represents the lowest order effects of 
4> ^ci’ anc * the diffuse multipath. The lowest order effect of the specular 
multipath component of s^ is contained in the second order term of 
Eq. (141), and is 


i < s k> ,2 -< s k) 2 / 2 ( n k) y(\) ms \ 

? 7£>r e A<\>/ 


(147) 


where is the variance of the specular multipath component of s^. 

This variance is worst-case estimated in Appendix E, where it is shown 
to be given approximately by 


(s 2 ) ~ a 2 d 2 x 2 (t -r ) (148) 

\ k/ms 2 ms A k o ms 

The variance ^nj^, as mentioned previously, is equal to the sum 
of the variances of the component of involving n v (t) and the diffuse 
multipath component. The former is essentially the same as for case (b) 
since the only difference is in the phase jitter on ci^(t - t Q + <t> c l(t)/ w CJ f )> 
and this phase jitter was shown to have negligible effect in the computation 
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of the variance in Appendix A. Consequently, ^n^yfor case (c) is equal 
to 



(149) 


where is the variance of the diffuse multipath component of nj., . 

This is computed in Appendix F to be approximately 



2 (2<y2) 

2 m 

4 Jn 2B IT 
v m 



(150) 


where a m and B m were defined by Eq. (131) and the superscript "o" in 
Xk ^ r ms^ indicates that the lowpass video filter was neglected in the com- 
putation. 


To proceed further, note that the multipath delays are small 
compared to the periods of thex £(t) functions for the larger values of k. 
Also the functionsx £( T ) are constant and equal to 1 for small T for k 2 7. 
A worst-case estimate of <fn? V may be made by assuming(x ?(t ))^ = 1 

for all k and for all T (likely to be encountered in the NAV/SAT appli- 
cation). Then 


<%) d 


ci < ~JTm 


2 2 
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The total variance 




is then less than 
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2 2 
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k m 


(152) 


With SL = 4, (p/-/ 2 a) at minimum(C/$) was shown in Sec. 3. 3 (case (a)) to 
be equal to 2.840. With 2aj^ equal to 0. 16, and with yj^ 1, the variance 
of the diffuse multipath component will be less than the variance of the 
noise if 


B 

m 


2 a 


> 


/ P 

wF IT \ J 2 O 


14. 3 Hz 


(153) 


As mentioned, it is expected that B will be at least several times as 
large as this. 
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3,5.6 Computation of (P c )-p, the Total Expected Probability of Correct 
BINQR Acquisition 


The total expected probability of correct BINOR acquisition 
(P c ).p is simply given by the product of the 12 conditional probabilities 
(P c )k- The conditional probability, (P c )k* i- s the average value of 
(P c { s k})k> and from Eq. ( 141 ) through {152) is approximately given by 



and g k is 



(156) 


Worst-case approximations may be made as follows. In f and g , 

K K. 



(157) 

(158) 


( 159 ) 
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and Ei [X (t - i> ,/u which is given approximately by 

9 cj 2\ k ' ° c £ cS -n 

E, (x ( 1 oh IXo B ) 
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(161) 
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(163) 
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In the last term in Eq. (98), X , (t 

k o 


a cid l°ct ) ■ (° 

n'max ' 


2 . 2 

ci ) 
d n'max 




(<?Jo 2 ) 

\ cd cn/r 


The ratios 

max' 

are obtained from 


max 


from Eqs. (122), (135), and (151). 


Four situations are investigated below in which different B m are 
considered as well as the presence or absence of a small specular com- 
ponent. In each case, a 6 dB aircraft antenna suppression of the multipath 
is assumed. In estimating the relative power in the specular multipath 
component, aj^g, an rms water surface displacement, ct w , of 10 cm and 
a minimum grazing angle, y, of 10° was introduced into the equation 
(given in Beckman and Spizzichino, p. 246) which relates the average 
relative power of the specularly scattered component to ct , sin y, and 
the wavelength A . 

( 4t 7CT siny\2 
— ^ .) 

f = 0. 06457 (164) 

ms '4' v ’ 

The factor (^) multiplying the exponential arises from the assumed 6 dB 
antenna suppression factor. 


Again, as above in Sec. 3. 5. 3, a 2 dB suppression of the. .diffuse 
multipath component by the scattering process itself is assumed, so 
that with the antenna suppression factor of 6 dB, the relative power in 
the diffuse 2 multipath component, is -8 dB, or 0. 16. 

The four situations are, specifically: 


1) Aircraft flying at approximately 600 mph at low 

grazing angles. Assumed values for B m and a mg 

are B„ ^ 50 Hz and a__ 8=1 0. 06457. 
m ms 


2) Aircraft flying at about 2 000 mph at low grazing 
angles. Then, B m is sufficiently large so that the 
diffuse multipath interference has almost negligible 
effects. B m is considered infinite, and a^ g is 
still 0. 06457. 


3) Aircraft is flying at speeds exceeding 600 mph and 

for grazing angles in excess of 20°. B m is assumed 

to be approximately 200 Hz and a^ is assumed zero. 

w m s 

4) Same as 1) but with £ = 5 instead of JL =4. 
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The probability of correct BINOR acquisition is evaluated at 
the worst-case (C/$) = 35. 7 dB-Hz. Tables of the values of the various 
factors inEqs.(154) to (163) are presented in Appendix H. For convenience, 
Table H-5 showing the total probability of correct BINOR acquisition is 
repeated below. 


Table H-5 

TOTAL PROBABILITY OF CORRECT 
BINOR ACQUISITION (P 

c T 


Situ ation 

12 

<p c >t = TT <p c > k 

cl k = 1 c k 

1 

0. 98 0 or 98.0% 

2 

0. 996 or 99. 6% 

3 

0. 998 or 99.8% 

4 

0. 993 or 99. 3% 


The results of Table H-5, which give the total average proba- 
bility of correct BINOR acquisition for the four situations considered 
indicate that 

a) even for the worst-case conditions considered, 
the average probability of correct BINOR acquisi- 
tion for the most difficult situation (i. e. , situation 
number 1) is still fairly good, there being only a 
2% chance of not acquiring the BINOR correctly 
for jC - 4. 

b) for the same situation as in 1, but with a 25% 
increase in integration time (i. e. , situation 
number 4) the probability of not acquiring the 
BINOR correctly drops to less than 1%. 

c) for the situation in which the grazing angle exceeds 
20° and specular multipath may be ignored (situa- 
tion 3), the average probability of correct BINOR 
acquisition is 99.8%. This is the largest of the 
four situations considered, and barely falls short 
of the desired 99. 9%, and 
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d) for the case where only the specular multipath was 
considered (situation number 2), the value of (P c )-p 
became 99. 6%. As aircraft speed is increased 
the effects of the diffuse multipath become less 
important as a result of the additional frequency 
spreading. This situation then represents an 
upper limit on how good the average probability of 
correct BINOR acquisition can get as higher aircraft 
speeds are considered. 

Examination of the table of conditional probabilities (P c ) k / 

(Table H-4 in Appendix H) reveals that the values of (P c )j ( for the k = 1 
case is the lowest of all the k's. This is primarily a result of the effects 
of clock jitter as may be seen from the portion of Table H-2 showing the 
clock jitter signal suppression fact EJ, ^(xtUo " ^ci/^ci^* This situation is 
compensated for somewhat by the factor in Eq. (99) which was set 
equal to 1 for all k (and is approximately equal to 1 for k ^ 3), but which 
is equal to 0. 9363 for k = 1. In addition, the factor X k ^( T ms ) which was 
set equal to 1 in estimating i s more likely to be less than 1 . 

for the k = 1 case than for the larger k, since T ms is likely to be at least 
as large as several periods of the first divided-down square wave com- 
ponent of the BINOR. Consequently, X k ^( T ms ) (for k = 1) can be anywhere, 
between 0 and 1, depending on T ms . This means that (P ) k for k = 1 will 
more than likely be somewhat larger than that indicated in Table H-4. 


An additional worst-case possibility must be considered in 
which, for one or more of the acquisitions of the square -wave components 
of the BINOR, the specular multipath part of the signal component at the 
output of the k^h integrate and dump filter, s k , is at, or very close to, 
its maximum possible value for a given a mg , and is of opposite sign 
from the desired signal part of s k . 

The conditions which must be satisfied simultaneously in order 
for this to occur for the k^h acquisition may be seen, from Eq. (137), to 
be 


1) T mS must be either small compared to the period of 
the k^ 1 divided-down square wave component (most 
easily satisfied for the larger values of k), or a 
multiple of one -half this period (most easily satisfied 
for the smaller values of k). 

2) The relative carrier Doppler frequency, f ms = oo mg / 27 t, 
must be small compared to the inverse of the integrate 
and dump filter integration time, iT. And finally. 
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3) The initial value of the relative multipath Doppler 
phase, <J) , (see Eq. (110) must be very close to 

0° or isfff whichever value causes the sign of the 
specular multipath part of s^ to be opposite that of 
the desired signal part of s^. 


It is much more probable that this worst-case possibility occurs 
for the large values of k since at the low grazing angles for which a 
specular multipath component is likely to be present, t is indeed small 
compared to the periods of these square -wave components, thereby 
satisfying the first alternative of condition 1). For the low values of k, 
this alternative is no longer satisfied, and the special values of 
required to satisfy the second alternative of condition 1) occur only 
occasionally. Consequently, attention will be focussed only on the large 
values of k. 0 


To assure that conditions 2) and 3) are satisfied, it is sufficient 
to require that during the integration time, jfT, the relative carrier 
Doppler phase, d mg (t), remain within ±30° of the required 180°. In order 
to estimate the probability of satisfying this requirement, the ±3 0° range 
may be considered to be divided equally between requirements on the 
initial phase (| mg and on the relative Doppler frequency f ms - In particular, 
the requirements chosen are ^ ms = 180° ± 15°, and |f g iT| = 15°/ 360°. 
The probability that <j> ms lie within this range (assuming a uniformly 
distributed probability distribution for <J> mg ) is equal to the ratio of 30° 
to 360°, or 0.0825. The relative Doppler frequency | f m | must be less 
than about 0.8 Hz for i = 4 and T = (78. 1 25)~* sec. For aircraft speeds 
in excees of 6 00 mph and for a carrier frequency of 1. 6 GHz, it is fairly 
improbable that the relative Doppler frequency will be as low as this. 

This is seen to be particularly true when the effects of unavoidable vertical 
aircraft motions on the relative Doppler frequency are considered. This 
was developed further in Section 2. 


If the required conditions are fulfilled for the k^ 1 square -wave 
component, the probability of its correct acquisition, (P c )^, becomes, 
to a good approximation, 


( p >, * T 

c k 2 


1 + e 


rf /"nF— f (l - a )\ 
\ v 2 ct k\ ms// 


where f is as given in Eq. (155). 

K 


(165) 


It is of interest to see how the results for the total probability 
of correct BINOR acquisition (P c )j. are modified for situation(l) when 
one or more of the values of ( p c )^ for k = 7 through 12 of Table H-4 are 
replaced by the values obtained from Eq.(l6 5). The value of (P ) k 
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obtained from Eq. (109) with a ms = 0. 06457 and (d/V 2 a)f^ = 2. 363 is 


fP c>k ■ 2 


1 + erf 


(v 2 


- f. (1 - a ) 
<7 k ms ) 


= 0. 993660 


( 166 ) 


When this value is used to replace n of the values of (P c )^ for k = 7 - 12 
in Table H-4, the values of (P c )j given below in Table 5 are obtained. 


Table 5 


(P ) FOR SITUATION (1) VS. n 
C T 


n 

(P ) 

c 

0 

98. 00% 

1 

97. 52% 

2 

97. 03% 

3 

96. 55% 

4 

96. 08% 

5 

95. 6 0% 

6 

95. 13% 


Consequently, even in the rather unlikely event of this additional 
worst-case possibility occurring for as many as 6 square -wave acquisi- 
tions, a condition requiring an even smaller relative Doppler carrier 
frequency than the 0.8 Hz considered above, the total probability of cor- 
rect acquisition of the BINOR code is reduced only to roughly 95%. For 
the more likely case where 1 or, at most, 2 acquisitions have this 
additional worst-case condition, the total probability of correct BINOR 
acquisition is reduced only to approximately 97%. These values may be 
entirely acceptable in view of the fact that the conditions under which 
they occur happen relatively infrequently. 

3. 6 Ranging Errors Caused by Multipath 

When the BINOR code is replaced by the clock square wave after 
acquisition, the signal power in the clock loop increases by a factor of 
l/pj 2 = 19.65 or approximately 13 dB, In the absence of multipath, and 
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with the loop designed so as to maintain a constant noise bandwidth, the 
loop SNR increases by the same factor. However, if multipath is present, 
the variance of the multipath component of the clock phase jitter on the 
extracted clock signal increases by the factor 1 / Pj 3 , (as showninEq. (63)) 
and may now exceed the variance of the noise component of the clock phase 
jitter as discussed in Sec. 3. 5. 5. In other words, with multitap present, 
the rms phase error on the clock during ranging may be determined 
primarily by the multipath rather than by the front end noise. 

The total phase jitter on the extracted clock may be written, as 
in Eq. (126), as 


h* (t) ‘ hi (t) + hi (t> + hi <*> 

n s d 


(167) 


The variance of this phase jitter, a is 


2 2 2 2 

rr =0 +0" +0 

ci cl cl cl. 

n s d 


(168) 


The variance of the specular multipath component may be seen, from 
Eq. (127), to be equal to 


c i 


, TT , .,2 2 1 . 2 r , 

= H .(w ) a — sin ho t 1 

1 ' — - 0 L c £ ms J 


cl' ms' 


ms 2 


(169) 


s 

The variance of the diffuse multipath component is the same as given in 
Eq. (132) and is 


2 2 

cr „ = sin [u> .t 1 . „ 

c£ . L cl ms J y 7 T B 

d v m 


2cr 2 -(f /B ) 2 
m m s m „ 
e B 


ci 


(170) 


The variance of the noise component is now p-^ times what it was during 
BINOR acquisition (when it was 0. 05) or 


°ci < 171 > 

n 

2 

As a consequence of the factor sin [ w CJ g T r ns ] • there will be certain 
ranges of multipath delays for which sin^[u> CJ gT S J wil1 be close to zero 
and the multipath will be unimportant, and certain ranges for which it 
will be close to one and the multipath will be an important factor in 
ranging performance. Since worst-case possibilities are of primary 
concern, s i n 2 [ u c _g r ms ] is set equal to 1 . In addition, I H c ^(w mS ) | ^ and 
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2 . 2 

e xp(-f m /B ) are also set equal to 1. The first three aircraft flight 
situations discussed in the previous section are considered. The results 
are presented in Table 6. 


Table 6 


RMS CLOCK PHASE ERROR DURING RANGING 


c £ 


cl 


n 



s/max 



d/m axj 


l / 2 


Situation 

2 

°cl 

n 

2 

°cl 

s max 

2 

°cl 

d max 

a c 1 

1 

0. 002545 

0. 032285 

0. 001555 

0. 21537 

2 

0. 002545 

0. 032285 

0. 000000 

0. 18663 

3 

0. 002545 

0. 000000 

0. 002889 

0. 07372 


The values of total clock error a c ^ given in Table 6 should be 
compared with the value of , which is the value of cr^ in the absence 
of multipath and is equal to 0. $5 045. It is seen that for the third 
situation is about 40% higher than ct c ^ . This is the situation which per- 
tains to most aircraft most of the time. When specular multipath is 
present it may dominate over the noise and the diffuse multipath as indi- 
cated in situations 1 and 2. Fortunately the relative multipath carrier 
Doppler frequency for high speed aircraft is most likely to lie outside the 
clock loop noise bandwidth, and |H CJ g(co ms ) | 2 may be sufficiently less 
than 1 to suppress the specular multipath variance, , to the point 
where its contribution to a c £ is on the same order as crg^, . For a clock 
loop noise bandwidth, B Cjf , of 6.4 Hz, and for a loop damping coefficient of 
0.707, a relative carrier -doppler frequency of only 5 Hz will cause | H ,(w ) 

to drop to 0 . 1 . C mS 

The rms range error is related to rr through the relationship 

C X/ 



(172) 


where c is the velocity of light (~10® ft/ sec), and f CJ g is the clock fre- 
quency (320 kHz). In terms of rms range error, the results are as 
given in Table 7. 
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T able 7 


RMS RANGE ERROR 


Situation 

ct r (feet) 

no multipath 

25. 1 

1 

107. 1 

2 

92. 8 

3 

36. 7 


Several possibilities exist for reducing <j r . The simplest of 
these are 1) to increase the clock frequency, f CJ g', and 2) to reduce the 
clock loop noise bandwidth during ranging. The latter possibility reduces 
°c ^5 by reducing |H c ^(w mg )[^ for relative carrier Doppler frequencies 
which lie outside the now reduced clock loop noise bandwidth. 

3.7 , Summary and Conclusions 

From the analyses conducted in Secs. 3. 5 and 3. 6, the following 
conclusions may be drawn: 

1) The presence of a small multipath interference will not cause 
a drastic reduction in the total probability of correct BINOR acquisition if 

a) at least a 6 dB aircraft antenna discrimination 
against the multipath is maintained. 

b) aircraft speeds are on the order of 6 00 mph or 
more so that the frequency spreading of the diffuse 
multipath component is on the order of 50 Hz or 
more . 

c) the grazing angle exceeds 10° and the rms displace- 
ment of the sea surface exceeds 10 cm, so that the 
average relative specular multipath power is no 
larger than approximately 0. 065. 

For the vast majority of the time the grazing angle exceeds 20° 
and the aircraft antenna multipath discrimination exceeds 6 dB. The 
specular multipath component is then negligible, and the frequency 
spreading of the diffuse multipath component is more on the order of 
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100 to 200 Hz. The effects of the multipath in this case become almost 
entirely negligible during BINOR acquisition. 

For grazing angles close to 10°, specular multipath may become 
a problem for calm seas at 1.6 GHz. However, the probability of correct 
BINOR acquisition is still on the order of 98% on the average for SI - 4, 
and may be improved to 99. 3% by making SL = 5. Analyses of additional 
worst-case possibilities that may occur in the presence of specular mul- 
tipath indicate that it is very unlikely that the probability of correct 
BINOR acquisition (with A - 4) will fall below 95% at the worst-case SNR 
used in all of the calculations. 

During the ranging measurement, in the present system, the 
multipath becomes an important consideration. This is so because the 
BINOR is removed and only square -wave clock is transmitted. This 
increases the signal and multipath powers in the clock loop to the point 
where the multipath, rather than the noise, becomes the limiting factor. 
This appears to be especially true when specular multipath is present. 
However, the differential doppler offset between the direct-path and 
multipath carrier frequencies may be sufficiently high compared to the 
clock loop noise bandwidth to cause considerable suppression of the 
specular multipath clock component. This would substantially reduce 
the limiting effects of the specular multipath component on ranging 
accuracy. The diffuse multipath component causes an increase in rms 
range error of approximately 50% over that caused by front-end noise 
alone. The effect of the diffuse component is reduced only when the fading 
bandwidth is much smaller than the differential doppler which is unlikely 
in the present application. 
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4. MODULATION DESIGN 


4. 1 , Introduction 


The discussion begins with an investigation of the possible 
modulation methods by which the ambiguity resolving square waves 
may be applied to the 320 kHz subcarrier. Low deviation PM is 
selected. Then the optimum method for encoding these square 
waves is considered. The BINOR technique (hard limited sum of 
square waves) is rejected in favor of a simple sum of square waves. 
Optimum deviations for the square waves on the subcarrier and the 
subcarrier on to the carrier are selected. Then, a description of 
the ambiguity resolving procedure is presented. The discussion 
concludes with a hardware design exercise. This effort results in 
a detailed Block Diagram of a range data demodulation and measure- 
ment system appropriate to the optimized ranging modulation. 

4. 2 Subcarrier Modulation Methods 

Two forms of modulation are analyzed and evaluated. 

The first consists of a PM ranging subcarrier which is partially 
amplitude modulated by the ambiguity resolving signal. The second 
technique consists of a PM ranging subcarrier which is partially 
phase modulated by the ambiguity resolving signal. The resulting 
spectrum of each is derived and evaluated in terms of convenience 
and efficiency in the NAV/SAT system. 

4. 2. 1 Partial AM 


This specific form of modulation is described by 


e L (t) = <v/25 


sin 


to t 
c 


+ 



+ P^b(t)j sin to^t 


( 173 ) 
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where 


S = signal power 

- carrier frequency - 2irxl. 6 GHz 


CO 


subcarrier phase deviation 


co = subcarrier frequency = 2 ttx 320 kHz 
s c 


0 . 


ambiguity signal smplitude modulation deviation 


b(t) = ambiguity signal (sum of square waves or BINOR) 

The various components of the spectrum of ej(t) can be found by a 
simple trigonometric expansion. Thus, 


e 1 (t) . 

= cos V cos ( 6 i f 1 + p i b,t| ) sin “sc‘ ) 


+ cos io^t) sin ^6^1 + (3^b(t)j sin 


U sc t / (174) 


Equation (174) can be further expanded as 
e l (t) r 

, = sin (co t) cos (6, sin co t) cos (6, (3, b(t) sin co t) 

\] 2S c [ I sc t 1 sc 

- sin (5, sin co t) sin (6 , p , b (t) sin co t) 

1 sc 11 sc J 

+ cos (co t) sin (6 , sin co t) cos (6 . 6, b(t) sin co t) 
cl 1 sc 11 sc 

+ cos (6, sin co t) sin (6 . 3 , b(t) sin co t) 1 
1 sc 11 sc J 


(175) 
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The characteristics of b(t), the ambiguity resolving signal component, 
can now be exploited to continue the expansion. As noted above (b(t) 
is either a BINOR code or a sum of square waves. For the B1NOR 
case, b(t) = ±1. Then Eq. (175) is rewritten as 

.. . ^ 

e i (t> , 

7 == = sin (to t) cos (6, sin to t) cos (5, 6. sin to t) 

V 2S cl I sc 11 sc 

- b (t ) sin (6, sin to t) sin (6, 6 sin to t) 

1 SC 11 sc J 


cos (to t) sin (6, sin to t) cos (6, 0, sin to t) 
cl 1 sc 11 sc 

+ b(t) cos (6 , sin to t) sin (6. (3 . sin to ) 

1 sc 11 sctj 

(176) 

Each of the four terms of Eq. (176) can be expanded and each will 
contribute a distinct component to the signal spectrum. The first 
coefficient of sin co c t in Eq, (176) is 


K, 


cos (6, sin co t) cos (6,3, 
i sc 11 

= ~ [cos (6 1 (1 + P 1 ) sin CO 

+ cos (6 , (1 - (3 ) sin co t) ] 
1 1 sc J 


sin <o 


t) 

sc 


t) 

sc 


(177) 
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Employing the Bessel expansion of cos (Asin $) = J D (A) + 2 V 

2k 

(A) cos (2k) $ and sin (Asin $) = 2 Y i J2k+1 ^ S: ' J1 (2k+l) 4>, the 

2k+l 

coefficient Kj may be approximated by 


K l ” | [ J o( 6 l (1 + P l>) + J o( 6 l (1 - h’)l (178) 


Similarly, the second coefficient of sin u c t, K2, may be expanded 
and approximated by 


K 2 * I [ J o(\ (1 * V) - J o( 6 l (1 - V)j b,t) 

(179) 

In the same manner, the coefficients of cos co c t, K3 and K4 are, 
approximately, 

K 3 ' ( J l( 6 l “ + P 1 >) + J l( 6 l - V)] Sin “ sc ‘ 

(180) 

and 

K 4 " ( J l( 6 l d + Pit) - J l(«l - ?!>)] b <‘) 

(181) 
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This completes the expansion of (t) to first order. With the aid 
of Eqs. (178), (179), (180), and (181), the modulated signal is 
approximated by 


e L (t) 
\[2S 


sin (w t) 
c 


+ sin gj t 
c 


+ cos u 


k ( J o( h (1 + V) + J c( 6 i 11 - V] 
o( 6 l 11 +P 1>)- J o( 5 l (1 - V)l 

c‘l J i( 6 i (i + t) + J i( h (1 - Ml 


+ COS GO t 

c 


J 1 ( 6 1 + T ) + J 1 ( ' 6 1 f 1 - m ). 


sin gj t 
sc 


b (t) sin (jJ 


(182) 

The first term of Eq. (182) is the carrier component of 
the signal. The second term indicates an AM modulation of the 
carrier induced by the ambiguity resolving code, b(t). Thus, in the 
vicinity of the carrier, the spectrum of the signal is that of a random 
square wave amplitude modulated sinewave with maximum amplitude 
proportional to J Q (fij (1 + j3j) and minimum amplitude proportional 
to J Q ( 6 1 (1 - /?i))*. 

The third term of Eq. (18^ is the subcarrier component 
and the last term is the AM induced on the subcarrier by the code. 
Again, the spectrum in 


It is noted that since J Q (x) and (x) are not monotonic, the roles 
of maximum and minimum may be interchanged depending on the 
values of 6i and /3j_. 


t 

s c 
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the vicinity of the subcarriers is that of a random square wave amp- 
litude modulated sinusoid with maximum amplitude proportional to 
Jl (6 1 (1 + /3 1 )) and minimum proportional to (6i (1 - 0j)). The 
approximate spectrum of (t) is shown in. Fig. 13. 



Fig. 13 Approximate Spectrum of a Carrier, Phase Modulated by 

a Subcarrier which is Amplitude Modulated by a BINOR Code 


Derivation of a similar spectrum for the case in which 
b{t) is the sum of square waves is substantially more complicated. 
However, under the assumptions that Cl^jCl and |b(t) | max - 1, 
Eq. (175) maybe approximated by 


e L (t) 
\TIS 


sin (w t) ll - (6 , sin oj t) (6 P b(t) sin cj t)l 

cl 1 sc 11 sc J 

+ cos (to t) lb. sin OJ t +6. (3. b (t ) sin 05 t] (183) 
C l 1 sc it scj 


sin (to t) 


6 ^ 1 ] 

1 - — - — \+ cos(oj t)[6 sin to t+ 6 B b(t) 
2 J c L 1 sc 11 


sin oj t]+ terms at twice subcarrier 
sc J 
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Eq. (183) is similar in form to Eq. ( 183. Again, the spectrum in 
the vicinity of the carrier is that of a sinusoid which is amplitude 
modulated by the sum of square waves, b (t), with deviation 
Similarly, the spectrum in the vicinity of the subcarrier is, again, 
that of a sinwave, amplitude modulated by b(t) with deviation, 
Therefore, the spectrum shown in Fig. 13 for the BINOR case is 
also an adequate description of the signal spectrum when b(t) is the 
sum of square waves. The attributes of the partial AM spectrum 
will be discussed below relative to those of the partial PM spectrum 
which is derived next. 

4.2.2 Partial PM 


The partial PM signal is defined as follows: 


e, (t) = \T2S sin u t + 6 _ sin t + S b (t)) 

£ l c 2 sc 2 J 

(184) 

where 6^ = sub carrier phase deviations 

(3 2 - ambiguity signal phase modulation deviation 

Again, the carrier, subcarrier and ambiguity code components can 
be found by straightforward expansion of Eq. (184). 


e 2 ^ 

2 g- = sin (“ c t) cos^S 2 sin (oj^t + (3 2 b (t))\ 

cos oj^t sin^6 2 sin (u> t + p b(t))j 

(185) 

But, using the Bessel expansions, 

cos (6 sin (to t + (3 b (t) ) ) = J (SJ + ... 

' 2 sc 2 / o 2 

(186) 
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and 


(S 

sin 

(<jj t + 

\ sc 

p 2 b (t ijj 

= 2 J (6 ) sin 

("sc* + P 2 b(tl ) + 

= 

2 J 1 

(6 ) cos 

(e 2 b <t>) 

sin (co t) 
sc 

(187) 

+ 

2 J 1 

(6 ) sin 

(p 2 b(t)J 

cos (go t) + . . . 
s c 



Eq. (187) can be further simplified when b(t) is a BINOR code. Then, 
cos (£2 b(t)) = cos $2 and sin ($2 b(t)) = b(t) sin £ 2 - hi this case 
e 2 (t) is approximated by 


e 2 ft) 


J (5 .) sin (go t) 
o 2 c 


+ 2J, (8 _) cos S_ sin (go t) cos (go t) 
12 2 sc c 


(188) 


+ 2J, (6 _) sin (3 b(t) cos (co t) cos (00 t) 
122 sc c 


The first term of Eq.(188) is the carrier component. The second 
term is the subcarrier component and the third term is the sub- 
carrier's phase modulation due to the code. The approximate spec- 
trum of the partial PM signal is shown in Fig. 14. 



Fig. 14 Approximate Spectrum of a Carrier, Phase Modulated 

by a Subcarrier which is Phase Modulated by a BINOR Code 
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A similar spectrum is derived for the case in which b(t) 
is the sum of square waves. Appendix I shows that when b(t) is the 
sum of four unity amplitude square wave each with deviation 
jS 2 /4<l, such that | (3 2 b(t) | max = then 


cos 


^P 2 b(t}J = cos 4 


and 


P 2 3 8 2 

sin l p 2 b(t)j = b (t) sin — cos — 


;in b(t)j = 


(189) 


(190) 


Substituting the approximations of Eq. (189) and (190) into Eq. (187) 
results in an approximation for e 2 (t) given by 


^^2s 


J (6 ) sin (u t) 
o 2 c 


+ 2J, (6_) cos — — sin (go t) cos (go t) 
12 4 sc c 


P 2 3 8 2 

+ 2J. (6 _) b (t) sin — — COS - 7 1 COS (go t) cos (co t) 

1 2 4 4 Sc c 

(191) 

Eq. (191) has the same form as Eq. (188), therefore, the approximate 
spectrum of Fig. 14 is also valid for the case in which b(t) is the 
sum of square waves. 
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4. 2. 3 


Evaluation of Spectra 


The most obvious difference between the spectra of 
Figs. 13 and 14 is the presence of AM sideband around the 
carrier in the partial AM case of Fig, 13. These sidebands are 
undesirable and form the basis for selecting partial PM over partial 
AM. First, the presence of the sidebands limits the flexibility of 
the waveform. In particular, if the ranging modulation provides a 
clear spectrum in the vicinity of the carrier, this region can be 
used for data communication purposes in which straightforward data 
modulation techniques applied to the carrier. Alternatively, if the 
region around the carrier is kept clear, the problem of carrier 
acquisition is that much simpler. More important, the presence of 
significant AM sidebands around the carrier indicates that the ambi- 
guity resolving code power is distributed inefficiently in the signal 
spectrum. Consequently, both the AM sidebands about the sub- 
carrier and about the carrier must be demodulated and combined in 
the partial AM case to provide overall performance which is com- 
petitive with that of the partial PM signal structure. It is concluded 
that more receiver hardware is required by the partial AM signal 
to provide equivalent performance. 

4. 3 Ambiguity Resolution Techniques 

The BINOR coding scheme was evolved for cases where a 
binary (two-level) signal is required or when the system is peak de- 
viation limited. The BINOR lends itself to a relatively simple 
demodulator implementation when each of the ambiguity resolving 
tones is measured in sequence. However, full utilization of the 
data energy in the BINOR signal requires a rather complex imple- 
mentation because for full utilization of data energy these operations 
should be done in parallel. For example, consider a BINOR code 
consisting of thirteen tones, i. e. , clock plus twelve ambiguity 
resolving tones. The correlation coefficient of each tone in this 
code is pig - 0. 225. Now is a measure of the available 

power in each tone but it is the energy which determines the 
accuracy derivable from each. Thus, if each tone is measured in 
sequence for 1/13 of the total measurement time, then accuracy is 
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proportional to (1/13) p-^g . One method suggested for reducing 
this effect without vastly increasing complexity is to transmit the 
entire code in 6 groups consisting of clock + 2 ambiguity resolving 
tones. The tone correlation coefficient during each group is 
p - 1/2, thus, while each tone is still measured for 1 / 1 3 of the 
total time, the energy per measurement is substantially increased. 
This idea can be carried to the limit by simply transmitting each 
of the tones separately in sequence. In this way the nonlinear 
interaction among components and the subsequent spectral spreading 
of component energy is avoided entirely. The major disadvantage of 
sequencing to this extent is the added complexity of the transmitter. 

In addition, this type of implementation does not provide for failure 
during the sequence. That is, failure during one step of the sequence 
constitutes a total failure which cannot be compensated until the 
sequence repeats itself. Parallel mechanizations fail more grace- 
fully in that poor operation during part of the total measurement time 
does not necessarily imply total failure. 

Parallel operation of the BINOR technique with twelve 
ambiguity resolving tones requires twelve phase locked loop and 
twelve phase measurements. Clearly, this involves a substantial 
increase in hardware complexity at the receiver over that required 
for sequential operation. The amount of hardware required can be 
reduced if some of the tones are deleted from the code. For 
example, every other tone can be sent while still providing the same 
performance if the phase of each of these tones is measured with 
twice the precision. Thus, the number of PLL's and phase measure- 
ments required for ambiguity resolution is reduced by a factor of 2. 
Similarly, every third tone can be sent while maintaining perform- 
ance if each of these tones is measured with 4 times the precision. 
The number of PLL’s and phase measurement devices may thus be 
reduced by a factor of 3. 

On the basis of the above discussion it is decided that 
parallel operation with a reduced number of tones should be investi- 
gated further. The number of tones and the method by which the 
tones are combined, i. e. , BINOR or straight sum, depends on 
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mission requirements. These latter requirements determine whether 
or not the system will be peak deviation limited, which in turn deter- 
mines the utility of the BINOR code vs. the sum of square waves. 

4. 4 Deviation Requirements 

At this point, a modulation optimization is carried out in 
order to determine, roughly, the various deviations which will meet 
mission requirements. If it is found that a b(t) in the form of a sum 
of square wave can provide the necessary performance while still 
in the linear phase modulation region, i. e. , sin « <f>, then it will 
be concluded that the sum of square waves provides performance 
superior to that of the BINOR code because nonlinear interactions 
are avoided. 

The ground rules assumed in the optimization process 

are : 

1) The modulation is partial PM as defined above, with 
b(t) = sum of square waves 

2) Total time to acquire including range measurement 
but excluding carrier acquisition is 1 sec. 

3) Total signal power -to-noise density ratio is 
S/N q = 36 dB = 4000. 

4.4. 1 Signal-to-Noise Ratio Expressions 

Assuming a modulated signal of the form of Eq. (184), 
the signal-to-noise density available for carrier tracking is 


(S/N ) 
o c 


2S 

N 


(6 ? } 
O £ 


sm 


(w t) = 
c 


s j (6 J 

o c 
N 


(192) 
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Thus, the r. m. s. phase error in a carrier loop with double side 
noise bandwidth, B mc , is 


A 4> 

c 


1 


V 


2S 

N B 
o me 


J 2 (6 ) 

o 2 


(193) 


The demodulated tone power at the subcarrier freauency, 
w gc , is found from Eq. (191). After demodulation by the carrier 
reference, yif cos w c t, the subcarrier is 


e 

sub 


*/S 




COS 


sin 


(w t) 
sc 


(194) 


The noise density after demodulation is again N 0 ; thus, the 
amplitude of the reference accounts for coherent addition of signal 
and incoherent addition of noise after demodulation so that the signal- 
to-noise density is given by 


(S/N o>t = ^ J 1 (6 2’ COE 4 (195) 

o 

Thus, the phase error of a subcarrier tracking loop with noise 
bandwidth B nt is 




1 



( 6 . 


4 

cos 



(196) 
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Finally, the demodulated code power is found from 
Eq. (191) as follows: the code term of Eq. (191) is multiplied by the 
carrier reference, yielding 


6 sub 

code 


(t) 


n/ST 


2J (6 ) b(t) sin -p 
(12 4 


3 

cos 



(w t) 
sc 


(197) 


This, in turn, is multiplied by the subcarrier reference J~2 
cos (w gc t) which yield the video code signal 


'code' 1 ' = — ( 2J ! < 6 2» b <‘> T cos ~ ) (19S) 

Each of the square waves in the code is filtered to extract its 
fundamental sinusoidal components. Each of these is described by. 


e amb (t) 


4 \/2S 

IT 2 



(6 2 ) sin — 



sin w.t 
i 


(199) 

where co^ is an ambiguity resolution tone frequency. Therefore, 
the available signal-to-noise density for these tones is given by 


(S/N 


o a 


16_ _S 
2 N 

TT O 


(6 2 > sin 
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The phase error in each of the ambiguity resolving tone loops of noise 
bandwidth, B na , is, thus, given by 


A 

a 


Equations (192) through (201) provide the quantitative data 
required to determine the utility of the modulation scheme. But, 
before establishing deviation requirements, the optimum number of 
tones to be transmitted must be determined. 

4. 4. 2 Number of Tones 

This section analyzes the relative merits of transmitting 

4 tones or 6 tones. The 4 tones would be 40 kHz, 5 kHz, 625 Hz and 
78. 12 5 Hz, i. e. , the tone ratio is 8. For the 6 tone case the ratio 
is 4 between tone frequencies and the tones would be 80 kHz, 20 kHz, 

5 kHz, 1250 Hz, 312. 5 Hz and 78. 125 Hz. The analysis determines 
the deviations required to provide 99. 9 % probability of acquisition 
P ac q in 1 second given that S/N 0 = 4000. 

A detailed description of the acquisition process for the 
type of system considered here is given in Section 4. 5. However, 
for the purpose of this analysis is need only be noted that ambiguity 
resolution will be completed successfully if each of the tones is 
measured to a specific accuracy, 

Consider the 4 tone configuration. The maximum error 
permitted per tone is A (6 a max = ± 22. 5°. Thus, the probability 
of exceeding an error of 22. 5° on each individual tone P e ^ must be 
small enough so that the probability of acquisition exceeds 99. 9%. 
That is. 


/ 32 

S 

/ e 2 

Fl (6 ) sin — 


/ 2 

N B 

CCS J 

1 * 

o na 

\ 

/ 


( 201 ) 


P 

acq 


(1 



& .999 


( 202 ) 
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but P aC q is approximately 

P « 1 - 4 p 

acq ei (203) 

Combining Eqs. ( 202 ) and (203) yields 

P . £ . 00025 (204) 

ei 

Such a probability of error will be achieved if 22. 5° is 3. 66 times 
the r. m. s. phase error per tone, i. e. , 


„ 22.5 , n 

A s ~ degrees = 


radians (205) 


Substituting this in Eq. (2 01) yields that 


P 2 3 P 2 

J i (6 2 } sin "T cos T " 


8 


(8 x 3. 66) it 


- \[2QQQ 


(206) 


3.66 

45 


= .0813 


where, it is assumed that B na = 1 and S/N Q B na = 4000. Eq. (2 06) 
may be approximated for small $2/4 as 


J 1 (6 2 ) P 2 " ' 325 (207) 

An entirely similar computation may be carried out for 
the case of tones. Then the maximum permitted error is 45°, but 
6 tones must be measured with a total probability of success of 
99. 9%. That is 


P s 1 - 6 P . (208) 

acq ei 
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and 


P e . 000166 (209) 

This probability will be achieved if 45° is 3. 77 times the r. m. s. 
error, i, e. , 


£ r 77 degre ' s = Wltrf) radians (210) 

This result cannot be substituted in Eq. (2 01) because Eq. (2 01) is 
appropriate specifically to the 4 tone system. The corresponding 
error for a 6 tone system is derived along the lines of Appendix I 
to be 


A <j> 

a 



Now, substituting Eq. (210) in (211) yields that 


W Sln .~ cos ~T 


3. 77 
89. 5 


= .0422 


T ‘ - 


(4) (3. 77) tt 


( 212 ) 
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Again, assuming that $ 2/6 is small, then, Eq. (212) may be approx- 
imated by 

J (6 ) /9 £ . 253 (213) 

12 2 

Comparison of Eq. (213) with Eq. (2 07) shows that to a 
first order approximation the 6 tone system requires less peak devia- 
tion on the subcarrier (with given deviation 62 ) than the 4 tone system. 
It is further concluded that a more exact calculation will not alter 
this result. The deviation criterion is not the only optimization 
criterion. In fact, equipment complexity could be a more important 
factor if the deviation requirements of both Eq. (213) and Eq. (207) 
can be met easily. The deviation requirements are met easily if 
values of §2 and @2 can be found which simultaneously provide the 
required ranging accuracy (such as defined by Eq. (205) for the case 
of 4 tones) while satisfying Eq. (2 07) or (213). 

Assuming that 1) the desired accuracy in ranging is 3 0 ft. , 
i. e. , must be less than 60. 4 mr, 2 ) S/N 0 is 4000, and 3) the 

noise bandwidth of the range tone is arbitrarily chosen to be 8 cps, 
then, Eq. (196) may be recast as 


(6 2 ) cos 


_2 

4 


( 6 . 04 x 10 ) 


J 


(4) (4000) 

( 8 ) 


= . 368 
(214) 


A similar equation can be generated for the case of 6 ambiguity 
resolving tones, namely. 


6 P 2 

J, (5 J cos — — £ . 368 

12 0 


(215) 
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At this point, let us summarize the optimization conditions 
for the 4 and 6 tone system in TABLE 8. 

TABLE 8. Optimization Conditions for Ambiguity Resolving Tones 


Condition 

4 Tones 

6 Tones 

Signal power-to-noise density 

4000 

4000 

Amb. res. bandwidth, B rla 

1 cps 

1 cps 

Range tone bandwidth, B ^ 

8 cps 

8 cps 

Total prob of acquisition 

99. 9% 

99. 9% 

Prob of error per tone 

. 00025 

. 000166 

Max permitted error per tone 

22. 5° 

45° 

Deviation product for 

ambiguity resolution j ^ ( s ^ n 

P 2 3 ^2 

— — cos — — 
4 4 

^ 2 5 * 

(6 ) sin — cos ■ 

> .0813 


> . 0422 

Deviation product for j ^ ^ cos 
ranging accuracy 1 2 

> 0. 368 

4L 

4 

T IS 1 5 P 2 

J 1 (6 2 )cOS T 

> 0. 368 


If the deviations selected just meet the conditions of 
TABLE 8 , then for the 4 tone case 


ty sin t 

. 0813 


cos 


3 ^2 


4 ^2 

J i ( 6 z ) CQS T 

0. 368 


= 1 


or 


(216) 


tan 



0813 

368 


. 221 


(217) 
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and, thus, 


= 50° 


the subcarrier deviation required is found from 


J, (6 ) cos (12. 5°) = .368 


or 


J 1 16 2 > 


. 368 
. 908 


.405 


or, finally 


6 2 = 0.9 


The deviations required to meet the conditions of 
TABLE 9 with 6 tones are found similarly, namely 


^ 2 5 P 2 

J (6 ) sin — cos — 
l £ o b__ 

.0422 


T ,, 6 P 2 

J i (5 2 ) cos T 

0.368 


or 


tan 


. 0422 
0. 368 


= 0.115 


and, thus, for this case. 


P 2 = 39 = 


110 


(218) 


(219) 


( 220 ) 


( 221 ) 


( 222 ) 

1 


(223) 


(224) 



( 225 ) 


Again, the subcarrier deviation required is found from 

J x (6 ) cos 6 (6. 5°) = . 368 


or 


J 1 < 6 2 > 


. 368 
.9616 


. 383 


(226) 


or, finally, 


6 = .84 (227) 

The deviation requirements are summarized in 
TABLE 9. It is noted that the 6 tone system has slightly lower 
suppression effect on the carrier. All other performance factors 
are essentially the same. The slight improvement in performance 
with the 6 tone system is obtained at considerable expense in term 
of equipment complexity. In particular, two more phase locked 
trackers and phase measurements and two more ambiguity resolution 
computations are required by the 6 tone system relative to the 4 tone 
system. For this reason, the 4 tone ambiguity resolution technique 
is selected. 

TABLE 9, Deviation Requirements 


deviation, mod. parameter 

4 tones 

6 tones 

carrier dev. 6^ 

0. 9 

0. 84 

carrier amplitude, (6^) 

0. 8075 

0. 8312 

peak subcarrier dev. P ^ 

50° 

39° 

individual square wave amp 

12. 5° 

6. 5° 

ranging accuracy (1 second ave.) 

10 ft 

10 ft 

prob. of acquisition 

99. 9% 

99. 9% 
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4. 4. 3 


Comparison with BINOR 


The discussion above shows that the requirements 
imposed by characteristics of the mission are such that rather low 
deviations may be employed in the modulation. It was suggested 
in Sect. 4. 3 that under these conditions the BINOR modulation tech- 
nique would not be as desirable as a straight linear combination of 
ambiguity resolving tones. This section will attempt a comparison 
of the performance obtainable with a BINOR code versus the square 
wave performance derived above, 

First, the performance obtainable with subcarrier tech- 
nique, in general, can be compared with a normal implementation of 
the BINOR technique in which a BINOR code is applied directly to the 
carrier. Basically, the subcarrier technique provides complete 
acquisition of range data in about 1 sec (excluding carrier acquisition) 
while avoiding sequencing of the transmitted signal. The carrier 
component amplitude is constant and is suppressed roughly 2 dB by 
the presence of the ranging and ambiguity resolving signals. The 
ranging signal provides the required accuracy in an 8 cps bandwidth. 
Therefore, if the ranging data (which is available during the entire 
acquisition internal) is averaged over 1 sec an r. m. s. error due to 
noise of roughly 10 ft. can be achieved. The probability of correct 
ambiguity resolution is the same, under the same conditions, for 
both techniques. The excess of ranging accuracy may be traded for 
a shorter acquisition time or increased carrier power, as desired. 

By avoiding sequencing, the subcarrier technique offers important 
advantages with respect to spacecraft reliability and total "mean 
time to acquire" as discussed in Sec. 4. 3. However, an important 
advantage of the BINOR is the simplicity of receiver hardware with 
the normal implementation of the technique. The relative complexity 
of the receiver implementation required for the subcarrier tech- 
nique will be discussed in Sec. 4. 6 below. 

The BINOR combining technique may be used in an alternate 
way to provide the ambiguity resolving tones using the modulation 
waveform of Eq. (184). For example, a BINOR may be constructed 


\ 
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by hard limiting the sum of 5 square waves at 80 kHz, 5 kHz, 1250 Hz, 
and 312. 5 Hz. It is instructive to compare the performance achievable 
with this 5 tone BINOR vs. a 5 tone linear combination. 

The number of linear tones is reduced to 5 to provide a 
fair comparison with the BINOR which is typically realized with an 
odd number of components when equal weighting of all components 
is employed. Thus, both techniques considered will provide ambi- 
guity resolution down to 5 00 miles instead of the full 2 000 miles. 

Since only 5 tones instead of 6 are employed, the probability 
of correct acquisition is given by 1-5 {. 000166) = 99.917%. This 
modification allows for much of the previous analysis of the 6 tone 
linear sum case to be used here. 

The implementation for demodulation of the 5 tone BINOR 
consists of 5 PLL's which employ square wave reference signals at 
the appropriate frequencies. The demodulated video BINOR is dis- 
tributed to the inputs of these PLL's without any filtering. It has 
been shown that such an arrangement leads to degraded performance 
because the harmonics of the square wave reference convert more 
noise than useful error signal down to video. Nonetheless, the 
PLL's are assumed to operate in the normal manner for the purposes 
of this analysis. It is thus assumed that the components of the 
BINOR code which are available for tracking by the PLL's are des- 
cribed 


e (t) = 

bi 


\f23~ 


(2 J. 


< 6 2 > 


p_ sin p 


e . (t) 
sqi 


(228) 


where p^ = .375 is the BINOR correlation coefficient, ^ is the 
deviation of the subcarrier by the BINOR signal, and e S qj(t) is a 
square wave at frequency Equation (228) corresponds to 
Eq. (199) for the sum of 4 squares. 
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to be 


Equation (199) is modified for the case of 5 square waves 


amb 


(t) 


4 \f2S 
TT 2 



(6 ) sin 




sin (oj.t) 
l 


(229) 


To form a valid basis for comparison, the carrier devia- 
tion, 62 , must be the same in both cases. Furthermore, the 
suppression of the subcarrier by the ambiguity resolution code must 
be the same in both cases. This latter requirement leads to a speci- 
fication for fi-fo. That is, from Eq. (188) and (191) (modified for 
5 tones) 


cos 

b 



(230) 


Since the probability of acquisition has been appropriately modified, 
Eq. (223) can be rewritten for 5 tones as 


tan 



0. 115 


or = 32. 5°. 


Thus, from Eq. (23 0) 


(231) 


cos p, = . 968 
b 

p b - 14.5- 


(232) 
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Under the assumed conditions for comparison, examination 
of Eq. {22 8) and (229) shows that the BINOR technique will provide 
superior performance if 


P sin 6 2 

5 b 


4 

W z 



(233) 


substituting from Eq. (2 3 0) yields 


P _ tan 
5 


> 


4 

TT \/"2 



(234) 


Inserting numerical values 


P c tan = (.375) (.258) = .09675 

b b 

and 

4 P 2 

—V tan = (. 9) (. 115) = . 1035 

it's 2 5 


(235) 


The results of Eq. (233 show that the BINOR code is 
slightly inferior to the sum of square wave. Actually, the results 
are not as close as they appear to be due to the optimistic assump- 
tions concerning the BINOR implementation. It may be argued that 
since it is shown in Sec. 4. 4. 2 that a 4 tone sum-of-squares system 
is inferior to a 6 tone system it is possible that a 3 tone BINOR 
would be competitive with a truncated version of the 4 tone wave- 
form. This assumption is not quite true, as shown in Appendix J. 
The question, in any case, seems academic because a 3 tone system 
(BINOR or linear sum) with a ratio of 8 between tone frequencies 
will not provide the required maximum unambiguous range. 
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In summary, it is noted that this comparative analysis 
has shown that a BINOR code can be used as indicated in Eq. (188) 
to provide ambiguity resolution in a subcarrier modulation system. 

It is also found that the BINOR code will provide performance 
slightly inferior to that obtainable with a linear sum of square waves. 

4. 5 Ambiguity Resolution 

It was assumed in Sec. 4. 4. 2 that if the phase measure- 
ment error on each tone of the 4 tone system was less than 22. 5°, 
then ambiguity would be resolved correctly. In other words, when 
the ratio between tone frequencies is 8, then the measurement must 
provide accuracy better than ± ir/B. Similarly, for the 6 tone system 
in which the ratio between tones if 4, the maximum permitted error 
is ± 7r / 4. This section discusses the ambiguity resolution calculation. 
In particular, a ratio-of-8 system considered as shown in Fig. 15. 

For the purposes of illustration, the true phase in binary 
form is written above each waveform. Examination of these numbers 
shows that ambiguity resolution in the absence of noise is trivial. 
Denoting the higher frequency 7-bit phase measurement by 
A - a^aga5a^a3a2ai and the lower frequency 7-bit phase measurement 
by B - bYbgb^b^c^bp the ambiguity resolved phase measurement 
is just b7bgb5a7aga5a^a3a2a2. Thus, the precision of the answer 
is contributed by the higher of the two frequencies while the most 
significant bits come from the lower of the two frequencies. The 
next step in the ambiguity resolution process consists of similar 
operation in which the 10-bit number just derived is combined with 
the first 3 bits of the phase measurement of the next lower square 
wave, etc. 


The noise induced errors complicate the process because 
then the lower frequency phase measurement may be off by 22. 5°, 
i. e., $> £ = b7bgbgb 4 b3b2b 1 ± 0001. Clearly, an error of this mag- 
nitude can change the value of b 4 and, more importantly, b^. 
Therefore, in the presence of noise ambiguity, resolution proceeds 
as follows: the most significant bit in the higher tone, a7, is 
compared with b^ If they agree they answer is given by 



In other words, if a 7 - b^, then b^bgbg must be correct or the phase 
measurement error is in excess of 22. 5°. (The modulation has been 
designed so that this does not happen frequently. ) 

When a 7 does not equal b^, it is possible that an error 
of less than 22, 5° has been induced in the measurement of the 
lower tone phase. Since a .7 is the most significant number of a 7 -bit 
measurement it is assumed to be correct. 

Lack of agreement between b^ and a 7 does not necessarily 
mean that bg is incorrect. However, from Fig. 15 it is clear that 
errors of less than 22. 5° denoted as "Type I" will increment bg, 
while errors- such as that denoted by "Type III" do not affect bg. 

In order to determine which type of error has occurred, a differ- 
ence is formed as follows: 


A= SVsh 


' b 4 b 3 b 2 b l 


then if A > 0‘ and ag = 0, then a Type III error has occurred and 
then the correct phas is b 7 b 6 b 5 a 7 aga 5 a 4 a 3 a 2 a 1 . If A > 0, but 
ag = 1 , then b-ybgbg must be decremented, since a Type I error has 
occurred. That is, the correct phas is b 1 7 b' gb’ 
where b'^b'gb’^ = (b-jbgbg - 001). Similarly, when A< 0 and 
ag = 0, a Type II error has occurred and b^bgbg must be incremented 
to form b* 7b* gb* 5 = (b 7 bgbg + 001 ) and the correct phase is 
b* 7 b* gb* 5 a 7 aga 5 a 4 a 3 a 2 ap Finally, when A< 0 and ag = 1, a 
Type III error has occurred and the correct phas is formed as 
b 7 bgb 5 a 7 agaga 4 a 3 a 2 a^ as before. The algorithm of resolving 
phase ambiguities is shown in Fig. 16. Fig. 16 show's that 
the ambiguity resolution process continues after the first cycle 
described above as follows: the next highest tone phase measure- 
ment, say c = C7cg. . . cj is read in and compared with 
b 7 bgbga 7 ag. . . aj. In particular, C 4 is compared with the most 
significant bit of the previously derived partially resolved phase 
(this may be b 7 b ’ 7 or b-' ; 7 , as defined above). Then the algorithm 
is repeated, for example, if C 4 = b 7 , the further resolved result 
is c 7 cgCgb 7 bgbga 7 . . . a^ , etc. Thus, the loop of Fig. 16 is 
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Fig. 15. Ambiguity Resolution 



traversed four times for a 4 tone ambiguity resolution systems. A 
similar operation involving six ambiguity resolution computations 
is required for a 6 tone system. In that case the result after two 
loops of computation would be of the form CYCgb^bga^aga^ . . . aj, etc. 

Implementation of the ambiguity resolution algorithm of 
Fig. 16 is straightforward. Very few digital integrated circuits 
are required and the time required for the entire operation is 
negligible relative to the 1 sec acquisition time. 

4. 6 Phase Measurement Implementation 

A complete block diagram of a possible implementation 
for the proposed modulation is shown in Fig. 17. Discussion of the 
carrier loop and in particular the carrier acquisition subsystem is 
deferred to Section 5. It need only be said that the constraints or the 
design parameters of the carrier loop are essentially the same as 
those previously considered. The exceptions are that the carrier 
component is slightly suppressed by the presence of the data which is 
always "ON” and that during range data acquisition the carrier compo- 
nent is substantially bigger than that previously considered. 


Similar components are appropriate to the subcarrier 
tracking loop. But at this point, a technique is introduced to 
simplify the phase tracking measurement hardware required for the 
subcarrier component and the ambiguity resolution square waves. 

The implementation simplification is based on a device which com- 
bines the phase -locked loop operation and the phase measurement 
operation. The device is realized using digital modules as shown 
in Fig. 18. The operation is as follows: an input tone and a 
reference tone are supplied to the device. The reference tone is 
at a frequency which is 12 8 (2^) times higher than the input tone to be 
measured. The reference tone is divided down in frequency in a 
7-stage digital divider (counter). It is recognized that the contents 
of the divider is equal to the phase of the reference tone divided 
down to the input tone's frequency. Thus, when the input tone has 
a positive going zero crossing, the contents of the divider constitutes 
the phase difference between input and reference. For example, the 
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the contents of the divider would be 000000 at the positive going 
zero crossing of the input if the two signals were in phase. The 
phase difference between the two is read out of the divider by using 
the positive going zero crossing of input as a trigger for a buffer 
register. The contents of the buffer register drive a digital filter. 

This may consist of merely a scaler (i. e. , gain), thus simulating a 
first order phase-locked loop. More likely, the digital filter will 
consist of a scaler and an accumulator to simulate the lead network 
compensation of the second order phase -locked loop. 

The output of the digital filter is used to modify the count 
in the digital divider so that the input and reference are in phase, thus 
closing the loop, The count in the divider may be altered conveniently 
at the instant when the contents of the divider are zero. At that 
moment a trigger pulse is produced which causes the output number of 
the digital filter to be loaded into the counter. It is noted that the 
divider acts as a digital VCO, in this instance. For example, if the 
digital filter output is constant, then the contents of the divider are 
continually incremented. As a result, the unit recycles more rapidly 
and simulates a positive frequency offset. When this input is removed 
the phase accumulated by this process is retained and the unit continues 
to recycle at its nominal center frequency. 
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A second divide down counter is attached to the high frequency 
reference source. The contents of this counter are never modified and 
define the system reference phase. After a time sufficient for the 
digital phase locked loop just described to stabilize, the contents of the 
tracking loop divider are read out when the reference divider's con- 
tents are zero. This readout constitutes a filtered digital phase measure- 
ment of the input tone relative to the system reference phase. 

A receiver demodulation system employing 7~bit digital phase 
tracking and measuring device is shown in Fig. 19. The carrier 
tracking and ambiguity resolution computation portions of the receiver 
are deleted. 

The implementation shown in Fig. 19 starts with the output 
of the carrier loop coherent demodulator and provides as output the 
5 digital phase measurements {relative to the system reference phase) 
required for unambiguous range measurement. The five 7 -bit relative 
phase measurements may be obtained at a maximum rate 78. 125 meas. 
per sec. Specifically, these are read when the contents of the reference 
oscillator divider (right-hand side of Fig. 19 are zero. The phase 
measurements may, in fact, be read out at this high rate provided each 
is averaged in an accumulator over the 1 sec acquisition interval. Then 
at the end of the interval the data is loaded into the ambiguity resolution 
computer and the resolved range is computed. 

4. 7 Summary 

A ranging modulation has been designed for NAV/SAT appli- 
cations which provides the following characteristics: 

1) No sequencing is required at the satellite transmitter. 

All modulation functions are carried simultaneously. 

2) Signal format provides for convenient generation at the 
satellite. 

3) Ranging accuracy on the order of 10 ft. 

4) 99. 9% prob. of correct amb. resolution in 1 sec. 
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5) Carrier suppression of 2 dB by modulation. 

6) Carrier and clock jitter have negligible effect on 
ambiguity resolution process. 

7) Improved performance in the presence of diffuse 
multipath because all data is averaged over the entire 
acquisition interval. 

8) Receiver demodulator complexity comparable to that 
required for the BINOR (when range measurement 
phase meter and system reference generator are 
included in the comparison). 

9) Receiver hardware lends itself readily to medium 
or large-scale integration. 

10) The signal format can provide simultaneous data 

transmission capability employing straightforward 
communication transmission and reception concepts 
and implementations. 

These characteristics are desirable for the NAV/SAT 
application. It is therefore recommended that the modulation tech- 
nique be given further consideration. In particular, modulation 
deviations may be readjusted to reduce ranging accuracy (which is 
3 times better than required) in trade for shorter acquisition time 
and/or higher probability of acquisition and/or higher carrier compo- 
nent power. It is noted, however, that these alterations are not 
necessary for successful implementation of the modulation technique. 
Further analysis of the demodulation and measurement technique will 
yield improvements in implementation beyond those of Fig. 19. ' 
Again, it is noted that the complexity of implementation as it now 
stands in Fig. 19 is on a par with other proposed systems. It is 
further recommended that the communication capabilities of the signal 
format be investigated in more detail to determine the optimum modu- 
lation method and maximum data rate capability. 
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5. DIVIDE DOWN PLL ACQUISITION 


5 . 1 Introduction 


Carrier acquisition is an important consideration in the design 
of every PLL. There are four basic ways to handle the acquisition 
problem: design the PLL pull-in range to cover the entire acquisition 
range, employ discriminator aiding, sweep the VCO until the frequency/"'’ 
offset is within the PLL pull-in range, or employ a divide down acqui- 
sition aid. This section is devoted to a discussion of the divide down 
acquisition approach. Design parameters are derived and computer 
simulation results are presented. 

5. 2 Design Description 

The block diagram for a PLL employing divide down acqui- 
sition is shown in Fig. 2 0. The diagram shows two loops. The outer 
loop comprises the main PLL and the inner loop implements the divide 
down acquisition aid. The input signal, r(t), is mixed with the output 
of the loop VCO to produce a convenient IF frequency. The result of 
the mixing operation is filtered by a bandpass filter whose width is 
sufficient to encompass the entire carrier uncertainty region. It is 
assumed that the SNR in this bandwidth may be less than unity. The 
output of the bandpass filter is applied to the divide down acquisition 
loop limiter input and to the input of a narrow bandpass filter in the 
main PLL. When the acquisition aid has reduced the carrier uncer- 
tainty to within the bandwidth of the narrow bandpass filter, a carrier 
component signal will be applied to the main PLL limiter and the main 
PLL will acquire, disconnecting the acquisition aid loop. The narrow 
bandpass filter bandwidth is selected so that the SNR at its output is 
high. 

The discussion below is concerned with the performance of 
the acquisition aid loop before acquisition is completed. Fig. 21 
shows a block diagram of just the relevant components of the system, 
i. e. , the divide down acquisition loop. The rest of this section will be 
devoted to an analysis of this loop. The output of the wide bandpass 
filter is limited, frequency divided by N, and compared against a 
similarly divided reference in a linear phase detector. The output of 
the phase detector is filtered and applied to the input of the loop VCO. 
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It is assumed that the input signal, r(t), is the transmitted signal plus 
noise: 

r(t) = s(t) + n(t) (236) 

where n(t) is a sample function from a white, Gaussian process of 
spectral height No/2, doubled sided. 

A phase model of the loop is shown in Fig. 22. It is derived 
in the following manner. The transmitted signal is taken as 

s(t) = A cos (w Q t + 0( t)) (237) 

and 

r(t) = A cos (oJ Q t + 6 (t) ) + n(t) (238) 

The input signal, r(t), is mixed with the VCO output to derive the IF 
signal. Since a low-pass phase model is desired, all operations will 
be performed using the complex envelope components of^the IF signal. 
Thus, assuming that the output of the VCO is sin (u v t + 6 (t) ), the 
complex envelope components going into the bandpass filter are 

e = sin (0 (t) - 0 (t) ) + n (t) (239) 

q 1 

and 

A 

e . = cos (0 (t) - 0 (t) ) + n« (t) (240) 

l 4 

where nj(t) and n2(t) are independent sample functions from a white, 
gaussian process of spectral height N 0 /2A^. The two envelope com- 
ponents are each filtered with low-pass equivalents of the bandpass 
filter. The limiter, divider, and linear phase detector are modeled 
as an arctangent function, modulo 2N7T. The arctangent function has 
unity gain, whereas the divider, phase detector in the loop has a phase 
attenuation of N. This is accounted for by introducing an attenuation 
of N in the gain of the loop filter. The rest of the phase model is 
equivalent to the actual loop implementation. This low -pass phase 
model will be used in subsequent discussions. 

5 . 3 Acquisition Loop Design 

The next step is the derivation of criteria for the design of 
the acquisition PLL. In particular, the parameter N and the loop 
filter characteristic G(s) must be determined. As an aid in making the 
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mathematics tractable, a linear PLL model shown in Fig. 23 will be 
used. The limiter suppresses the signal when the input signal-to-noise 
ratio is less than unity. This introduces an attenuation inside the loop 
which is dependent upon the SNR. This attenuation will be designated k. 
In addition, the noise introduces cycle skips at the output of the limiter. 
Every time a cycle skip occurs, the phase undergoes a step change of 
2ir radians. The cycle skip events are Poisson distributed and equally 
likely, positive and negative. For the purpose of the linear model, an 
additional phase noise contribution due to cycle skips, CS(t), is intro- 
duced. CS(t) is a random walk, staircase function with step height 2i r. 

It is important to remember that the variance of such a function 
approaches infinity as time approaches infinity. Thus, it is important 
that acquisition occur before too many cycle skips occur. This criter- 
ion will establish the loop bandwidth. In particular, the characteristics 
of CS(t) will dictate the form for G(s). 


By inspection of Fig. 23, the transient frequency error (as 
observed at the VCO input) due to a cycle skip step is expressible as 


E f (s, - f 


K kG(s) 

v 

1 + K kG(s) 


v 


1 

s 


For the purpose of this discussion, it is assumed that G(s) will 
either result in a first order or second order PLL. If a second order 
PLL with half damping is desired (for k = 1) then 

to 

G(s) = 1 + — (241) 

s 

Since the bandwidth of the bandpass filter will be wide compared with 
the loop bandwidth, F(s) will be ignored. Figure 24 shows plots of the 
frequency error in the loop for a single cycle skip step, parametric in 
the suppression attenuation factor k. Note that for small k, disastrous 
errors can' result. This particular characteristic makes a second order 
loop unacceptable, since, in fact, k will be minimum during maximum 
cycle skip activity. 


The first order implementation is considered next. For a 
first order loop, the frequency error produced by a cycle skip never 
exceeds w n , and the damping is a function of the suppression factor k. 
This is definitely superior to the performance of a second order loop. 
But, a first order loop has a phase error proportional to the carrier 
frequency offset. If the loop bandwidth is relatively narrow, the 
resulting phase error can be hundreds of degrees. To avoid difficulty 
at large offset frequencies, the divider ratio N is selected so that the 
total phase error at the phase detector never exceeds the linear range 
of the phase detector which is ± 180°. To permit ample safety margin, 
the loop should be designed not to exceed ±9 0°. 
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Summarizing, G{s) is determined to be a pure gain, thus 
producing a first order PLL. The gain will be K v k and the loop band- 
width will be K v /2. The maximum carrier frequency offset in Hertz is 
equivalent to the single-sided bandwidth of the bandpass filter, BW / 2. 
Thus, the maximum phase error in the loop will be 

Phase Error = 7rBW/K k radians (242) 

v 

The divider ratio N should be selected with a significant safety margin. 
In particular. 


N ^ 2 (Phase Error)/ir (243) 

where N is assumed to be an integer power of 2, K v k is selected so 
that a small number of cycle skips (less than 10) occur within the 
expected loop time constant, l/K v k. The suppression factor k and 
the average rate of cycle skips is determined from the worst case 
expected SNR. 


.5. 4 Computer Simulation 


As a means of verifying the expected performance of the divide - 
down acquisition technique, a digital computer simulation was per- 
formed. The simulation included one additional effect not previously 
discussed. In particular, it was assumed that the center frequency of 
the bandpass filter was offset from the nominal IF frequency. This 
would introduce a frequency offset noise bias. The block diagram used 
for the simulation, including this offset effect, is shown in Fig. 2 5. 

The frequency offset is simulated low -pass. This is accomplished by 
performing a frequency offset operation on e q and e-, to produce the 
offset quantities e q ' and .e^ . 


e T = e^ cos t + e- sin a> 0 ^ t 


(244) 


ej' = e i cos uj 0f t - e q sin Q 0f t (245) 

where co 0 ^ is the offset in the filter. After the filtering operation, the 
offset is removed by performing a complementary operation. These 
operations are depicted on the figure. The simulation program is 
listed in Appendix K. It is a straightforward implementation of the 
operations shown in Fig. 25. 
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Fig. 25 Phase-Lock Loop Simulation Block Diagram 











The first simulation performed was to verify the expected 
failure of the second order loop. The bandpass filter implemented 
was a single pole filter with noise bandwidth of 2 0 kHz, double sided. 
The divider ratio, N, was 16 and the double -sided loop bandwidth was 
500 Hz. The loop would not acquire for signal-to-noise ratios less 
than unity in a 20 kHz bandwidth. 

The second simulation performed was to verify the expected 
performance of the first order loop. At each noise level several 
different noise sample functions were employed. This was accomplished 
by initializing the random number generator with different starting 
numbers. Acquisition was achieved within less than 10 milliseconds 
for signal-to-noise ratios down to -10 dB in the 2 0 kHz IF bandpass 
filter. It was also found that the loop would acquire for signal-to- 
noise ratios down to -30 dB. This excellent performance is caused 
by the nature of the first order loop. The loop always tends to reduce 
the frequency error and eventually, the random noise, in combination 
with the action of the loop, produces a condition where the IF signal 
lies within the bandwidth of the narrow tracking loop bandpass filter. 

The hand-off to track is then immediately accomplished. For all 
simulation runs, the loop bandwidth was set to 1000 radians/sec. This 
is consistent with a reasonable narrow bandpass filter design for the 
track loop. 
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6. STEERABLE NULL CONTROL FOR THE 
DIFFUSE MULTIPATH CASE 


6 . 1 Introduction 

In the Interim Report , ADCOM demonstrated how specular 
multipath interference in an aircraft -navigational satellite link could 
be reduced by the use of a steerable null control system. The major 
results of that problem are reviewed briefly in the following section, 
and the modifications required for the inclusion of both coherently 
fading and incoherently fading diffuse multipath are indicated. ADCOM' S 
conclusions concerning the effectiveness of the steerable null control 
system for combatting multipath interference at 1. 6 GHz in a NAV/SAT 
application are summarized as follows: 

1) Since diffuse multipath interference arises from scattering 
off a much larger area of tne sea surface than for the case of specular 
scattering, the multipath interference in the diffuse case arrives at 
the aircraft from a range of different directions centered on the spec- 
ular direction. Consequently, placing a sharp null in the direction of 
the specular point will not eliminate all of the multipath interference. 

In addition, slight variations of antenna gains with directions of 
arrival of the diffuse multipath components will preclude the possibility 
of achieving the required amplitude gain balancing necessary for 
achieving a true null. At best, only the coherently fading diffuse mul- 
tipath component can be nulled, with the incoherently fading diffuse 
multipath interference being simply reduced somewhat. 

2) The presence of incoherently fading diffuse multipath 
interference interferes with proper operation of the steerable null 
control system and may prevent the null from being placed on the 
direction of arrival of the coherently fading diffuse multipath component. 


Study of a Multipath Rejection Technique applied to Aircraft 
Navigation by Satellite, 15 April 1968. 
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3) Steerable null control at 1. 6 GHz does not work as well 
at low grazing angles as at intermediate grazing angles (i. e, , around 
3 0°), even for the case of specular scattering. For the 1. 6 GHz 
NAV/SAT channel grazing angles as low as 10° are being considered. 

4) Steerable null control may not be needed at 1. 6 GHz, since 
antennas can be built at these frequencies which may adequately dis- 
criminate against the multipath if the angular separation between the 
directions of arrival of the direct path and multipath signals at the 
aircraft is greater than about 2 0°. 

A modification of the steerable null simulator (constructed 
previously under this contract) to simulate the randomly fluctuating 
carrier phase of a coherently fading diffuse multipath component was 
carried out. The results indicated that the steerable null system could 
work well under such circumstances in that the rapidly fluctuating 
multipath carrier phase prevented the phase of the multipath carrier 
from dwelling too long in a region where operation of the steerable 
null control system became ineffective. 

6. 2 Review of Steerable Null Control for the Case 

of Specular Multipath Interferences 

Specular multipath interference implies a mirror-like reflec- 
tion of the transmitted signal off the sea surface. Such is the case for 
VHF frequencies at grazing angles below about 30° and for normal seas. 

In the Interim Report, a direct-path signal of the form 
e d <t) = cos^t + 0(t)j (246) 

was considered. The received multipath signal was written as 

e (t) = a cos /u t - w r , + 0 (t - T H )) (247) 

m \ c c d d / 

for the case of pure specular scattering, and as 

e (t) ~ a(t) cos fcj t + 0 (t-r ,) + $ (t ) ) (24 8) 

. m \ c; d c / 

for pure diffuse scattering. Here o> c is the direct -path carrier fre- 
quency, including Doppler, and r d is the relative multipath time delay, 
and is time varying. The signal in Eq.(248) is actually a coherently 
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fading diffuse multipath component in that only a single relative mul- 
tipath time delay was considered. If the range of relative multipath 
delays received at any given instant is a significant fraction of a period 
of the highest modulating frequency in 9 (t), it may no longer be appro- 
priate to consider only a single modulation delay. In this case, the 
received multipath interference may be represented by the sum of a 
coherently fading diffuse multipath component, as in Eq. (248), and an 
incoherently fading diffuse multipath component involving the sum over 
a range of relative multipath time delays. The received multipath 
signal then may be written as 

e (t) = a(t)cos/cj t + 0(t-T.) + (t)) 

m \ c d c / 


(249) 


+ V a^(t) cos (« c t + 6 (t - r^) + ^(t)) 


In the steerable null system two essentially omnidirectional 
antennas, which are vertically separated by one-half wavelength, are 
used to receive the direct path and multipath signals. The separation 
of the two antennas results in phase shifts between the carrier-phases 
of the received direct path and multipath signals in the upper and lower 
antennas. An additional controlled phase shift between the signals 
received by the upper and lower antenna is introduced at IF so as to 
result in the multipath components of the signals received by the two 
antennas being 180° out of phase, and the direct path components 
being approximately in phase. Summation of the signals from the 
upper and lower antennas at this point can result, if the multipath 
component amplitudes are equalized, in a nulling of the multipath 
components and an enhancement of the direct path components. The 
purpose of the steerable null control system is to electronically main- 
tain the proper phasing and amplitude balance required for multipath 
null as the aircraft flies along its course and the directions of arrival 
of the multipath and direct path signals change. 


To do this, use is made of the fact that the presence of a 
residual multipath component at the output of the summer will result 
in amplitude as well as phase fluctuations in the resultant signal. If 
these amplitude fluctuations are detected by a square law detector 
followed by a high-pass filter, a signal fluctuating at the modulation 
rate and containing two components results. One of these components 
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is proportional to the phasing imbalance and the other to the amplitude 
imbalance. Since the fluctuating signals multiplying the two error 
signals are approximately orthogonal under certain circumstances, 
extraction of the error signals with negligible cross-talk error can be 
achieved if the fluctuating signals multiplying the error signals could 
be generated from the signals at the inputs to the summer and used as 
reference signals in a coherent detection scheme. Appropriate refer- 
ence signal extractors for obtaining these fluctuating signals can pro- 
duce good reference signals for the case where the errors are already 
small. 


An experimental simulation of the multipath signals, the 
antennas, and the steerable null control system for phasing control 
was built which demonstrated the basic operation of the steerable null 
system. 


The pertinent equations for description of the operation of the 
phase control portion of the steerable null system are now reviewed in 
order that the modifications required for the case where incoherently 
fading diffuse multipath components are present can be more easily 
understood. The multipath signal is assumed to be a coherently fading 
diffuse multipath interference. 

The signals at the inputs to the summer are written as indicated 
below. The subscripts "u" and " refer to M upper" and ’’lower" 
antenna received signals, while "d" and "m" refer to the "direct path" 
and "multipath" components. 

Upper antenna direct path signal component: 

e ^(t) = cos yw^t + 6 (t)j (250) 

Lower antenna direct path signal component: 

(t) = (1 + r\ ) cos L t - £ . , + 6 + 9 (t)) (251) 

d \ o id j 
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Lower antenna multipath signal component: 

e. (t) = (1 + T ) a (t) cos t - <{> (t) + 6 + 0 (t - t )) 
i-m m \ o c d / 

Upper antenna multipath signal component: (-52) 

e (t) = a(t) cos (u t - ^ (t) - |3 + 0 (t - r,)l 

um V o c um d ) (253) 

Here I'd and r m are amplitude imbalances in the direct path and multi- 
path components respectively, and (3 um are the carrier phase shifts 
introduced by the antennas, and 6 is the carrier phase shift introduced 
by the electronic phase shifter. The phase shift 6 is written as the 
sum of the value required for null, 5 n , and an error component, . 

The value of 5 n may readily be seen to be 

■ ^n “ 71 ^um (254) 

Defining the carrier phase difference between the direct 
path components at null, y d , by 


7d = 6 n ’ £ld 


(255) 


The output of the summer becomes 


e (t) = e d (t) + e (t) 

sum sum, sum, m 


= cos ^ uo t + </> o + 6 (t)j 

- j r m | a(t) cos too t + <f>^ + 6 (t) - (£(t) + A 6 (t, r d ) )j 

V ° ° 7 (256) 

+ ] 6 e | a(t) sin^co 0 t + <f> 0 + 0(t) - (£(t) + £8 (t, r d ))j 
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where 


A d = l 

| (2 + r d ) 2 cos 2 [f (y d +6 £ )| + r d 2 sin 2 |^(y d+ 6 E )j } 

(257) 

4 

O 

T (^d + 6 f> + ?d 

(258) 

tan C 

p 

d = (2 + r d ) tan 1 2 (y d + 6 c’l 

(259) 

( r mj 

p cos <-j- if) r m 

(260) 

{•*} 

= (2 + T ) sin( 6 f ) 

m 2 e 

(261) 

C(t) ■ 

- — y, + $ (t) + B + C 

2 c ' p um M 

(262) 

A 8 (t, - 

r d ) = 0 (t) - 0(t - Td ) 

(263) 


These results are obtained by straightforward vector addition and 
trigonometric manipulation. 

If this signal is square -law detected and high-pass filtered, 
the resulting signal is 

e iAM (t)= -{ r m } a(t) A d (cos[f<t) + A9(t, Td ,]) (264) 

HP 

6 e [a(t) A^sin (t)+ A 9 (t,r d )]j 

HP 

where the subscript "HP 1 ' denotes high-pass filtering of the bracketed 
expression. If the delay 73 is large enough to produce a sufficiently 
large fluctuating signal A 0 (t, T d ), the fluctuating signals 


142 



j^cos £(t) + A 9 (t, r d )H and 



will be approximately orthogonal. This will be particularly true if 
0 (t) is a narrowband signal process, such as a cluster of range tones. 

In this case, generation of signals proportional to these fluctuating 
signals can permit coherent detection techniques to be used to recover 
the error signals j r m j and j 6 e j. Appropriate signal processing of 
the input signals to the summer (as described in the Interim Report) 
result in the generation of a phase error reference signal e§(t) given by 


e_(t) = - 
o 


a(t) A, 


(sin[(0 £ 


c d ) + at) + Ad (t, r d ) 


HP 


where 


a = (2 + r + r ) cos 

5 dm 




1 + tan 6 



(265) 


(266) 


tan fig 


( r d - r m ) 

( 2 + r d + r m) 


tan 



y d + h 


(267) 


and tan£ d was given in Eq. (14), Note that as T m "* 0, tan0^ 
approaches TanC d and (^5 - Cd) approach zero. The quantity (^ - Cd) 
is also reduced as 6 f - 0 and 1/2 y d becomes small. That is, if the 
amplitude and phase imbalances are small, and the direct path signals 
add approximately in phase at the summer, the phase ({6 g - Cd) will 
become small and the phase error reference signal will be precisely 
what is required for coherent detection of the phase error signal | 6 e j. 
The steerable null system works best therefore when the errors are 
already small. 


6 . 3 Modifications in Steerable Null Operation for the 
Case of Incoherently Fading Diffuse Multipath 

For the case where the multipath interference contains sig- 
nificant incoherently fading diffuse multipath components, the multi - 
path components at the input to the summer must be modified to 
include both the range of relative multipath time delays that may be 
present and the fact that even for components having the same relative 
multipath time delays different directions of arrival may produce 
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different antenna phase shifts and gain imbalances. In this case then 
Eq. (252) becomes 

e. (t) - (1 + r ) a(t) cos U t - (t) + 6+ 0 (t - r d )) 

£m m \ o o a / 


+ E a kj (t) (1 + r k .) cos (u Q t - + 6 + 0 (t - T k )j (268) 

k. J 


e (t) = 


um 


= a(t) cos /w t - (t) - 0 + 0 (t - r , ) ^ 

\ o o um d J 


+ J a k . (t) cos (o> Q t - ^ kj (t) - + 0 (t - r k )) ' (269) 

k, j 

If A/3 k j is the additional carrier phase shift caused by the antenna from 
the kj incoherently fading diffuse multipath component, i. e. , 


AS . = 8 - 0 

kj kj um 

Then the output of the summer is modified to 

e (t) = A cos (u, t + (6 + 6 (t)\ 
sum d o o j 


-l av 

k, j k j 


. (t) 


+ E a ki w 

k.j kj 


(270) 


-a(t){r m [cos |E Q t + + 0 (t) - (f(t)+ A 0 (t, r d ))) 

+ a(t)| 6 C |sin (u Q t + ^ + 0 (t) - ( £(t) + A0 (t, r d ))j 

cos (4" 6f + T" A Aj) r kj ] cos (v + K * 6 (t) ~ (271) 

(W*) + Ae<t ' T k ))) 

<2 + r k3 )sin (T^ + T A A 0 )] sin (“o t+ C + S(t) - 

(« kj <t> + A«(t,T k >)) 
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where 

5 kj = T y d + ^ka <t)+ ‘T A V *um + «d 12 ' ™ 

and 

A 0 (t, r k ) = 0 (t) - 6 (t-r k ) (273) 


The terms in square brackets in Eq. (271) may be expanded as 



+ ( 2+r k 3 ) COS (T 6 «) Sin ("2 _A 3 kj) 

(275) 


With and I\j small for all kj, the second term in Eq.(274) may 

be neglected compared to the first term. These terms may be written 
in abbreviated form as 

cos (“2~ 6 e + i" A ^kj) r kj ” J r kj j 

and 

2 + r kj) 3in (i-h + -r A %)] = [ 4 *kJ + | A \j 



(276) 


(277) 
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With the above notation the output of the square-law detector and high- 
pass filter combination becomes to first approximation 



The approximation involves neglecting products of incoherently fading 
diffuse multipath components with each other and with the coherently 
fading multipath components. In the above note that when and 
r m - 0, j l 6 ekj [goes to zer °, but j r^j | and j A/? k -j j do not. These rep- 
resent the amplitude and phase imbalances of the kj incoherently fading 
diffuse multipath component. Also note that the presence of these 
components will prevent 6^ and T m from going to zero when the loop 
is closed. Since the rms value of any a k j (t) is much less than that of 
a(t), the coherently fading diffuse multipath component should dominate, 
with the incoherently fading diffuse multipath components acting more 
or less like noise. 

As might be expected, the phase error reference signal is 
also corrupted by the presence of incoherently fading diffuse multi- 
path components. This may cause the system to fail more readily than 
the presence of these components in the output of the square -law 
detector high-pass filter combination in that a good reference signal is 
necessary for proper system operation. The phase error reference 
signal is to a first approximation 
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e 6 (t) = 


where 


and 


- a (t) f sin 


“A \j <‘> A k ^sin 


(V T c) +f (‘> + a «('- T d ) 

(^- T d) + «ki' t,+ A# (‘- T k -) 


( 279 ) 


kj 


'd kj, 


1 cos 

-y—y, + g 

L , + 2, | 

) 

2 d e 

1 k J| 


1/2 


( 280 ) 


tan * kj = (2+r +r.) tan 

d k] 


-r- y. +6 
2 d c 


( 281 ) 


Again products of incoherently fading diffuse multipath components with 
each other and with the coherently fading multipath components have 
been neglected. 


If most of the received multipath power is in the coherently 
fading diffuse multipath component, the steerable null control system 
may function to place the null in the direction of arrival of the coher- 
ently fading diffuse multipath component. On the other hand, if the 
incoherently fading diffuse multipath component dominates the coher- 
ently fading multipath component, the steerable null control system 
will probably not function properly. The non-linear nature of the 
steerable null control system and the complex nature of the various 
signals precludes an exact analysis of the behavior of the steerable 
null control system in the presence of large incoherently fading diffuse 
multipath interference in much the same way as a phase -locked loop 
becomes difficult to analyze near and below threshold. However, since 
incoherently fading diffuse multipath interference is much more noise - 
like than specular or coherently fading diffuse multipath insofar as 
its effects on range tone phase errors are concerned, an antenna which 
suppresses the multipath sufficiently so that only first order interfer- 
ence effects are significant, may be all that is needed to insure 
acceptable ranging performance. 
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In summary, if the incoherently fading diffuse multipath com- 
ponent can be made smaller than the direct path signal by the use of a 
single directive antenna, that is all that may be needed for acceptable 
ranging performance. Whereas for proper steerable null control sys- 
tem performance, the incoherently fading diffuse multipath component 
must be smaller than the coherently fading diffuse multipath component. 
This is so because in this control system the coherently fading diffuse 
multipath component plays the role of the desired signal and the 
incoherently fading diffuse multipath signal plays the role of interfer- 
ence. The expression for the phase error reference signal e§(t), 
demonstrates this point quite nicely. 

6. 4 Simulation Experiment for the Steerable Null Phase 

Control System for the Case of Coherently Fading 

Diffuse Multipath 

The steerable null phase control system simulator described 
in the Interim Report was modified to simulate one of the major 
effects of a coherently fading diffuse multipath interference. As indi- 
cated previously, the difference between the specular and coherently 
fading diffuse multipath signals was that the amplitude and carrier 
phase of the coherently fading diffuse multipath signal are time varying 
random processes. The randomly varying amplitude and phase results 
from the fact that the scattering area from which multipath radiation 
is received by the aircraft is much larger for the diffuse case than 
for the specular case. Although all the multipath components received 
from this area have essentially the same modulation delay (this is the 
definition of coherent fading), their carrier phases may differ by many 
radians. The fading of the resultant multipath signal occurs primarily 
as a result of the differential doppler across the scattering area. 

The randomly varying carrier phase is by far the more 
important difference so far as steerable null performance is concerned. 
This is so because fluctuations in the amplitude are readily compensated 
for by the AGC system. If the multipath amplitude gets too small for 
the steerable null control system to track properly, the multipath will 
be too small to cause any problems anyway. On the other hand, as 
the multipath amplitude gets larger, the easier it is for the steerable 
null control system to track and null the multipath. The behavior of 
the multipath carrier phase has a much more significant effect on the 
behavior of the steerable null control system, since if this phase 
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dwells too long at an unfavorable value the steerable null system 
may lose lock momentarily. This effect was noticed with the 
simulator when, specular multipath interference at 136 MHz was 
being simulated. With a randomly varying multipath carrier phase, 
as would be experienced for a coherently fading diffuse multipath 
component at 1. 6 GHz, the likelihood that the phase will dwell for 
too long at an unfavorable value is much reduced. This was indeed 
found to be the case. 

In the simulation experiment conducted under Task V of 
the present contract, the phase of the multipath interference was 
randomly varied by driving a General Radio Frequency Synthesizer 
(GR1164-A) by Gaussian random noise of 2 0-2 00 Hz bandwidth in a 
VCO mode. The output frequency of this synthesizer was set at 
1 MHz (nominal) and the synthesizer was used to drive the multipath 
multiplier chain. The phasing in the antenna simulator was adjusted 
so that the multipath would cancel and the direct path would add in 
phase at the output of the summer. The phasing of the multipath 
component in the antenna simulator was then changed so that the 
multipath components at the input to the summer would be approx- 
imately 30° away from their null phasing. This was done with the 
loop open and closed and the waveform of the signal at the output of 
the summer observed on the scope. Scope photographs of this sig- 
nal are shown in Fig. 26 . The upper waveforms in both photographs 
are for the case where the multipath is nulled by adjustment of the 
phasing in the antenna simulator, and the lower waveforms in both 
photographs are for the case where the multipath phase is shifted 
away from null in the antenna simulator as indicated above. The 
effect of the steerable null control system is the nulling of the 
multipath component at the output of the summer by a feedback 
compensation for the improper phasing at the inputs to the summer. 
This is evidenced by the greatly reduced ripple on the signal at the 
output of the summer when the loop is closed compared with that 
when the loop is open. For convenience and facility in interpretting 
the various waveforms in the system, a simple 100 kHz sinusoidal 
signal was utilized as the modulation. 
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Fig. 26(a) Input Signal to Summer (upper) and Output 
Signal of Summer (lower) for the Antenna 
Simulator Adjusted for Multipath Cancella- 
tion (Loop is Closed). 



Fig. 26(b) Same as 26(a), but for the Antenna Simulator 
Adjusted so that Phases of the Multipath Com- 
ponent at the Input to the Summer are 3 0° Away 
from Null. 
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Addendum 


Computation of Relative Multipath Carrier Doppler 
for an Aircraft - NAV/SAT Link 


The following addendum contains a revised derivation of 
the differential carrier doppler f ms which avoids a potentially 
invalid approximation in Sec. 2.2. 1. This is so, because the 
approximations used to compute T mgJ the relative specular multi - 
path time delay would, in fact, yield zero differential carrier 
doppler. 

Conditions : 

1) Aircraft is at altitude h. 

2) Multipath grazing angle is y. 

3) Component of aircraft velocity parallel to the earth's 
surface and in the plane of the plane of incidence of 
the NAV/SAT signal is v^ . 

4) Component of aircraft velocity along the local vertical 
is v h . 

Find : The relative contributions of v^ and v^ to the relative 

multipath doppler. 

Geometry; See Fig.l . (Note that the notation used for certain 
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The phase of the received direct-path signal is 

<J> , = CO t - CO “ ( 1 ) 

d c c c 

whereas, the phase of the received multipath signal is 

r (t) r (t) 

A - co t — co — co — — — (2) 

ms c c c c c 


where r, rj and ^ are the direct-path range and the range seg- 
ments of the multipath range. The relative multipath doppler is 
taken to be 

« 

■ ( r 2 does not enter into the 

” r i* problem as maybe verified 1 

by direct calculation) 


This is equal to the difference between the projections of the air- 
craft's velocity vector along r and along r^. (The satellite is 
presumed stationary with respect to the earth, although for com- 
pleteness, the effects of satellite motion should be computed. ) 


• • 

In terms of v h and v^, (r - r^ is 
( -cos(h) - coso^) + v^ (cos(-^ — 0 ) - cos (~ - 0^)) 


or 

(r - r, ) = - v n (cos (h) + cos a . ) + v 
1 h v 1 

where 0 is the angle between the local vertical and r, and is 
the angle between the local vertical and rj. We wish to express 
this in terms of the aircraft altitude h and the grazing angle y. 


4 > 


(sin 0 


sin Oi i ) 


(4) 
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s 



E 

R-5960 

The following relationships may be derived from Fig. 2. 


a 

sin a, - , , 

1 a + h 

cos y 

(5) 

a 

sin a n - . , 

2 a + d 

cos y 

(6) 

1 7T 

<b = - 0! 

M 2 1 

-y 

(7) 

1 7T 

^2 2 " a 2 

- y 

(8) 

<l> = ^ + <t> 2 

sin <j> 

1 

(9) 

r = (a + h) 


(10) 

1 

cos y 
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( 11 ) 


r 


2 


sin <j> 

(a + d) 

cos y 


r - j (a + h) 2 + (a + d) 2 - 2 (a + h) (a + d) cos <{> 


2 2 

r = I r, + r„ + 2 r, r n cos 2 y 
\ 1 2 12 


sin A 6 = sin 2y 

r 


sin 6 - sin 2 y 

r 


5 + A 6 - 2 y = 0 


(h^tt-6-o^ - ?r - 2 y - a^ + A 6 


Let us examine the factor 

(sin (h) - sin a ) 

(0)= 7r - 2 y - + A e 


use 


and 


a 


7 r , 

12 7 " ^1 


sin (S) = sin (a ^ + 2 y - A0 ) = cos (y - <(> - A0 ) 


sin = cos (y + <j>j) 


( 12 ) 

(13) 

(14) 

(15) 

(16) 
(17) 
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(18) 


so that 

(sin (h) — sin a) = cos (y - <J> - A 0) - cos (y + <fi ) 

If y is not too small, and h/a is small, then <j)^ (the angle 
at the center of the earth between the scattering point and the 
airplane) will be small. In the flat earth approximation, <j)^ = 0. 

If the satellite is very far away (synchronous orbit), and 
if h is small, A 0 (the angle between the directions of propagation 
of the direct-path and multipath signals from the satellite) will be 
small. In the parallel ray approximation, A 9 = 0. 

If ^ is very small, and A0 is very small (compared to y), 
(sin (h) - sin o^) will be small, and will be zero in the parallel ray, 
flat earth approximation. In previous work, the parallel ray, flat 
earth approximations were made in order to relate the relative 
multipath delay T ms , to h and y. The real spherical earth, non- 
parallel ray situation was used to express y in terms of h and the 
angle <f> (see Figs. 1 and 2) (and the earth radius, a, and the NAY/ 
SAT altitude, d). The time delay r mg was then differentiated with 
respect to time, and expressed in terms of <}) (v^ = (a + h) cj>) and 
the other pertinent parameters. Since the expression for y in terms 
of h and <|> was only approximate, the validity of the value of T ms 
must be questioned, particularly, since T ms was expressed in 
terms of y and h using the flat earth, parallel ray approximation. 
The correct (exact) derivation (sin ffij - sina^) reveals that the 
contribution to the relative multipath .doppler from the horizontal 
aircraft velocity component vl is non-zero only because the values 
°f 4> i and A 9 deviate slightly from zero. Consequently we shall 
examine (sin (0) - sin a) more carefully. 
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cos (y - ^ - A0) - cos (y + (jjj) = cos (y - + -~ )) 


A Q 


i A0 , /i , A0 u 
- cos (y — + 


. A0 . . ,i AO . 

2 sin (y — ) sin ((J > 1 + ) 


(19) 


, A 6. . , Ad . 

sin y cos ( — — ) - cos y sin (— — ) 


sin 


(j>^ cos ( 


AO 


i . ,A9 \ 

+ cos 9 sin ( — — ) 
X z 


( 20 ) 


= 2 


2 AO 1 

sin y sin ^ cos (-^— ) + sin y cos ^ sin (AO ) 


2 Aft 

cos y sin cf> sin (A 6 ) - 2 cos y cos <J> sin" ( — — ) 
1 12 


( 21 ) 


2 sin y sin <|> + sin y cos <|> sin (A0 


- cos y sin <j» sin (A0) - 2 cos y cos <f> sin 2 

1 1 Cl 


2 A 0 

- 2 sin y sin <f> sin ( — — ) ( 22 ) 

For y not too small, the above is written in order of 
decreasing importance (i. e. , sin A0 < sin i^). 

To show that sin A0 < sin (j)^, express sin A 0 and sin <f>^ as 


r i 2r i 

sin A 0 = sin 2 y cos y sin y (23) 

r r 
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(24) 


sin +i = sfe) C03 r 


or 


which is 


sin A 0 = 


(a+h) 


cos y 


2(a+h) 


sin y 


sin A 0 = sin 


2(a*h) 


sm y 


since for y not too small. 


2 (a+h) 


sin y < 1 


Consequently, the first two terms in the expansion of 
(sin (h) - sina^) are most important. 

o/K \ 

With cos i =1-2 sin I 1 

1 \ 2 / 

(sin (h) - sin oij) = 2 sin y (sin <}> + ~ sin A0) 

- cos y sin sin (A0) - 2 cos y cos <j) ^ sin 2 

2 / A0 \ 2 / ^1 \ 

- 2 sin y sin <f>^ sin“ f — — j - 2 sin y sin (A0 ) sin ( — — j 

sa 2 sin y (sin <{> + — sin A0 ) 

1 Ci 


(25) 

(26) 

(27) 

(28) 

(29) 

( 30 ) 
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Now sin <j)^ = sin - yj - cos (a^ + y) = 
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If (2h/a) « sin 2 y, 

I For example: If h = 5 miles 

a « 4 000 miles 

2h = _1_ 
a " 400 

sin y » — , or y » . 05 radians j 
2 0 I 

then the square root may be expanded in a Taylor series to obtain 


(sin ^ - sin a ) =» 2 sin y 


, i a+h 
1 + [ I sin y 


a \ . 

— ' sin y cos y 




(t) 


2 

sin y 


> + 4 - 

2 \ a 


1 - 


. 2 
sin yj 


(44 cos 4 4 (t) 


.2 

COt y 


„ , / a+h , 

1 + [ J sin y 


(36) 


which for h/a «1 may be further approximated by 


1 

! h \ 2 

a 

1 " T 

— cot y 

1 + — siny 

2 

^ a / 

r 


(37) 


Additional algebraic manipulation or automatic computation 
may be used to reduce this result still further, Computation of 
(cos ^ + cos a^) follows similar lines, except that this factor will 
not be small. Consequently, the relative multipath doppler is much 
more sensitive to vertical aircraft motion than to horizontal air- 
craft motion. Because of this unanticipated result, the effects of 
satellite motion should also be considered, since it may contribute 
an unexpectedly large proportion to the total relative multipath 
doppler. # Both changes in cjj and in d should be considered. Note 
also that ^ will not be zero in this case. 
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Appendix A 


CALCULATION OF n FOR CASE (b) 

k 


o o 


+CJ V' 


^ dt„ E, ( c £ . (t. - t + 1 


2 <{, „ \ k\"l 'o co 


cC lt„ - t + 


c St 


n 


ci 


Kn ( V' 

n 


k\ 2 o to 


c S. 




( A -1 ) 


where E, and E denote statistical averaging over the clock phase 
9 „ n 

cjC v 

n 

jitter and the video Gaussian noise. 

Now, since n v is assumed to be "white" (i. e. , much larger band- 

•>: 


width than the bandwidth for h-, p) and of double-sided noise power density 
4-/2, E n /(n v ®h Lp y (n v ®h Lp ) \ is 

\ G ‘2, 




(A - 2) 


where 


r oo 

(T) = I h LP (t)h LP (t + T)dt 


LP -« 


(A -3 ) 


Introducing the function P^^,(t) defined by 


p iT <‘> 


1 0 < t < iT 

0 elsewhere 


(A-4) 


cr^ may be written as 


CT S = ° 2 f°° dT K (T >‘V (T) = 17 6 y 


where 


‘LP 


2 

a --n 


2 2 

k 


(A -5) 


(A-6) 
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and 


V T> = C d ‘ ^ ' *° 

n 


+ C 1 <th 

n 

cj „ 
c SL 




+ci " + T) ' 


n 


* C-Mt - t + T + ■ 

k\ o w 


cl! 


lp i T (t + T) ' 


(A-7) 


To proceed further, a specific lowpass filter function is chosen. For the 
sake of simplicity, a single -pole lowpass filter impulse response is 
chosen. The 3 dB filter bandwidth is taken to be equal to five times the 
clock frequency. That is 

h Lp (t) = ae~ ot u_ : (t) (A -8) 

where 


a - 5oj , 
c £ 


In this case, 0: (T) is 

h LP 


‘LP 


, a -or 

(T) = 2 e 


(A-9) 


Since is a very narrow pulse compared with the pulse 

width of ^i c ( r ). which is 2iT, it is only necessary to determine i/> k (T) for 
small r. To determine just how small a range of |t| need be considered, 
notice that for |t| = T c ^ = 2-nj^ c ^, 


* 


. (T ) 
h LP Ci 


1 1> 


h 


( 0 ) 


LP 


- 1 On 
e 


2. 27 x 



(A-10) 


which is entirely negligible. Consequently ^{T) need only be considered 
in the range |t| < T In this range 0 k (T) may be shown, for the case 
of no clock jitter, to be given by 


<£ k (T) - 1 


T b 

|_l 

M 

! 

2jtj) 

(T 

u/ 2 ) 


for t 



(A -11) 
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where 


T 


k 



T - 2 N T . (N-12) 

ci 


The factor T n / XT is 
k 



XT N 4, 096 I 

u X. 


<A-12) 


which is much less than 2 for k small compared to N. 

For <J> CJ g ^ 0, = 1, but will be slightly less than 

b k (T) | for <t> c ^ n (tj = 0, since the effect of the clock jitter is to further 
decorrelate the functions 



and cX It 


t + r + 

o 


n 

CJ . 

c X 


The amount of decorrelation increases as the ratio of the bandwidth of 
the clock loop to the frequency of the clock increases. Since this 
ratio is approximately 6. 4 to 320, 000, the amount of decorrelation, as 
a result of a small non-zero clock phase jitter, in the range of |t| < T ^ 
is entirely negligible. 

2 

Consequently a is approximately 

JC 


2 r 

CT J 


cX 


a 


- T c. 2 


“M/j . ( 2 .JTj It I \ dT 

l \ 2lT/(T k /2) I 


2 H " X /l 

- ) dx e 1 


oT a XT 
k 


(A-13) 


2 A 1 

a 1 + 


5rX2 5 tt2^ 


N 


2 2 
= a \ 


13 2 

Since 5rX2 » 1, the function y is approximately 

1C 


2 _ , 1/1 
y 1 / 

k 5r\ 0 (k-l) 


(A-14) 
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A table of values of y for k = 1 , 2, ••• 3 1 2 is given below 
in T able A. 1 . 


Table A. 1 


k 

2 

y k 

1 

0. 936338 

2 

0. 968169 

3 

0. 984085 

4 

0. 992042 

5 

0. 996021 

6 

0. 998011 

7 

0. 999005 

8 

0. 999503 

9 

0. 999751 

10 

0. 999876 

11 

0. 999938 

12 

0. 999969 
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Appendix B 


COMPUTATION OF <s k ) FOR CASE (b) 


For case (b) s is given by 
K 

s k = / iT(C)1 ' ,2sir ^ cos i c <t) ( b i3© h L p) C ' f k(’ 


♦ci T 


t -t + 


n 


O 00 


dt 


ci 


- op'j /” d« h Lp (S) ^ / * 1 b (* - 5)ci It - t o + -f—jcos A c (t)dt 

N ci { 


-00 


1 f CT 

o 


Ue <tK 


where 


F - (C) 1; 2 siniS (iT)p 13 


and use has been made of 


( b 13 ® Up) = d? h LP (?)b 13 U ~ ° 
\ —co 


The expectation value of s is 

k 


<\> =U 


1 r 00 


13 -oo 


p uu 1 rJfT 

f _ h LP^ It f b i3^ " ^ 

o 


<!> 


ci 


X E , ci, It - t + 
9 c£ ' k ' ° °° 


n' 


)E^ (cos <j> c )dt 


ciV ‘ c 

The expectation value E , (cos <J> ) is 

9 C 

E, (cos + c ) = /cos(+ ct + + cn l) = COS + ct E. (cos 4> 

T r> r r> ' ' T r>n 


cn 


cn 


cos 4> ct e 


1 2 

'9 CT 

2 cn 


(B-l ) 

(B-2) 

(B-3) 

(B-4) 


(B-5) 


(B-6) 
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where a 

cn 



The cross-correlation function between b^(t) and ci^(t) is 
defined as 

♦c/b < r > = Y / T b 13 Ct)C - S k (t + (B ‘ 7) 

k 13 o 


It is a property of BINOR codes that 


b (T) 

C k 13 


*ci b ( - T) 
k 13 


(B-8) 


Interchanging the order of statistical averaging over 
integration over t, (s ) becomes 



and 


1 2 

/ \ , , 2 cn i „ 

< s k> = M cos 't’ct E 4, 


c i 


n 


4> 


1 r 00 


ci 




13 -oo 


k 1 3 x ci 


(B-9) 


where use has been made of Eqs. (B-6), (B-7), and (B-8). Now 

/ f 


ci 


f * h LP («>* c p „ (‘o - T7 - y -(h.P®* 

k 13' ci 


-00 


T p 'O' r n u 

\ 13) 


(B-10) 


t - 


^ci 


n 


o oj 


ci 


is the convolution of ip. . . and h T _ evaluated atft - <b „ lo 1 Defining 

ci, b, n LP \o ci / ci/ 

/ k„\ k13 


x„ t 


n 


k\ o u> 


by 


ci ' 


ci 


X, It 


n = 1 


k V° w c J P 13 l LP’-'Xb 


(B-ll) 


k 1 3 J 


t - 


^ci 


n 


o w ci 


(s ) may finally be written as 


1 2 
-— CT 
2 cn 


ci 
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(s k )^e cos * ct ^ \Xk\fo 


n 


Cl) 


(B-12) 


ci ' 


n 



Appendix C 

COMPUTATION OF jX R | 


n’ 


7^cx 
7 ~^-t 
l °> 


A 


cl 


r cl ’ 
n 


The expectation value of X, ( 0- _ t ) is given by 

°f 

\ - € k U - ln2) - 


v\ c.0 / J n ' n/ 


— in2 
5 


+ ^„ P (*=0 


2f {1 - in2) 
k 


5(j> 


cl 


“V 


n 




cl 


n J 


with 

_< t’ 

P (*'0 = e 

n 

The first integral in Eq. .(C-l) is simply 

h = 1 - V 1 - ln2) ' 5 t k ff ci 

n 

The second integral may be evaluated by using the following trick, 
the integral 



Ki 5 "tci 

G(o) = / ° e n P 


M d * 


cl 
n' n 


In terms of this integral, the second integral in Eq. (C-l) is 

dG 


I 2 ‘ 


2 e (1 - in2)G(0) - e. G(l) + 2e 
k k k da 


o = 0 


(C-l) 

(C-2) 

(C-3) 

Consider 

(C-4) 

(C-5) 
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Consequently, we need only evaluate G{a) and specialize the variable 
a to obtain the various terms in 1^. 


o<«>-/ ° 


*cf 54 cf 


-00 


-t 2 2a 2 , 

, / cH 

n 1 n n , i 

Tsrrr e ^ 

c SL n 

n 


(C-6) 


Completing the square in the exponent. 








(C-7) 


G(n) becomes 



G(0) 



+ erf I 


Sin 2 


V~2 5ct 


CJf 


n 


(C-10) 
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2 


G ( 1 ) - e ' 


2 1 - 0rf 


_n in2 

2 y/2 5 a 


(C-ll) 


1 2 


2\ v 2 w? 


V2 5a 


( 012 ) 


Introducing the above into (C-5) and rearranging terms, E . ( xj — 

becomes mi \ \ g o 




ci 

C n 2 

+ “vT7v e 


V 2 5a 


/5a , \ r ,5a „ 

( Ci n) 1 / Ci n in2 

+ e VvT7 2 1 ^"TT" " 7T^~ 


(C-13) 


C k 5r (1 " V 


171 




preceding page blank. N©T...Piy^ti3- 

Appendix D 


i 





This is expanded in a power series about cf> ^ = 0. To do this the follow- 
ing Taylor series expansions are used: s 
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( D— 4 ) 



In Eq. (D-5) use has been made of the fact that 



<D-5) 


(D-6) 


Introducing the above series in Eq. (D-l), and grouping terms of like 

that E 


s 



powers of <j> c ,g g , and realizing 
to second order in <j)^ 


V c 
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(D-7) 


Examination of Eq. (D-7) shows that in the expression multiplying 
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and sixth terms cancel, the second part of the second term and the 
second part of the third term cancel, and the first part of the third 
term and the second part of the sixth term cancel. This leaves only 
the fourth term. Consequently Eq. (D-7) becomes 




ci 


o CO 


ci 


' Vmax ‘ E k r <V 


- € 


TI / . | 2 1 2 .2. ,1 2 

H (to ) -a sin ho r ] -(o) 
ci ms 2. ms <- r-o y 


c S. ms J 2 


(D-8) 


iir) l[ I + erf (v? 

x e . vv 


in2 

42 


5a 


ci 


where use has been made of Eq. (50), which defines I (a .), and of 

c ii 

the relation 


E 





(to ) 

ms 


2 1 2 

— cl 

2 ms 



(D-9) 
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Appendix E 


COMPUTATION OF 


<\> 

\ k/ms 


In estimating ^ s kXns' 1 ^ ie sma ^ effects of <j> c and <f> cJ g may be 
neglected. That is 


(sf) = a 2 C sin 2 # E 
\ k/ms ms < 


ms 


/^ T dt i I** dt 9 coS d <V C0S d (t ) 

J 1 J 2 ms 1 ms 2 


iT 

o 


X ( b 13® h Lp) _ t ( b 13® h Lp) T 'VVV'VVV 

1 ms 2 ms 


To proceed, we write 

cosfd (t)l « cosho t+d (0)1 
L ms J L ms ms J 


(E-l) 


(E-2) 


for 0 < t < i T. That is, it is assumed that co mS , the relative carrier 
Doppler frequency of the specular multipath component, is essentially 
constant over a time interval on the order of iT. Then, averaging over 
d ms corresponds to averaging over d ms (0) which is assumed to be 
uniformly distributed over 2n . When this is'done it is found that 

E (cosfd (t )] cos[d (t )j\ = \ cosfco (t - t )] 
d l L ms 1 J L ms 2 J 2 L ms 2 1 1 

ms v 


1 . 2 “ ms (t 2 ' V 

° 2 2 


Then is approximately 

\ k/ms 


(E-3) 


( S k>ms= a L C ^|(/ iT dt(b 13 ®h 


cje (t-t ) 

Lp / t T k ° 

1 ms 


> 2 CsinV^sin 2 
ms -22T 


co r 

ms 




t -T 

1 ms 


cje (t -t )p._(t ) 
k 1 o JLT 1 


X ( b 13® h LPj t _ T +7 . C ^k (t l " *0 + T '* P ^T (t l +T) 
'1 ms 


(E-4) 


177 



where 


P^(t) is defined by 


P 


HT 


(t) 


1 

0 


0 < t < IT 
elsewhere 


(E-5) 


An upper bound on \ s ^/ ms is obtained by neglecting the second 
portion of Eq. (E-4). If u> mg is small compared to 2irj2lT, this term 
will be much smaller than the first term. The first term is simply 



a 2 (C 
ms 


sin 2 |S(i T) 2 


fa) 


xf(t 

k o 


- T ) 

ms 


1 2 

2 a ms 


2 2 
M X, <t 
k o 


- T ) 
ms 


(E-6) 


which is less than or equal to 
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COMPUTATION OF (njy d 


The properties of x (t) and x (t) discussed in the beginning of 


- 7 9\ H 

Sec. 6 may be used to write approximately as 


("k)d =Csin ^ E d U dt l/ *2 Ub ' V 

ms I o o 

cosfd (t )l cosfd (t 9 )l + sinfd (t )lsin[d (t )] ] 
L ms I J L ms 2 1 L ms 1 J L ms 2 J i 

P )‘C, 


x fb J3 <g) h L 


ci (t - t ) b 
k 1 o 13 


LP ‘VW 

2 ms 


(F-I) 


Now 


cosfd (t )1 cosfd (t )l+sinfd (t )] sinfd. (t )] = cosfd (t )-d (t )] 
L ms 1 J L ms 2 J L ms 1 J L ms 2 1 1 ms 2 ms 1 J 

(F-2) 


so that if 


d (t) « u t + d (0) 
ms ms ms 


over the range 0 < t < 



becomes 


( n k)d = C sta2(3 /" ®l/" dt 2 t m < l 2 - V COS [“ ms (t 2 - tp] 


*1^13®^, T C VVV( b 13® h Lp) t T ‘VVV 

1 ms 2 ms (F . 3) 


To proceed further, the effects of the lowpass video filter and 
the phase shifter time delay are neglected as being of only minor 
importance. In this case (n^ ^becomes 



C sin f dt f dt <() (t - t )cosfu> (t, 
J 1 J 2 T m 2 1 L ms < 

o o 


x b (t - r )cl (t )b „(t„ - t w (tj 
13 1 ms k 1 13 2 ms k 2 



(F-4) 
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The autocorrelation function J> (t„ - t ) is now expressed in terms of 
., „ . , . m 2 1 

its Fourier transform 


<t> (t - t ) = J ~ ( 

m 2 1 J 2n m 


j*><t ,-t. ) 


u)e 


(F-5) 


-oo 


where <£ (w) is the power spectrum of x c (t) and x^(t). It is useful to work 

with since, although its precise form may not be known, an 

estimate of its bandwidth can be made. Writing cos [ CJ ms ^ t 2 - as 


j“> <t 0 -t ) -jco <t -t ) 

ms 2 I , ms 2 1 
e + e 


COS l“ms tt 2 ' Tb = 5 
and making use of Eq. (F-5), <\> d may be written as 

/ 2\ r 00 dw rfi T rS . T ^^2^1 ) 

\ n k>d" Csin ^/ 2^*m (0j) f dt l J dt 2 6 


(F-6) 


-00 


g m 1 + e ms 1 lb (t -T )ci (t )b (t -T )ci (t ) 

13 1 ms k 1 13 2 ms k 2 

(F-7) 

The usefulness of this procedure may now be more fully appreciated in 
that examination of Eq. (F-7) reveals that the variables tj_ and t2 may now 
be separated as follows. 


/ 2 \ 1 . 2- 2 2 r 00 dw . . 

\ n k)d~ 2 C sm ^ WT) p l 3 / co S* m ( “ ) 


1 

1 

, P 13 

' iT 

1 

1 

P 13 

iT 

1 

1 

P 13 

iT 

1 

1 

P 13 

iT ■ 


r C T 

/ b (t - T )ci (t )e 
J 13 1 ms k 1 


-j(oo+w )t ■ 
ms 1 


« / b 13 (, 2 ' T m .>°VV 


j(co+w )t- 
ms 2 


f ' T )ci, (t. )e 

iT J 13 1 ms k 1 


-j(u-w )t 
ms 1 


b (t - r )ci (t )e 
J 13 2 ms k 2 


m s <-i 


(F-8) 
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Since b 13 (t - T ms )ci^(t) is periodic with the period of the BINOR 
code, T, it may be expanded in a Fourier series as 


b,„(t-T M, (t) = V o (t )< 


Jp “l3 l 


13 ms k 


: -00 


p ms 


(F-9) 


a 


(T ) | iT b (t - T )cJL (t)e 

ms ■ IT J 13 ms k 
o 


■^is 1 dt 


(F-10) 


Since b-^t - T m s) c ^k^ is rea1 ’ “.p( T ms ) = °p^ T ms^ In particular. 


a o< T ms } is 


a (t ) - — f b (t - t )ci (t dt 
o ms SLT J 13 ms k 
o 


= +ci„b, *- T ms ) * W < T m s> = hsVW 
k 13 k 13 


(F-U) 


o 


The superscript "o"in X r s ) indicates that there is no video filtering 
being considered. The effects of the other coefficients c <p( T ms ) for p f 0 
are shown in Appendix G to be of negligible consequence. 


Introducing Eq. (F-9) into Eq. (F-8), and making use of the 
relationship 


]W t 00 

1 rlT dte ° b (t - t )ci (t) = Yj a (T )j- f 
13 ms k ^ p ms HT J 


— r J 

j?T J 


1 r & T 


, = -00 


dt e 


j(d +w )t 
O p 


o 


j8T 


00 j(w +w ) 

= y a (r > e 0 p 2 

LJ _ p ms 

p =-00 ^ 


sine 



+ w 

p 



where w = pu in , 
p 13 


and 


sme x 


_ sin x 


(F-12) 


(F -13) 
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^ become s 

/ 2 \ _1 2 r °o dto ^ .. 

W)d"i M / 5tV“> 


CO 


Si x 

j(w -to-co )— 


V v ' p ms 2 . r/ 

> a (T )e y sine { 

l—l n m s l 


i Pi o p ms 

L 13 p=-*> ^ 


JT-, 

CO “CO “CO ) 

p ms 2 J 


00 


i T 

, i(co -Ho+to )— — 

1 v /_ x J r ms y 2 r/ , 

— ) o (r )e sine (co + to + to )-r- 

r m s *- r ms 2 J 


Pin t ~‘ ^ r ms 
*13 r = -°° 


+ 


00 


iT 

. i(oo -<o-ho )— — 

-J- V a (T )e p ms 2 

P 13 p = -» P mS 


iT 

sincRto - go + go )— — 1 

L p ms 2 J 


00 


iT 

i(co -ho -to ) 

y 1 v r ms 2 . r - v iT, 

* a (T )e sinc|(oj + w - u )tH 

p, „ _ r ms 1 r ms 2 J 

13 r = -°° 


Examination of the products 

$.1 i. T 

sine ( io to to )-— lsincr(<o -f-to+to )—— I 

p ms 2 J . 1 r ms 2 J 

i T . & T 

- sinc[(to+co -to sinc[(co +to +co )— — 1 

t mS p 2 J L rn ft r* 2 -* 


ms r 2 


and 


:[ (to - to + co 

)“] sinc[ 

(to + to - to )-— 1 

L P 

ms 2 J L 

r ms 2 J 


(F-14) 


(F-15) 




x i T 

= sinc[(to - to - co sincf (co - to + to )-— 1 

L ms p 2 J L ms r 2 J 


(F-16) 
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indicates that they are essentially negligible except for r = -p. Retaining 
only the r = -p terms, Eq. (F-14) becomes 


K>, 


2 \ 1 2 /-oo dw . . 

U F ^ / T" * < w 

d 2 J 2-n m 
oo 


00 


)' ^ ot (T )/* (T ) sinc"[(w-w +co ) yy ] 

2 p ms -p ms " ° 


p^_oo p 


2 ■ XT 

p ms' 2 


13 


2 i T 

+ sine r (to - w - w )-r-l 
L p ms 2 J 


This may be written as 
2 


/ 2\ 1 2 1 2. poo dw . , . 

\n ), % -5- o' < T ) / “ $ <<*>) 

\k/d 2 2 o ms 2 tt m 

P- - 00 


13 


(F-17) 


2 j£T 2 j?T 

sine f(w + w )-yl + sinc“f(w - w )-yl 
L m 2 J L m 2 J 


CO 

2 1 . . . v r 00 dw , . 

+ ^ X - 5 - a (r )a (T )/ — $ (w) 

2 p ms -1 ms J 2ir m 
P=1 p 13 o 


jsinc [(» -P“ 13 + “ ms >-rl 


PT 2 

v — ] + sine [(w - pw 


viT, 

13 " “ms* 2 1 
(F-18) 


Assuming a Gaussian shaped spectrum for <i> (co) 


a 2 -(f/B ) 2 

$ ( W ) = --g- e m 

m Jir B 
v m 


(F-19) 


where ctT is the total power associated with the spectrum $ (w), 
m m 


a 2 = cf) (0) = f $ (w) 

m m J 2, ir m 


-00 


00 dw 

2tt T m’ 


(F-20) 
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it is assumed that B m is sufficiently large so that over the bandwidth of 
the functions 


■ 2 r, JT, , • 2 .iT, 

smc |(u + u )— + sine f(u - u L 

l m .q 2 J l m s z J 


<£ m (u>) may be replaced by $ m (0). The p = 0 term in (r ijT) ^ is then 


/ 2\ , 1 2 1 2,_ , 1 CT m 

\ n k/d P = ° = ^ _ ^ _Q ' (T 


2 2 o ms 2 t: Jtt B 

P 13 


< 2 >W* 


(F-21) 


where use has been made of the definite integral 


r 00 . 2 

J sm x 

-« 2 
x 


r oo 2 

dx = ir = J sine x dx 
.00 


2 / 2 \ 

Introducing the expression for o (T g ) into Eq. (F-21), yi ^ (p = 


(F-22) 


0) 


become s 


(n l) (p = 0) = p 2 o 2 r 1 (x °{T )\ 

\ k/ d m i it B iT\ k ms I 

m 


(F -23) 


This is by far the most important contribution to and is 

the only one which will be. considered in the computation of (P^)^,. 
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Appendix G 


ESTIMATE OF EFFECTS OF tt (T ) a (T ) FOR p i 0 

p ms -p ms 


til 

In Appendix F a (T ) for the k frequency-divided down 
wave was defined as r 


a. 


JIT 

<T ms ) = rr / b !3 (t - T ms )c V t,e 




dt 


Now b, „(t - T ) and c £ , (t) were written as 
13 ms k 


>, Q (t - T ) = )' 0 e 

13 ms t-j n 


j n “is' - jn “l3 T ms 


n = -co 


and 


A °° 1 

ci k (t) % y sin t ( 2 q+i) “k‘] 

q = 0 


4 ™ 1 1 ^ j(2q + 1 )“ k t -j( 2 q + l)o> k t 

~ h ^ (2q + 1) ~2j \ 
q = 0 


-e 


Introducing these expressions into Eq. (G-l), a (t ) may be shown 

pm 


. -ipu, „ T oo , 

a (r ) = i e 13 ms Y .7— i — 

- - Zj (2q + 1 2i 


P mS 7T 


q = 0 


0 e 

“k 1 
( 2 q+ 1 ) — — -p 

“l3 ■ 


j(2q+l)u, T 
J k ms 


- 0 


-j(2q+l) w T a 

k ms 


(2q+l ) +p 

“l3 


square - 


(G-l) 


(G-2) 


(G-3) 
to be 


(G-4) 
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Now /9 is 
n 


0 n = k / b i ,<‘> e 




1 r 

-j T J b l 3 (t)sin(nw l 3 t)dt = 


6 

. n 

-jy 


since b, „{-t) = - b „(t). From Eq. (G- 5 ) it is also seen that /3 = - j 9 , 

1 o lu “ n n 

so that a (T ) is 
p ms 


, . 4 jP a, 13 r mS /] L \ y 1 

“p ms ) ~ir e UJ Li ( 2 q+l) 

q =0 

L j( 2 q+ 1 V ms . „ 


-j(2q+l)w k T i 


( 2 q+l) p 

CO 

13 


( 2 q+l) +p 

W 13 


and o (r )a, (r ) = a. (r )o> (r ) is 
p ms -p ms p ms p ms 


1 « 00 1 1 
3 ^ms^-p T ms^ 2 ^ ( 2 q + 1 ) ( 2 s + 1 ) 


7 r q =0 s =0 


* 0 


i (2q+l >v ms , . 

e + p 


[L pq + D- -p] [( 

*[,* ' 


( 2 q+l) — +p 

co 

13 


-j( 2 q+l)u T 

k ms\ 


+j( 2 s+l)o) T ' 
k ms 


( 2 s+l) — -p 

OJ 

13 


r co 

( ( 2 s+l)— +p 

l 
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<M | CO 


which is 


00 CO 


a (r )a (r ) = )' Y 77; — 7 

p ms p ms 4 2 Li U (2q + 


1 


“ -n (2(1* 1) (2s + 1) 

7r q = 0 s = 0 


j(2q+l)u T 

5 e k ms + « 

^ Cl) * 

(2q+l )— ■— -pi 
13 J 


-j( 2 q + i) w ^ T mg i 


"k 

(2q+l ) +p 

to 

L 13 


x 


( 2 s+ 1 ) — 

_[ W 1 3 


-i(2s+l)co, T 
J k ms 


+ 6 


j( 2s+l)co T 


k ms 


co 


“P 


( 2 s+1 ) — — — ip 

L 


_k_ 

J 1 3 




rto, 1 

■ v 

r 3to. 

k 


k +n 

CO ^ 

13 J 


+ P 

L 13 


6. 


to. 




— +pf /c,os(2co r ) 
' ‘ k ms 


L w 13 V 



CO. 


to 


13 


-p I icos(4co T ) 
^ 1 ' k ms 



• 2 /o 

sin (3co. 


T 

ms 


)! 


j 


(G-8) 


+ . . • 

(G-9) 
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The dominant contribution to Q (t )or (t ) comes from the 

p ms p ms 

q = 0, s = 0 term. This is the only term which will be considered. 


From Eq. (G-5), 


“ k 

, ~ p 
L 13 


rp 

’I / b l3 (t)sin [“k t * p “l3 t ] dt 


4 r Tl2 [ 

= T J b i3 (t>sl Tk t * p “i3‘ 


dt 


(G-l 0) 


since b^f-t) = " b 13 (t). 


By utilizing the relationships 
b 13 (t) = 1 6 n 


and 


n=l 

00 


b 13 (t) “ 2 8 n < ““l3 )cos [ n “l3 t J 
n=l 


Eq. (G-10) may be written 


T/2 


“k 


W k . 

±p 


±p 

G) 

L 13 


. W 13 . 


/ b^ 3 {t)cos[w t t ± po> 1 ^tjdt 


13 


w T 
13 


(G-ll ) 


(G- 12) 


2 1 


7r r <»>, 

-a ±p 

“13 


T/2 i 

/ b 13 (t)cos[pw 13 t]dt 
0 


T/2 


n I 

J b 13 (t)sin[w w t]sin[pu 1 ^t]dt 


13 


(G-13) 
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This is convenient because b^t) consists only of impulses of areas 
+ or - 2 at the zero-crossings of b (t), and this greatly simplifies the 

X O 

evaluation of the integrals in Eq. (G-13). The sum 

b 13 (t)cos[w k t]cos[pw 13 t]dt 

(G-14) 



For those values of k less than about 4, the above is quite 

negligible because u /cj is equal to or greater than 256, and the inte- 

gral is much smaller than that. For values of k above about 8. sin (oj, T ) 

k ms 

is quite small for the values of T which would be encountered in the 

ms 

NAV/SAT application {see Eq. (G-9)), so that for these values of k the 

* 

dependence of a (T )a (r ) on T is essentially negligible. The 
p ms p ms ms J 

largest values of Eq. (G-14) occurs for values of k and p for which 

[co /co _ ± p] is a small integer. Since only small integer values of p 

are of interest, the values of k for which ± p] is small are 

approximately k = 9,10,11, and 12. Since a a (p^O) is negligible for 

P P 

all but the largest values of k, only the values of (P c ) for these k are 

affected by the terms involving a a (pyO). Since the effects of these 

P P 

few values of(P ). are small, the total effect on{P ) may be assumed 

C K ' C 

to be quite negligible in relation to the effects of the p = 0 component. 
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Appendix H 


TABLES OF FACTORS INVOLVED IN THE COMPUTATION OF 
THE PROBABILITY OF CORRECT BINOR ACQUISITION, (P ) 

C T 


Table H. 1 Factors in f Which are Independent of k 


Situation 


J (a 

o 


ms 




I 

0. 984 

0. 9582 

0. 8826 

2 

0. 984 

0. 9786 

1. 0000 

3 

1. 000 

0. 9745 

0. 9357 

4 

0. 984 

0. 9582 

0. 8826 


Table H. 2 Factors in f Which are Dependent on k 


Situation (X . ) -e, T(a ) 

k max k cj£ 



k=l 

2 

3 

4 

5 

6 

7-12 

1 

0. 9210 

0. 9636 

0. 9855 

0. 9953 

0. 9987 

0. 9999 

1. 0000 

2 

0. 9255 

0. 9660 

0. 9865 

0. 9955 

0. 9988 

0. 9999 

1. 0000 

3 

0. 9242 

0. 9655 

0. 9862 

0. 9955 

0. 9988 

0. 9999 

1.0000 

4 

0. 9210 

0. 9636 

0. 9855 

0, 9953 

0. 9987 

0. 9999 

1. 0000 
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Table H. 2 (Continued) 


2 


Situation 


1 2 
€. r a 
k 2 mg 

.(*?•' 

( ia cl\ 

( n ) i 

2 

1 - erfj 

' 5C7 c i 

in2 \) 

~ rt 5a ci /J 


k=l 

2 

3 

4 

5 

6 

7-12 

1 

0. 0137 

0. 0062 

0. 0025 

0. 0008 

0. 0002 

0. 0000 

0 

2 

0. 0148 

0. 0067 

0. 0027 

0. 0009 

0. 0002 

0. 0000 

0 

3 

0 

0 

0 

0 

0 

0 

0 

4 

0. 0137 

0. 0062 

0. 0025 

0. 0008 

0. 0002 

0. 0000 

0 


Situation 




k=l 

2 

3 

4 

5 

6 

7-12 

1 

0. 9073 

0. 9574 

0. 9830 

0. 9945 

0. 9985 

0. 9999 

1. 000 

2 

0. 9107 

0. 9593 

0. 9838 

0. 9946 

0. 9986 

0. 9999 

1. 000 

3 

0. 9242 

0. 9655 

0. 9862 

0. 9955 

0. 9988 

0. 9999 

1. 000 

4 

0. 9073 

0. 9574 

0. 9830 

0. 9945 

0. 9985 

0. 9999 

1. 000 
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Situation 


Table H. 3 


Factors in 



f (cos <j> is between 0. 989 and 1. 0 for Ir/gl v 1, 

K. Cl 

so it is assumed approximately equal to 1. ) 



k=l 

2 

3 

4 

5 

6 

7-12 

1 

0. 7550 

0. 79 67 

0. 8180 

0. 8276 

0. 8309 

0. 8321 

0. 8322 

2 

0. 8770 

0. 9233 

0. 9474 

0. 9578 

0. 9616 

0. 9629 

0. 9630 

3 

0. 8428 

0. 8804 

0. 8993 

0. 9078 

0. 9108 

0. 9118 

0. 9119 

4 

0. 7550 

0. 7967 

0. 8180 

0. 8276 

0. 8309 

0. 8321 

0. 8322 


Situation 


( 




r— 4 

II 

1A 

2 

3 

4 

5 

6 

7-12 

1 

2. 144 

4 

2. 263 

2. 323 

2. 350 

2. 360 

2. 363 

2. 363 

O 

2.491 

2. 624 

2. 691 

2. 720 

2. 731 

2. 735 

2. 735 

3 

2. 394 

2. 500 

2. 526 

2. 578 

2. 587 

2. 590 

2. 590 

4 

2. 397 

2. 530 

2. 597 

2. 627 

2. 639 

2. 642 

2. 642 
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Table H.3 (Continued) 



Situation 


1 

2 

3 

4 


Situation 


k = l 

0. 998786 
0. 999786 
0. 999645 
0. 999650 


2 

0. 999314 
0. 999897 
0. 999796 
0. 999827 


3 

0. 999490 
0. 999929 
0. 999823 
0. 999880 


4 

0. 999555 
0. 999940 
0. 999867 
0. 999898 


5 

0. 999577 
0. 999944 
0. 999873 
0. 999905 


6 

0. 999584 
0. 999945 
0. 999875 
0. 999907 


7-12 

0. 999584 
0. 999945 
0. 999875 
0. 999907 




k=l 

2 

3 

4 

5 

6 

7-12 

1 

0. 002395 

0. 001497 

0. 001167 

0. 001040 

0. 000997 

0. 000984 

0. 000984 

2 

0. 000716 

0. 000382 

0. 000274 

0. 000237 

0. 000224 

0. 000220 

0. 000220 

3 

0 

0 

0 

0 

0 

0 

0 

4 

0. 001061 

0. 000582 

0. 00424 

0. 000366 

0. 000346 

0. 000340 

0. 000340 
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Table H. 4 


Conditional Probabilities 



Situation 



k=l 2 


3 4 


5 6 


1 

2 

3 

4 


0. 996391 
0. 999070 
0. 999645 
0. 998589 


0. 997817 
0. 999515 
0. 999796 
0. 999245 


0. 998323 
0. 999655 
0. 999823 
0. 999456 


0. 998515 
0. 999703 
0. 999867 
0. 999532 


0. 998580 
0. 999720 
0. 999873 
0. 999559 


0. 998600 
0. 999725 
0. 999875 
0. 999567 


Table H. 5 Total Average Probability of 
Correct BUsOR Acquisition 


Situation 



1 0. 980 or 98. 0% 

2 0. 996 or 99. 6% 


7-12 

0. 998600 
0. 999725 
0. 999875 
0. 999567 
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Appendix I 

EXPANSION of cos (^bft) and sin (/? 9 b(t)) 
The expansion of 


co s 


(p 2 b ( t)) = 


V 

^ 2 

j = COS 

~4 b l{ t ) +b 2 ( t ) +b 3 ( t )+ b 4 (t) 


(I. 1) 


wh ere 


b (t), b^ft), b^(t) and b^(t) are unity amplitude square waves i. e. , 


b. (t) - ±1 and R. is the maximum deviation. 

i 2 


C0S l ~4 ( b l (t) +b 2 ( t ) +b 3 ( t ) +b 4 (t)) 


f ^2 

. , \] 



cos[ — 


cosj T 

( b 3 (t) + b 4<t))j 


- sin [T ( b i (t ) + b 2 (t f] [ sin -J ( b 3 (t)+b 4 (t) )J 


= { cos 


cos 


b 2 


r ^2 i i 

■^2 | 

It b 2 (t) 

- sin 

[ T b i (t, | sin 

T V‘>l 

I b 2 l 


13 2 l 

r^2 1 

It b 4 (t) ] 

. 

- sin 

T b 3 (t) ) sin 

It b 4 (t »J 


r 2 

sin|— b^t) 

I \ ^2 

x | sin l "T b 3 (t ^ 


COS 


. p 2 „ , 

f P 2 


(t V‘> 

+ cos — b x (t) 

sin 

[ b 2 . 1 

f ^2 

1 

It b 4 (t) l 

+ cos[ — b 3 (t) 

sin 


f b 4< 


(I. 2) 
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collecting terms and recognizing that b. (t ) = ±1, yields 


cos 


4 b 2 


b l (t)b 2 (t)+b 1 (t)b 3 (t)+b 2 (t)b 4 (t)+b 2 (t)b 3 (t)+b 2 (t)b (t) 


+ b 3 (t)b 4 (t) 


2 P 2 2 ^2 4 ^2 

cos — sin — + b 1 (t)b 2 (t)b 3 (t)b 4 (t) sin — 


4 ^2 

= cos — — ~ 1 

4 


(I. 3) 


for P 2/4 « 1 
The expansion of 


/ P 2 / 

sin(p2 b ( t )/ = sin T ( b 1 ( t ) +b 2 (t)+ b 3 (t)+ b 4 (t) j 


(1.4) 


is performed similarly 

P 


sin 


_2 

4 


(b i (t)+b z (t)+ b 3 (t)+b 4 (t)) = 


rp{ / 

r/ Sn V -I 

(I. 5) 

sin T n ^ + k 2 (t)J cos 

[(T <b 3 (t > +b 4 (t >)] 


+ co ‘[(t b l ,t)+b 2 (t) )j 

| Sin |!V b 3 <t)+b 4 (t) )j 



sin 


^2 

TV** 


COS 


13 2 r 13 2 i r P ? 

— b 2 ft) + cos — b 1 (t)l sin — b 2 (t) 


x {cos 


+ { cos 


T V‘> 


T V‘> 


sin 


( ^ 2 

T b 3 (t) 



' b 2 1 


’ ^2 

cos 

T b 4 (t) 

sin 

T b 3 (t) _ 


V 


P-2 

cos 

— v*’ 

- sin 

T b i (t) 


^2 


P 2 

cos 

T b 4 (t) 

+ cos 

T b 3 (t) 


ft 


sin \ — b 4 (t) 


sin 


sin 


^2 

T b 2 (t) 


T b 4 {t) 
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collecting terms and recognizing that b.(t) = ±1, yields 

/ \ ^2 3 ^2 

| b L ( t ) + b z (t)+b 3 (t)+b 4 (t) J sin— cos — 

-('b i (t)b 3 (t)b 4 (t)+b z (t)b 3 (t)b 4 (t)+b i (t)b 2 (t)b 3 (t) 

\ 3 ^2 ^2 

+ b 1 ( t ) b 2 (t) b 4 ( t ) ) sin — cos — 


E b. (t) 


sin 


for '2/4 « 1 


(I. 6) 
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Appendix J 

COMPARISON OF A 3-TONE BINOR CODE 
WITH A 3 -TONE LINEAR SUM 

For the 3 tone case, equal subcarrier amplitudes are achieved when 


cos (3 ~ cos ~ (J - 1 ) 

but from Equation 45 modified for 3 tones 
P 2 

tan — = . 221 (J -2) 


or 


^2 = 12.5° 

3 


(J -3) 


as before substituting {J -3) into (J -1) yields 
cos p = .93 


and 


p b = 21.5' 


J-4) 


Assuming optimum demodulation of the BINOR code, the BINOR 
system performance will be superior if 


p tan (3 ^ jr- 

3 b rr \j 2 


P 2 2 P 2 
sm — cos — 


(J-5) 


where p^ = ^ is the BINOR coefficient. Substituting from (J -1) and noting that 
P 3 tan p = (. 5 ) ( . 3939) = • 197 


and 


(J-6) 


P- 


4 r Z 

WS tan T 


9) (. 221) = . 1989 


it is concluded that the BINOR is not superior in performance to the sum of 
square waves for this case. 
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PRECEDING 


page BLANK not HEME Appendix K 

COMPUTER SIMULATION PROGRAM 
PHASE LOCK LOOP SIMULATION! G-106 


TYPE 20 

20 FORM AT ( 1 H/# 1 IMPUT THE FOLLOWING PARAMETERS V 

1 ' M-ZERO# OMEGA, SIMT# BDV.T DTH# DIV M# TGC# TGD# A# VCO '# / > 
ACCEPT 1 2# STDI V#OMO FF# FO# 8W# EM# TGC# TGD# A# T1 

12 FORMAT ( 9 F ) 

THETA= 0 

I RAM I T= 3088 3623 

CALL GAUSS < I RAM I T# STDI V# 0.0# RAM DOM ) 

WO F T = 0 
F00 = 0 
A 3 = 0 
R 3= 1.0 
T 1 = T 1 * FO 
PHI 2=0 
PH0 = 3 
DPH IW = 3 
I NJDEX = 0 
T2 = T 1 
E M P = 0 
EMM = 3 
ATCO S= 3 
A T S I M = 0 
TY?E 23 

23 FORMAT C 1 H #' STEP '# 2X# ’ VCO FREQ G AX# 'VCO 01JTPUT’#2X# 
1 'DIV 0 UTPUT ' . 2X# VM/PLUS *# 6X# 'M/MINJUS') 

STDI V=SORT C STDI V/FO > 

BW= 1 . 0- EXP ( - RW* FO ) 

I THETA=C-T2+DPHIM+THETA> 

THETA= AMO D( THETA# 6. 2531 8 > 

IFCTH ETA-3. 14] 59) 1 0# 5# 5 

5 THETA=THETA- 6. 2831 3 

GO TO 1 5 

13 IFCTHETA+3. 1 41 59) 1 1 # 1 5# 1 5 

II THETA=THETA+6. 28318 
15 COMTI.NJUE 

TSIM=SIM (THETA) 

TCO S= CO S ( THETA ) 

A TCO S = ( T CO S - A T CO S ) * A+ AT CO S 
A T S I M = ( TS I M- A TS I N )*A+ ATS IM 
CALL GAUSSf I RAM IT# STDI V# 0. 0# RANDOM ) 

ACOMP= T S I M + RAM DOM 

CALL GAUSS C I RAM I T# ST D I V# 3 . 0# RAM DOM) 

BCOM P = TCO S + RAM DOM 
WO FT=OMQ FF+FO + WO FT 
CO WT=COS( WO FT) 

SI WT=SIM( WO FT) 

A 2 COM ° = A COM P* CO WT+BCOMP*SI WT 
B2 COM P = BCOM P* CO WT - ACOMP* S I WT 
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A 3= ( A2C0MP-A3)*BW+A3 
03= ( B2C0MP-83) +BW+B3 
A 4= A 3* CO WT-B3* SI WT 
9 4=B3*C0 WT+43*SI WT 
PHI =ATAM2( A 4* B4> 

I FCABSCPHI-PHI2)-3* 1 41 5926) 103* 100* 101 

10! I F ( PH I - PH I 2 > 1 02* 1 30* 1 03 

102 P HQ =P HO +6*2831 8 
GO TO 100 

103 PH0 = PH0-6* 2831 8 
100 PH I 2= PH I 

PHA= (PHI+PHO )/EM 
PHA=AMOD(PHA* 6. 283 1 3 ) 

I FCPHA- 3* 1 41 59) 105*104*104 

104 PHA=PHA- 6. 2831 8 
PHO = PHO- 6. 2831 8*F9 
GO TO 10 7 

105 I FCPHA+3* 1 41 59) 106*107*107 

106 PHA=PHA+6. 2331 8 
PH0=PH0+6- 2831 8* EM 

107 PHA = PHA* EM 
GOO = PH A- FOB 

I F C GO 0 - 3 • 14159* EM) 25* 30* 30 
30 G00= GOO- 6* 23 3 1 8* EM 

GO TO 3 5 

25 IF < G00 + 3* 1 41 59* BM ) 26* 35* 35 

26 GOO = GO 0 + 6 • 2 8 3 1 8 * EM 

35 I F( GOO-AMQDC 600, 6. 2831 8) ) 36* 39, 33 

36 EMM = EMM + 1 * 3 
GO TO 37 

38 EMP= EMP+ 1.0 

37 F00 = PH A -AMO DC PH A* 6* 28318) 

39 CONTINUE 
T1=CPHA*TGC>+T1 
T2 = T 1 +PHA* T GD 

L I P = MO D ( IMDEX* I 00 ) 

IFCLIP * EO. 0) GO TO 500 
IJ9DEX = INDEX* 1 
GO TO 1 

500 T 3= T24 FO 
PHA=PKA/ EM 

TYPE 50 1 * IMDEX* T3, THETA* PHA* EMP* EMM 

501 F0RMATC1H * I 4* 5 C 2 X * 1 PE 1 0* 3 ) ) 

IMDFX= IMDEX+ 1 

I F < C C A TOO S* AT CO S ) + C AT SIM*ATSIM) ) . GT* 0*6) GO TO 50 
GO TO 1 
50 TYPE 502 

502 F0RMATC1H-, 'HAM DOFF TO TRACK*) 

EM D 


2 04 



Appendix L 


NEW TECHNOLOGY APPENDIX 


After a diligent review of the work performed under this contract, 
no new innovation, discovery, improvement or invention was made. 
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